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Revolutionary 
direct-reading 
altimeter... 

by 
Bulova 


Test pilots flying high-performance craft enthusiastically responded to the new 
Bulova Altimeter. No difficulty was experienced in making accurate, split-second 
readings even when altitude changes exceeded 25,000 fpm. 


Bulova tape presentation ends ambiguous readings...ends pointer scanning... 
increases reading speed and accuracy. These features, coupled with an instrument 
error of less than 10 ft. at sea level, 50 ft. at 40,000, make the Bulova Altimeter 
an important advance in flight safety. 

Experience in precision design, in precision manufacture, is the Bulova tradition, 
the Bulova capability. Has been for over 80 years. For more information write — 


Industrial & Defense Sales, Bulova, 62-10 Woodside Ave., Woodside 77, N.Y. 
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knowledge and basic applications of suc}: 
knowledge. The sciences of astronautics ar 
understood here to embrace selected aspect. 
of jet and rocket propulsion, space fligh: 
mechanics, high-speed aerodynamics, fligh: 
guidance, space communications, atmospheri« 
and outer space physics, materials and struc 
tures, human engineering, overall syster, 
analysis, and possibly certain other scientifi 
areas. The selection of papers to be printe:! 
will be governed by the pertinence of the topi 
to the field of astronautics, by the current o 
probable future significance of the researc! 
and by the importance of distributing the in 
formation to the members of the Society an: 
to the profession at large. 


Information for Authors 


Manuscripts must be as brief as the prop« 
presentation of the ideas will allow. Ex 
clusion of dispensable material and concis¢ 
ness of expression will influence the Editors 
acceptance of a manuscript. In terms o 
standard-size double-spaced typed pages, 
typical maximum length is 22 pages of tex 
(including equations), 1 page of references 

page of abstract and 12 illustrations 
Fewer illustrations permit more text, and vic: 
versa. Greater length will be acceptabl: 
only in exceptional cases. 

Short manuscripts, not more than on 
quarter of the maximum length stated for ful! 
articles, may qualify for publication a 
Technical Notes or Technical Comments 
They may be devoted to new development 
requiring prompt disclosure or to comment 
on previously published papers. Such manu 
scripts are usually published within tw: 
months of the date of receipt. 

Sponsored manuscripts are published 
occasionally as an ARS service to the indus- 
try. A manuscript that does not qualify fo: 
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city, edition and page or chapter numbers. 
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Manuscripts must be accompanied by 
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ment sponsorship. Full responsibility rests 
with the author. 
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Storable and Powertul 
Rocket Fuel Oxidizer 


Chlorine Trifluoride (cir.) 


Much is expected of Chlorine Trifluoride as a rocket fuel 
oxidizer. That’s because it combines good handling and 
storage properties with very high performance. To be 
specific: 


CTF is not difficult to handle! Boiling point of 
Chlorine Trifluoride is 53.15°F. Freezing point is 
—105.4°F. Vapor pressure at 140°F is 80 psia. 


CTF may be stored conveniently! Chlorine Tri- 
fluoride is storable over a very wide range of temperatures. 
You can count on years of storage life! 


CTF has excellent stability! Shock resistance of 
pure CTF is very high. And it is thermally stable to high 
temperatures. 


CTF offers high performance! The high density of 
Chlorine Trifluoride (1.825 at 68°F) leads to outstanding 
density impulse values. An important plus—CTF is hyper- 
golic with hydrogenic fuels over a wide range of pressures 
and temperatures. 


Fluorine and other Fluorine-Based Oxidizers 


CTF is one of several fluorine-based chemicals produced 
by General Chemical which are now considered as excel- 
lent rocket-fuel oxidizers for various missions. Another is 
Bromine Pentafluoride which shares many of CTF’s desir- 
able physical and chemical properties. General Chemical, 
the sole producer of elemental liquid fluorine, has exten- 
sive technical data available on Fluorine and Halogen 
Fluorides. Write today for further information on these 
high-energy oxidizers. 


BAKER & ADAMSON® 
Products 


Typical test stand of Rocketdyne, a division of North American Aviation, Inc. 


GENERAL CHEMICAL DIVISION 


40 Rector Street, New York 6, N. Y. 
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control... 

without 
contact... 


output 


Now—measure pressure, temperature, vibra- 
tion... any physical quantity you can convert 
to a change in capacitance . . . without fuss, 
muss, or bother. Decker’s Delta Unit takes the 
output from the simplest A capacitance sensors 
you can devise, converts it to an analagous 
phase sensitive differential DC voltage for in- 
dication and recording. 

The new 902-1 Delta Unit is based on Deck- 
er’s patented T-42 Ionization Transducer. It 
provides all the advantages of the T-42 ready 
to use in a simple, compact package as versa- 
tile and flexible as your own imagination. 

Accepting initial capacitances from 1 to 50 
uF, the 902-1 provides maximum signals of 
+30 VDC. Sensitivity in the 5-20 uuF initial 
capacitance range is approximately 0.2V/% 
A C. Unit is provided with a single measuring 
probe and cable assembly; however, two meas- 
uring probes may be operated simultaneously. 
Complete information is provided in Instru- 
ment Data Sheet 902-1, available on request. 
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Recent Advances in 
Cryogenic Engineering 


ROBERT B. JACOBS 


National Bureau of Standards 
Boulder, Colo. 


Dr. Jacobs is Chief of the Cryogenic Equipment Section, National Bureau of Standards. He re- 
ceived his B.S. from the University of California and his M.S. from MIT, both in Mechanical Engi- 
neering, and his Ph.D. in Applied Science from Harvard. The author has taught in the Mechanical 
Engineering Department at MIT where he did research in cryogenic engineering. Since 1951, he 
has been at the Cryogenic Engineering Laboratory working on problems of handling, storage and 
transportation of liquefied gases; his work includes research on two-phase fluid flow, pumping and 
cavitation, boiling heat transfer, flow measurement, pressurization and stratification, insulation 


and equipment optimization. 


IMPORTANCE of cryogenic engineering to rocketry 

stems from the desirability of using liquefied gases as pro- 
yellants. This importance will probably increase because 
the use of liquefied gases appears mandatory for the more 
ambitious missions in the not-too-distant future. It is 
well to realize that the large effort put into, and the conse- 
quent rapid growth of, the field of cryogenic engineering is in 
large part due to the importance of rocket propulsion. Be- 
cause of the relatively short time in which the major growth 
of the field has occurred, an exposition of the recent ad- 
vances of cryogenic engineering should cover the whole field 
of eryogenics which is of importance to propulsion. In ad- 
dition, because of security restrictions, proprietary considera- 
tions and a definitely limited interest, a large part of the 
literature in cryogenic engineering has had very limited dis- 
tribution. Finally, there has been a tendency to ignore the 
cryogenic aspects of the problems facing the rocket engineer; 
in some cases the problems resolve themselves, but too fre- 
quently the neglected problem remains to plague the engineer 
with malfunction and failure. 

Perhaps the most important point to stress is that cryogenic 
problems can generally be solved through application of the 
same basic principles which have been successfully employed 
for many years in mechanical, chemical and electrical engi- 
neering. Situations in which new principles (e.g., super con- 
ductivity) arise are very rare. Normally, the significant devia- 
tions of cryogenic problems from the more familiar ones are 
caused by changes in emphasis or differences in order of mag- 
nitude; design difficulties are usually the result of insufficient 
knowledge of properties of materials. For example, if the 
basic phenomena occurring in a centrifugal pump are ade- 
quately understood, prediction of the performance and cavita- 
tion characteristics of a liquid hydrogen pump will be just as 
reliable as the predictions for a cold water pump. Asa matter 
of fact, because of the gaseous and solid contaminants usually 
present in cold water, the hydrogen problem may be easier. 
An engineer who has experience designing valves for operation 
at ordinary temperatures will notice that the relative im- 
portance of some of his problems will change; whereas hy- 
draulic efficiency and leak tightness across the seat and 
through stem packing may have been primary considerations 
in the past, problems concerned with differential contraction, 
freezing of packing materials, icing of movable parts, impact 
strength, seat materials with satisfactory room and low tem- 
perature properties, and insulation may be more important 
in the cryogenic application. Perhaps the most striking ex- 
ample of a specification which is of different order of magni- 
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tude in cryogenic engineering compared with ordinary prac- 
tice is permissible heat leak; insulations with mean effective 
thermal conductivities of less than 0.0002 Btu/hr ft F are 
common in cryogenic equipment, whereas for the same tem- 
perature range one of the best conventional insulations has a 
mean thermal conductivity of about 0.02 Btu/hr ft F. 

The cryogenic problems encountered in rocket propulsion 
can arbitrarily be segregated into those associated with 
ground handling and those associated with equipment and 
processes aboard the vehicle. The differences in these two 
groups are not due solely to weight and space limitations im- 
posed upon the one, but also to requirements pertaining to 
permissible propellant loss, acceptable flow rates, readiness 
conditions, cooldown requirements, ete. (The secondary 
uses of liquefied gases for rocket propelled vehicles, such as 
refrigeration of instrumentation and power for pneumatic con- 
trols, will not be discussed here; they do not involve any 
basic problems not encountered in the propellant system. 
Also, the field of gas liquefaction will not be discussed here.) 

Ground handling equipment is that required to convey the 
propellants from the liquefiers to the final vehicle: Storage 
and transport containers; piping systems including valves, 
couplings, filters, ete.; both low and high capacity, moderate 
head transfer pumps; pressurization equipment, such as tank 
farms, gas generators and compressors; instrumentation for 
measuring, controlling and recording flow rate, temperature, 
pressure, liquid-vapor composition, and liquid level. Ex- 
amples of equipment found aboard the vehicle are: Storage 
tanks; piping, valves and other fittings; high capacity, high 
head pumps; heat exchangers involving boiling as well as 
convective heat transfer; instrumentation. In addition to 
the problems associated with the equipment, there are prob- 
lems associated with procedures and techniques. Should pro- 
pellants be stored at low or high pressure? Should pres- 
surization or pumping be used for a certain transfer operation? 
Should transport containers or piping systems be used? 
Should equipment be evacuated and purged, or just purged? 
Adequate solutions to these problems require information 
concerning, for example, properties of materials, insulation 
techniques, pumping, tank and piping system design, single- 
and two-phase fluid flow, heat transfer and instrumentation 
techniques. 


Properties of Materials 


It is obvious from the preceding discussion that data con- 
cerning all of the properties required for equipment design 
at room temperature and above are also required for the de- 
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sign of equipment to operate at low temperatures. In order 
to make all pertinent information and data available to 
those concerned with cryogenic problems, the Cryogenic En- 
gineering Laboratory, National Bureau of Standards, Boulder, 
Colo., maintains the Cryogenic Data Center. Although it is 
estimated that the Center will not be completely up-to-date 
for another two years, its services are currently available. 

In addition to the mechanical properties commonly used 
in the design of equipment, the design of cryogenic equipment 
requires a particular awareness of thermal expansion co- 
efficients and specific heats. A. knowledge of the former is re- 
quired so that, after equipment is cooled to its operating tem- 
perature, the required clearances are attained and thermal 
stress limitations are obeyed; specific heat data are required 
for prediction of cooldown characteristics of equipment. 
Thermal conductivity is another property which the cryogenic 
designer should bear in mind: If heat leak through metallic 
members is to be minimized, a stainless steel, monel or other 
poorly conducting alloy should be used; if temperature strati- 
fication in liquid containers is desirable, the inner tanks should 
be stainless steel, whereas if stratification is to be inhibited, 
aluminum or copper alloys may be desirable. 

A series of articles by Corruccini (1)! presents data and 
general discussion concerning thermodynamic, transport and 
mechanical properties of solids, including metals, plastics, 
glass and insulating materials, at low temperatures; a 
selected list of useful references pertaining to the properties is 
given in part one of the series. Chelton and Mann (2) have 
assembled a very useful compilation of thermodynamic, 
transport, mechanical, electrical and miscellaneous properties 
of solids, liquids, gases, insulations, etc. Besides (1 and 2) 
and the Data Center, useful collections of cryogenic data ap- 
parently exist only in the private files of those who have been 
working in the field. 

A rather complete bibliography of references on mechanical 
properties of solids below zero deg C has been assembled by 
Kropschot and McColm (3). This report covers the litera- 
ture from 1945 to 1955, and includes an abstract for each 
reference. As with other materials, the strengths of austenitic 
stainless steels, copper alloys and some aluminum alloys in- 
crease greatly with decreasing temperature; however, even at 
liquid hydrogen temperature these specific materials possess 
sufficient ductility to permit their use under essentially the 
same impact and vibration specifications as at normal tem- 
peratures. It follows that equipment fabricated of these ma- 
terials, which is to be operated at all temperatures from room 
temperature down, should be designed using room tem- 
perature strength data; it is also apparent that if weight 
and space specifications are not severe, equipment can be con- 
servatively designed by utilizing room temperature data. If 
weight limitations are severe and equipment is to be sub- 
jected to working stresses only after cooldown, appreciable 
weight saving can be realized by designing to the mechanical 
strength available at the operating temperatures. Because 
of the wide use of welding as a method of joining in cryogenic 
equipment, information concerning the mechanical properties 
of weldments at low temperatures is required by the designer. 
Unfortunately, work in this area is very limited (see Mikesell 
and Reed (4)). 

Because of their light weight, low thermal conductivity and 
adequate mechanical properties, plastics have found wide 
use in cryogenic equipment as insulating supports, gasketing 
material, valve seats, insulation, tank material, ete. How- 
ever, adequate information is not available for optimum de- 
sign in this area; results have been reported by Reed and 
Mikesell (5), McClintock (6) and Corruccini (1). Epoxy 
resins have been successfully used as structural adhesives in 
cryogenic equipment for use at temperatures down to 20 K; 
McClintock and Hiza (7) have reported work in this area. 
Information pertaining to the electrical and magnetic proper- 
ties of solids at low temperatures is meager, and consequently 
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optimal design of electric and magnetic equipment to operate 
at these temperatures is difficult. As most measurements 
are made with the purpose of investigating basic behavior, ma- 
terials which are not of interest to the designer are studied. 
Examples of the type of work which is useful to the engineer 
are afforded by Hall et al. (8), and Brown and Brennand (9). 

Although much information about the properties of fluids 
is still needed, especially in the high pressure liquid 
and the high density gas regions, sufficient information is 
available to produce nonoptimum design. Knowledge of es- 
sentially all of the thermodynamic (pressure, density, tem- 
perature, specific heat, etc.) and transport (viscosity, ther- 
mal conductivity, ete.) properties is necessary for the desired 
fluid flow, heat transfer, pumping, etc., calculations which 
must be made by the rocket engineer. Although the effort is 
apparently not yet intensive enough, programs to determine 
some of the necessary information are now under way. The 
paper of Woolley, Scott and Brickwedde (10) is probably stil! 
the best compilation of the properties of hydrogen. Althougl: 
comparable compendia for the other propellants do not exist, 
it is the intention of the Cryogenic Data Center to issue such 
documents as soon as possible. As there is insufficient space 
here to list even a token sample of the available references, the 
reader is advised to obtain necessary information from (2), 
from the Cryogenic Data Center or directly from the litera- 
ture. 


Insulation Techniques 


The effectiveness of thermal insulations required for cryo- 
genic equipment runs the gamut from no insulation whatso- 
ever (except that afforded by condensates which may form) 
to the best vacuum-type insulations. Perhaps the major 
reason for using insulation is to prevent the vaporization, and 
consequent loss, of liquefied gas. However, there are other 
reasons for using insulation which, for certain applications, 
may be more important; two of these are to prevent cooling 
of adjacent equipment and to minimize pressure drops by 
preventing two-phase flow. 


Uninsulated Equipment 


The greatest, and perhaps only widespread, use of unin- 
sulated ground handling equipment is at rocket engine test 
sites. The storage tanks are usually well-insulated, but the 
transfer systems and run tanks are not; consequently, as much 
as 75 per cent of the liquid oxygen purchased by a site may be 
lost. For a transfer system of any permanence, there seems 
to be no economic justification for not using insulation; in 
many cases, high vacuum insulation can be justified. There 
may be good reasons for leaving run tanks uninsulated; how- 
ever, the one usually given is that insulation would alter the 
dynamic characteristics of the test stand and thus render the 
test results less applicable. Inasmuch as run tanks are not 
flight tanks anyway, the reasoning appears to be faulty. It 
therefore appears that the use of uninsulated equipment at 
test sites could be reduced. However, as the loss of liquefied 
gas due to heat leak is a function of the heat leak per unit 
mass of fluid flowing (11), transfer systems with high flow 
rates can tolerate poor insulation. In fact, high capacity pro- 
pellant loading systems should be uninsulated, since the load- 
ing operation would be complete before an insulation could be 
effective. 

The use of insulation aboard the vehicle, where the weight 
and space limitations are severe, may be determined primarily 
by these limitations. It is probable that in-flight insulation 
of liquid oxygen equipment in present missions is not justified 
except in situations where piping passes through fuel tanks, 
etc. 

The use of uninsulated equipment, including tankage, with 
liquid hydrogen is not such a clear-cut issue. It is probable 
that for very rapid loading, uninsulated transfer systems can 
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be used. The problems of using uninsulated equipment 
with liquid hydrogen are currently being investigated. 
Possibly the most challenging insulation problem in cryo- 
genic engineering is that of developing tankage which will 
carry liquid hydrogen on the more ambitious missions. 

The factors determining the steady-state heat transfer rate 
to liquefied gases in uninsulated equipment depend upon the 
liquid as well as ambient conditions. With liquid oxygen a 
loose frost or snow normally forms. The presence of this 
lxayer can halve the heat transfer rate (12); however, since it 
is easily blown off, it should not be relied upon to provide in- 
sulation. Some investigators have reported that rain bounces 
off uninsulated liquid oxygen containers, while others have 
reported that the rain adheres, building up thick layers of ice. 
It is obvious that controlled investigation should be carried 
cut; some work of this nature is being done. With liquid 
nitrogen, air (as well as water) condenses on the equipment. 
isecause the normal boiling point of liquid nitrogen is above 
the triple point temperature of air, liquid air continuously 
1uns off of the equipment, washing the ice with it. Conse- 
quently, the outside film coefficient, which may be the rate de- 
iermining factor, should be computed from condensing in- 
iormation. Because the normal boiling point of fluorine is 
below the condensing temperature of air, we would expect 
uninsulated fluorine equipment to condense air, and thus the 
leat transfer computations would be similar to those for ni- 
trogen equipment. Experience with uninsulated liquid hy- 
drogen equipment has shown that solid air does not adhere 
to the exposed surfaces to form a layer of solid, but is con- 
tinuously washed off by liquid air. Therefore, heat transfer 
computations should involve film condensation considera- 
tions. Available information indicates that, if the hydrogen is 
boiling, in some cases the condensing heat transfer coefficient 
may be the rate controlling factor, while in others the boiling 
coefficient may be controlling. If the liquefied gas in unin- 
sulated equipment is sufficiently pressurized (or subcooled), 
boiling will not oceur, and although the phenomena occurring 
on the outside surfaces may not change appreciably, the heat 
transfer rate in these cases may be controlled by the film co- 
efficient in the liquefied gas. 


Conventional Insulations 


Whether or not insulation, and especially non-vacuum type 
insulation, is used depends upon the length of time that the 
equipment is in use for each operation. It is obvious that if 
the use time is extremely short, the equipment will not have 
time to cool, much less transfer energy from the environment. 
Also, if complete cooldown of metallic parts does occur, the 
thermal diffusivities of conventional insulations (e.g., fiber 
glass) are usually small enough so that a fraction of an inch 
insulation thickness is adequate. For example, a 1}2-in. 
thickness of fiber glass insulation on a 4-in. liquid oxygen line 
does not reach steady state for almost 2 hr (13). The 
use of conventional insulations, with liquid oxygen transfer 
lines for example, can reduce the steady-state heat link by 96 
per cent (12). 

The primary problem connected with the use of con- 
ventional insulations is the vapor barrier; this problem is 
automatically circumvented with vacuum insulations. The 
objective is, of course, to prevent water, air, etc., from con- 
densing in the insulation, thereby materially decreasing the 
insulating quality; either the condensate can solidify and 
produce an insulating layer with the thermal conductivity 
of the solid, or it can liquefy, flow to warmer parts of the in- 
sulation and evaporate, causing a high rate, regenerative- 
type heat transfer mechanism. There are plastic and glass 
foam-type insulations which are unicellular and can be pre- 
formed; with these it is only necessary to seal along the joints 
with a suitable mastic. Programs to develop satisfactory 
“‘foamed-in-place” plastic insulations are under way; these 
will not have any joints to be sealed; the primary difficulty 
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with these materials is failure due to thermal stresses. With 
the open, or porous, insulations, a covering sheath is generally 
used. This covering can operate in one of three ways: 

1 Itcan be gas tight; because the pressure in the insulation 
space will drop when the equipment is cooled, this type of 
installation must be able to withstand external pressure. 

2 The covering can be gas tight but will be allowed to 
“breathe” gas that does not condense at the operating tem- 
perature; in this way the insulation need not withstand ex- 
ternal pressure. 

3 The covering need not be gas tight, but a sufficient 

quantity of suitable gas is fed into the insulation so that it 
maintains a positive gage pressure therein. 
The second method is most desirable. Besides the difficulty 
experienced in truly sealing a vapor barrier, especially around 
valve bodies, supports, etc., there is the problem of ruggedness 
with respect to both mechanical and environmental abuse. 
The designer should be aware of the possibility that, in 
providing a vapor barrier with the necessary ruggedness, he 
may end with equipment that is as expensive as and heavier 
than vacuum-insulated equipment. 

The thermal design of cryogenic equipment with conven- 
tional insulation involves conventional heat transfer calcula- 
tions. The determination of the transient characteristics 
will require a knowledge of specific heat, density and thermal 
conductivity, whereas a determination of the steady-state 
characteristics requires a knowledge of thermal conductivity 
only (12 to 16). 


Vacuum Insulations 


Although vacuum-type insulations have been used for 
many years to insulate tanks, there has been a strong reaction 
to their use with transfer systems and other components. 
The reason has been the almost universal lack of knowledge 
and technique concerning high vacuum systems; however, 
the situation is changing. Once the suspicion of unre- 
liability had been overcome and enough fabrication experi- 
ence had been obtained so that realistic cost estimates 
could be made, it was possible to demonstrate that, in many 
cases, vacuum insulated systems are economically as well as 
technically desirable (17). With the use of liquid hydrogen 
as a propellant, and the increased handling of helium in liquid 
form (18, 19), the place of vacuum-type insulation in rocketry 
‘annot be denied. 

Vacuum-type insulations can be divided into three classes: 
High-vacuum, evacuated powder, and so-called ‘super- 
insulations”’ or “breakthroughs” which have ‘apparent mean 
thermal conductivities” that are an “order of magnitude” 
lower than the usual designs in the first two classes. The 
details of these insulations are rather well-guarded proprietary 
information. Vacuum-type insulations can be permanently 
sealed (static), continuously pumped (dynamic), periodically 
pumped and sealed, or of the condensing type. In the last, 
the vacuum jacket is filled with a gas (e.g., carbon dioxide) 
which solidifies on the cold surfaces when the equipment is 
put into operation, producing a high vacuum; this system 
has an advantage with respect to long storage life. 

The advantages of evacuated powder over nonshielded high 
vacuum are that, for larger equipment, the heat transfer is 
less and only rough vacua are required; the disadvantages 
are greater insulation thickness, some added fabrication 
problems and more cooldown mass. In the past, high 
vacuum has been made very effective in larger equipment for 
use with hydrogen and helium by employing liquid nitrogen 
cooled radiation shields. It appears that the “superinsula- 
tions’ mentioned above will replace liquid nitrogen shielding 
in many applications. The “superinsulations” of which the 
author is aware are all vacuum-type. It is probable that at 
least some of them are high vacuum, but with clever methods 
of shielding. 

The methods for calculating heat transfer through vacuum- 
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type insulations are straightforward. In high vacuum, only 
radiant energy is significant, and so only information con- 
cerning geometry and optical absorptivities is needed (11, 20). 
If, for some reason, the pressure in a high vacuum insulation 
becomes too high, the energy transport due to gas conduction 
will have to be computed (21). The energy transport 
through evacuated powder occurs by means of radiation, gas 
conduction and solid conduction; the powder absorbs and 
reflects radiation and reduces the gas conduction. Although 
investigators are currently studying the three mechanisms 
separately, engineering design computations are based upon 
mean apparent thermal conductivities (11). As these num- 
bers are derived from the experiments by means of the Fourier 
conduction equation, conventional heat conduction calcula- 
tions are used; the values used for the conductivities must, 
of course, apply to the temperature range in question. Values 
of apparent mean thermal conductivities have been reported 
by Fulk et al. (22); however, more data are needed, especially 
down to liquid hydrogen and helium temperatures. 


Flow Considerations 


The fluid mechanical design of transfer systems for liquefied 
gases does not require the use of new principles. However, 
it will be useful to put down some thoughts and references. 
(Pumping of liquefied gases will be discussed in a later section.) 

Systems which involve single-phase (liquid) flow and sys- 
tems which involve two-phase (liquid and vapor), single- 
component flow are encountered. It is generally desirable 
to have single-phase flow; pressure drops, equipment size 
and costs can be smaller; design computations are simpler and 
more accurate; flow measurement is easier; undesirable 
compressiblity effects, such as choking are reduced (or 
eliminated), etc. Single-phase flow can be achieved by 
assuring that the static pressure is always greater than the 
vapor pressure. A fairly thorough presentation of the single- 
phase problem has been given by Jacobs (11). 

The design procedures for two-phase systems are still 
relatively crude and inaccurate. Although the results may 
have an error of more than 50 per cent, the most widely used 
empirical correlation is that of Lockhart and Martinelli 
(23). This correlation is based upon data obtained for 
isothermal, two-component flow with no mass transfer be- 
tween the phases; it has been shown, however, that the cor- 
relation can be used with single-component systems, with 
appreciable vaporization, if a simple correction is made to 
account for the variation in the axial component of momentum 
(24, 25, 26). Harvey and Foust (27) set up a system of 
equations for one-dimensional flow of two-phase, single- 
component fluids utilizing the correlation of Lockhart and 
Martinelli to account for the frictional pressure drop. How- 
ever, since the labor connected with numerical computations 
based upon these equations is very great and the errors in- 
herent in the equations (due to shortcomings of the Lockhart 
and Martinelli correlation) may be large, design calcula- 
tions based upon the simple momentum correction (26) may 
be preferred. An insight into the complexity of the method 
of Harvey and Foust can be obtained from the work of Rogers 
(28), who made some calculations for the two-phase flow of 
hydrogen; Rogers’ results have not been compared with 
experimental results. 

The extreme complexity of the two-phase flow problem has 
prevented appreciable progress through more _ theoretical 
approaches; although the work of Linning (56) is less em- 
pirical than that in the preceding references, numerical 
calculation depends upon an appreciable amount of experi- 
mental information. There are at least two groups carrying 
out theoretical work on annular-type, two-phase flow, but 
results have not yet been published. 

The flow problem about which least is quantitatively known 
is that concerned with the transient phenomena that occur 
when cryogenic equipment is cooled from atmospheric tem- 
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perature down to operating conditions. The two questions 
which arise are: How fast can systems be cooled, and what 
pressures occur during the various possible cooldown pro- 
cedures. In many cases, sufficient time is available and the 
equipment can be sufficiently well designed so that cooldown 
offers no difficulty. However, if one is not aware of potential! 
pitfalls, he can lay out a transfer system that could never be 
cooled. 

In their most detailed aspects, cooldown problems involve 
consideration of unsteady flows of highly compressible fluids. 
two-phase flows, convective and boiling heat transfer, etc 
Because of the extremely low critical velocities of two-phase 
systems, the flow is readily choked; because of the rapic 
vaporization which occurs when liquid contacts a relativel) 
warm wall, there are large pressure surges which can revers« 
the direction of liquid flow. A common technique for decreas- 
ing the cooldown time of relatively long transfer systems is 
to provide vent valves through which the large amounts of{ 
vapor can be bypassed, thereby eliminating the necessit; 
of having to push large volumes through the whole system. 

Apparently all of the investigations, to date, of the cool- 
down of eryogenic systems have been oriented toward the 
prediction of cooldown time. Only the work of Burke et al. 
(29) appears in the open literature; they report reasonable 
correlation between theoretical predictions and experimenta! 
results for transfer systems in which the flow is controlled by 
a choking orifice near the exit; the experiments were carried 
out with liquid oxygen. Macinko (30) has developed a 
method for the prediction of the cooldown time for transfe1 
systems without restrictions and has obtained good agreement 
between his predicted cooldown times and those actually 
measured with both hydrogen and nitrogen. 


Pumping Equipment 


There is a wide range of specifications for pumping equip- 
ment which may be used in connection with cryogenic pro- 
pellants. Besides the obvious differences in weight and bulk 
between ground handling and spaceborne equipment, head 
specifications can cover the whole range up to 50,000 ft and 
capacity specifications can cover the whole range to more than 
10,000 gpm; high shaft speeds will be necessary. It is ap- 
parent that positive displacement, centrifugal and axial 
pumps have their places. Also, with hydrogen and helium, 
compressibility effects are appreciable. It should be stressed 
that the fundamental concepts and methods, which are 
adequate for the fluid-mechanical design of pumps for common 
fluids, are just as applicable to liquefied gases. The per- 
formance characteristics will differ because of such things as 
greater leakage due to lower viscosity; the efficiency may be 
lower due to the fact that a smaller density will require less 
total power input, but the fixed losses (such as bearing) will 
not decrease appreciably. The cavitation characteristics will 
differ because of differences in thermodynamic and transport 
properties, but the same fundamental considerations apply 
in all cases. 

The centrifugal-type pump has been most widely used with 
liquefied gases although turbine-type transfer pumps are com- 
mon. Reciprocating pumps have long been used to compress 
liquefied gases for the production of high pressure gases; of 
course, this type of equipment has low capacity, but pressures 
up to 10,000 psia are achieved, and development of equip- 
ment for pressures to over 13,000 psia is under way (81). 
Some work has been done on vane-type cryogenic pumps, but 
no results have been published. Work is also progressing on 
airborne multicylinder pumps (32) for providing high pressure 
gas for purging, pneumatic control, cooling, etc. (33, 34). 

The work done on the pumping of liquefied gases with 
centrifugal pumps indicates that, as far as performance 
characteristics are concerned, no extraordinary behavior is to 
be expected. The years of experience with liquid oxygen 
ground handling and rocket engine pumps, as well as the work 
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of Jacobs, Martin et al. (35) with water, oxygen, nitrogen and 
hydrogen, and the work of Caine et al. (36) with hydrogen, 
substantiate this. However, studies into cavitation char- 
acteristics are still needed. Substantial progress has been 
made in this area recently by analyzing cavitation char- 
acteristics through consideration of the thermodynamic 
properties of the fluids (35, 36, 37, 38). For example, it 
was predicted theoretically and shown experimentally (35) 
that the minimum required net positive suction head (NPSH) 
for hydrogen is only about 75 that required for nitrogen in 
the same pump; in a water pump with a specific speed of 
1930 (at 3450 rpm) only 2} in. NPSH were required to 
prevent cavitation with hydrogen (35), whereas a booster- 
type pump that was designed for hydrogen and had a specific 
speed of about 800 (at 12,000 rpm) required an NPSH of 
only 6 in. (36). The measurement of these NPSH’s with 
»dequate accuracy requires vapor-pressure instrumentation 
39). 

Besides the low NPSH requirement, hydrogen pumps are 
different from others because the low density requires very 
high heads for high pressures. Where weight and space limi- 
tations are severe, it is necessary to use high speeds. One 

levelopment along this line is the radial-vane pump described 
»y Bosco and Carter (40). 

Two problem associated with pumps and which have 
1ot been mentioned previously concern bearings and seals. 
[he bearing problem is frequently solved by using shafts 
which are so long that frictional energy dissipation, heat 
transfer from the atmosphere, or heaters can keep the bear- 
ings warm; in these cases, no exceptional bearing problems 
rise, but the disadvantages of long shafts and the possibility 
of appreciable heat transfer into the pump are present. 

The advantages of submerging the bearings in the liquid 
which is being pumped have stimulated work along this line. 
This can be achieved by either using a sump configuration or 
by bleeding some of the liquid from the pump discharge 
through the bearings. All of the work on submerged bearings 
reported to date has been with thoroughly cleansed and 
degreased ball bearings; this work indicates that the 
separator is the critical item in the performance of ball 
bearings submerged in liquefied gases. Bearings have been 
evaluated in pumps operating from 3000 to 22,000 rpm in 
hydrogen and nitrogen and in a test rig from 3000 to 10,000 
rpm in nitrogen (41). The most successful separator mate- 
rial has been fabric-reinforced phenolic. Satisfactory per- 
formance has been achieved with angular contact bearings 
with balls and races made of both 52100 steel and 400-C 
stainless; the stainless is preterable, however, because of its 
greater resistance to corrosion in the event that water con- 
denses on the bearings. Running times of hundreds of hours 
have been achieved in tests up to 12,300 rpm. Long lives 
may also be attainable to higher speeds, but tests have not 
yet been run; the tests at 22,000 rpm in hydrogen were ter- 
minated after about 6 hr, but not because of bearing failure. 

Small ball bearings have been run successfully in gaseous 
hydrogen from room temperature down to liquid hydrogen 
temperature at 10,000 rpm (42). These bearings had 440 
stainless steel balls and races, used micarta ribbon retainers, 
and ran successfully both with and without oil. 

As in the case of bearings, the problems arising from the 
operation of shaft seals at low temperatures have been cir- 
cumvented by putting the seals in ambient temperature 
regions. The author is not aware of any literature on low 
temperature liquid seals except for the limited and unde- 
tailed comments in (35 and 40). One reason low temper- 
ature shaft seal development has lagged is because sump-type 
pump configurations do not require cold shaft seals; if the 
pump is driven by an extension shaft with the drive in the 
atmosphere, only a seal to hold the sump pressure is needed; 
if the driver is in the sump as an integral part of the pumping 
unit, no seals whatsover are needed. It has been thoroughly 
demonstrated that electric motors can be designed for opera- 
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tion submerged in any liquefied gas (36). There is no ap- 
parent reason why submerged turbine drives cannot be 
developed for the larger power requirements. 


Heat Transfer 


Although heat transfer considerations in cryogenic equip- 
ment do not involve new principles, it should be pointed out 
that there is a definite lack of heat transfer data and that 
some of the published results are contradictory. Many 
convective heat transfer calculations for cryogenic equip- 
ment are based on the accepted correlations obtained from 
data with more common fluids. Of course, when new geome- 
tries are used, the designer should be most skeptical in ex- 
trapolating data. Much heat exchange in cryogenic equip- 
ment involves boiling and/or condensing. Most of the in- 
formation available has been obtained with simple geometries 
(48, 44, 45) (wires, cylinders, etc.) or with special geometries 
(46), and so is not directly applicable to most heat exchanger 
designs. Although investigations concerned with heat trans- 
fer with liquefied gases flowing through actual heat exchanger 
configurations are in progress, no results have appeared in 
the open literature. 


Instrumentation 


The instrumentation with which the cryogenic engineer is 
concerned and which is important in the evaluation and opera- 
tion of propulsion systems is related to pressure, temperature, 
flow, liquid level and static mass measurement. As is gener- 
ally the case with other equipment, instrumentation in 
eryogenics is not significantly different from conventional 
instrumentation except that the low temperature, and pos- 
sibly high vacuum, requirements may effect the details of 
some of the sensing elements. 


Pressure 


Static (i.e., steady-state) pressure measurement offers no 
difficulty in cryogenic equipment. The geometry of the pres- 
sure tubing should be such that liquid cannot continuously 
drain or percolate into it, and if heat conduction is to be 
minimized, the tubing should be stainless steel. All of the 
normal procedures, such as the use of piezometer rings and 
prevention of unknown hydrostatic heads should be followed. 
Of course, all contaminants which can freeze and block the 
pressure taps and connecting tubing must be removed before 
cooldown. 

Dynamic (i.e., unsteady) pressure measurement offers no 
difficulty if the transducer can be located so that its temper- 
ature isknown. This can usually be accomplished by having 
the transducer in the atmosphere and connecting it to the 
pressure tap with a piece of stainless steel tubing; the only 
reason it may be inadvisable to do this is that the dynamic 
response of the system, with the connecting tubing, may be 
unsatisfactory. There are transducers on the market which 
are made and calibrated for use down to about —350 F and 
possibly lower. As the zero shift with temperature of these 
units is quite large, the temperature of the transducer must 
be known. Thus transducers can be connected directly to 
equipment if the pressure data are needed only when the 
transducer temperature is known. However, if pressure 
transients are to be measured during the cooldown of equip- 
ment, a transducer connected directly to the equipment will 
not yield accurate pressure information directly. The signals 
or records will have to be manipulated to account for an ef- 
fective transducer temperature at each instant. 


Temperature 


The problems of temperature measurement can also be 
divided into those associated with static and those associated 
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with dynamic situations. The former can be adequately 
handled, but for the latter better response characteristics 
are desired. For a comprehensive exposition of the field of 
thermometry, the reader is referred to (47 and 48). 

Five types of thermometers are commonly used in cryogenic 
applications: Thermocouples, vapor-pressure bulbs, gas 
thermometers, platinum resistance thermometers and semi- 
conductor thermometers. The choice of system to be used 
depends upon the conventional considerations, such as ac- 
curacy and sensitivity, cost and dynamic characteristics. 

Thermocouples have wide use in cryogenic engineering, 
their main assets being adequate accuracy, good response 
time and moderate cost. A primary disadvantage seems to 
to be a drift in calibration (or instability) due to thermal 
cycling. Of course, as is true with all thermometers, the 
installations must be such that the junctions are adequately 
tempered; that is, heat conduction through the leads and well 
to the junction must be made negligible. As the thermo- 
electric power of thermocouples decreases at low temperatures 
(for copper-constantan it is only 6 uv per deg K at 20 K), 
junctions with higher powers are continually being sought; 
a junction which has found wide use, and which has a power 
of 15 uv per deg K at 20 K, is gold with a 2.1 atomic per cent 
cobalt impurity vs. copper (49). Because the instability 
is believed to be due to the copper, work is now being done 
on a junction of the gold-cobalt, mentioned above, vs. silver 
with an 0.37 per cent gold impurity. 

Vapor bulbs have the advantages of extremely high ac- 
curacy and simplicity, but have limited ranges. Tests in- 
dicate that vapor bulb systems can have time responses that 
are about as fast as those of thermocouples made with 30 gage 
wire. In equipment where the measured temperature range 
is between the triple and critical points of a suitable substance, 
the use of vapor pressure thermometers should be seriously 
considered. 

Gas thermometers have the serious disadvantages of being 
difficult to fabricate and difficult to use. They are therefore 
probably undesirable for use in engineering equipment. 

Platinum resistance thermometers have the advantages of 
stability of calibration, accuracy, linearity (down to 40 K) and 
a large temperature range. Their primary disadvantages 
are cost and procurement time; the author is not aware of any 
data concerning the response characteristics of platinum 
resistance thermometers, but they are relatively sluggish. 

Semiconductor thermometers are resistance thermometers, 
employing semiconductors, such as carbon, germanium, 
silicon or metallic oxides as the sensing elements. Their 
advantages include relative simplicity, accuracy and, in some 
cases (50), stability when subjected to thermal cycling. They 
have the disadvantage of nonlinear calibration curves, and 
in some cases (e.g., carbon) they are unstable. The large 
changes in the resistance of the units which occur over ap- 
preciable temperature ranges require that special considera- 
tion be given to the circuitry if the measurements are to be 
optimized (51). Because these sensors can be made rather 
small, their response characteristics can be fast. 


Flow 


One of the major instrumentation problems in the cryo- 
genic aspects of rocketry is flow measurement. Once again, 
it can be said that the basic considerations with liquefied 
gases do not differ from those with the more common fluids. 
Some development work of mass flow meters for use with 
liquefied gases is now in progress, but essentially all flow 
measurement has been volumetric. Turbine-type meters 
have been most extensively used; however, their superiority 
over other types of meters is not apparent. One of their 
major difficulties is that if two-phase flow inadvertently 
occurs, the turbines can be wrecked. Based upon work with 
hydrogen, nitrogen and water, Richards (52) concludes that 
sharp edged orifices are reliable flow measuring devices for use 
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with liquefied gases; the results were reliable except when two- 
phase flow was known to enter the orifice. 

The primary difficulty in flow measurement is encountered 
if two-phase flow exists. This problem, which is common to 
all volumetric flow meters, will disappear if a mass flow meter 
is developed, but until then, the development of an instrument 
to determine the liquid-vapor composition of a flowing fluid 
is needed. There are several programs under way which 
aim to perfect equipment for detecting the presence of vapor 
in liquid (53), but the adequacy of the instruments so far 
developed for the quantitative measurement of the compos '- 
tion of two-phase flows is questionable. 


Tank Contents 


Another type of instrumentation important in handliny 
liquefied gases is that required for determining the amount «/ 
liquid in a tank. Several adequate systems are available for 
use with ground handling equipment, but attainment of th: 
desired accuracy in uninsulated tanks still presents problem: 
The instruments can be classified as direct weighing, hydro- 
static and liquid level. The direct weighing of cryogeni: 
tanks is not different from other tanks; load cells using straii: 
gages, inductance, capacitance or resistance principles, o 
some other sensing device can be used along with conventiona! 
circuitry. 

The use of hydrostatic pressure systems is widespread 
depending upon what the pressure signal is to be used for 
either electrical or mechanical sensors, recorders, etc., may be 
used. The location of the pressure taps and their layout is 
important. In many installations the taps have been placed 
in the vent and liquid lines; consequently, when there is liquid 
or vapor flow, the instrumentation cannot be used. The sys- 
tem must be such that the location of the liquid-vapor inter- 
faces in the pressure lines is known. In some systems, the 
high pressure tap has been taken directly from the bottom of 
the tank, permitting percolation; this produces a disturbance 
in the signal. Ina relatively well-insulated tank, the mass of 
liquid can be obtained from a hydrostatic system if the tank 
geometry is known. However, if there are appreciable con- 
vection currents, there may be appreciable error in a mass 
determination. 

Liquid level devices are widely used for control as well as 
indication of liquid level. As these are simply height measur- 
ing devices, information concerning tank geometry as well as 
fluid density must be known if fluid volumes and masses are 
to be determined. The capacitance liquid level device (54) 
has been used successfully for many years for the indication, 
control and recording of the liquid level of liquefied gases; 
it is a continuously indicating device and has a fast time re- 
sponse. For point (and multipoint) indication, control and 
recording, the carbon resistor liquid level device has been 
successfully used for years. Because the properties of carbon 
resistors drift under thermal cycling, the associated circuity 
should be made relatively insensitive to this. With proper 
circuit design, the carbon-resistor level indicator can respond 
in 2 to 3 sec. Some work has been done on ultrasonic-type 
indicators, vibrating paddle-type indicators and radiation- 
type indicators, but no information is available in the open 
literature. 


Summary 


The point to be emphasized is that cryogenic engineering 
is a synthesis of the older branches of engineering (mechanical, 
chemical, electrical), and therefore does not involve new 
basic principles. The discussion presented here has avoided 
those areas where usual engineering experience is directly 
applicable. An attempt has been made to discuss situations 
where the unwary may go astray and where the cryogenic 
engineer has evolved new techniques. In addition, an at- 
tempt has been made to indicate the location of data and 
information applicable to new materials and situations, as 
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well as that applicable to the familiar materials and situations 
in the low temperature regions. Finally, it should be pointed 
out that a comprehensive treatise on cryogenic engineering 
(55) is to be published early in 1959. 
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Atomization in High Velocity 
Airstreams 


MALCOLM A. WEISS! and 
CHARLES H. WORSHAM? 


Esso Research and Engineering Co. 
Linden, N. J. 


An experimental study was made of the drop sizes obtained on injecting a liquid into large hoi 


airstreams of sustained high velocity. 


The liquid, a molten synthetic wax, was injected contra- 


stream or costream through simple cylindrical tubes. Downstream, a traversing probe withdrew 
a representative sample of the stream, cooled it and froze the droplets. The collected solid par- 
ticles were analyzed by sedimentation and by sieving. The results were correlated empirically by 


the dimensionless equation 


2 
OL 


OL PL 


Mass median diameter (X), air density (94), relative velocity (V), liquid viscosity (u,) and mass 
injection rate (W) were changed over 4- to 25-fold ranges. Surface tension (¢,), liquid density 
(o,) and air viscosity (u4) were not changed significantly. 


N EXPERIMENTAL study was made of the drop size 
distributions resulting from spraying liquids into large 
airstreams of sustained high velocity. The intent of this 
study was to simulate fuel injection in turbojet afterburners 
and in ramjets. Although there are much data on drop size 
for liquids atomized by small, rapidly expanding airjets, 
there are very little data on atomization in large, steady flow 
airstreams of high velocity. The latter is the physical sys- 
tem in jet engines, and it is different enough from the systems 
studied to warrant investigation. 

A survey of possible techniques for measuring drop size 
distributions led to the choice of a modified ‘“substitute- 
fluid” method for this work. The substitute-fluid method 
injects a liquid which is solid under normal conditions. How- 
ever, when molten under injection conditions, the liquid has 
physical properties similar to those of the real liquid fuels of 
interest. The molten drops can be frozen after atomization, 
collected as solid particles and analyzed for drop size 
distribution. 

In order to obtain meaningful results, Longwell (1)? (and 
preliminary experiments here) showed that there must be 
no marked change in temperature until atomization is com- 
pleted. Ifa hot molten material is sprayed into cold air, liquid 
temperature can drop, and freezing can begin while drop 
breakup is still occurring. At the least, rapid cooling of the 
liquid leads to an ambiguity in temperature (and thus in 
liquid physical properties) during atomization. Therefore, in 
all of this work the molten liquid was sprayed into hot air at 
the same temperature as the liquid, and the drops were not 
cooled and frozen until drop breakup was completed. 


Experimental Methods 


The liquid injected in all tests was molten Acrawax C, a 
synthetic wax manufactured by the Glyco Products Co. It 
melted at 284 to 290 F and was not volatile; the vapor pres- 
sure was estimated to be 10-> mm Hg at 300 F. Added to 
the wax was 4 per cent of alcohol-soluble nigrosine black dye; 
the dye made the solid drops opaque for photomicrographic 
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calibrations. The dyed wax had the following physical 
properties 


300 F 350F 400 F 


specific gravity 0.828 0.811 0.806 
surface tension (dyne/em) 22.0 20.1 18.2 
viscosity (centipoise) 11.3 5.86 3.25 


Except for the study of liquid viscosity effects, all runs were 
made at 300 F. 

Dyed wax was melted in a heated blowcase containing a 
large spherical float (to minimize diffusion of the pressuring 
inert gas into the liquid). The liquid was piped through 
heated lines to a turbine flowmeter, then through a high 
surface heat exchanger (to match air and liquid temperatures), 
and finally to the injector. 

The injectors themselves consisted of cylindrical tubes, 
¢i-, gs- and ;%;-in. ID. Each tube was supported and fed 
by a thin strut of flattened tubing, welded in a ring of 6-in. 
pipe (see Fig. 1). Injector axis always coincided with the 
air duct axis, but the orientation could be either costream 
or contrastream. 

The air duct consisted of 6-in. pipe. Two perforated plates 
and a straight run of about 14 ft preceded the flanges between 
which the injector ring was clamped. The length of pipe 
downstream of the injector could be varied as the sampling 
distance required. The air compressor available allowed 
velocities up to about 1000 ft per sec in the duct at atmos- 
pheric pressure; at higher pressures, maximum velocities 
were proportionately lower. 

For tests at atmospheric pressure, the length of duct down- 
stream of the injector simply exhausted to the atmosphere. 
The entrance to the sampling probe then was located in the 
exit plane of the duct, and samples could be taken at fixed 
points or while the probe was moving in a traverse path 
across the exit plane. For tests at pressures above atmos- 
pheric, the duct exhausted into a 12-in. pipe section with a 
variable exit area. The sampling probe was contained 
within this section, and the sampling probe entrance was 
again located in the duct exit plane. However, the probe 
could not be moved during high pressure tests, and all samples 
were taken at a single point 0.5-in. above the duct axis. 
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The sampling probe used is shown in Fig. 2. This design 
was developed in order to: 

1 Sample at stream velocity to prevent size discrimination 
while collecting. 

2 Collect and freeze droplets without distortion. 

3 Collect a reasonable amount of sample in a reasonable 
length of time. 

4 Enable measurements at fixed locations or across the 
traverse path. 

In order to traverse, the entire probe was mounted in a 
counterbalanced assembly with a follower which tracked 
in a slotted template. The path described by the probe 
axis is shown in Fig. 3. During normal operation, air was 
withdrawn through the probe filter at a rate which equal- 
ized static pressure in the probe tip and static pressure 
in the airstream; this procedure is the most convenient 
way to obtain sampling at approximately stream velocity. 
The solid particles collected in the probe were flushed out 
after a run, dried and analyzed. 


Evaluation of Analysis and Sampling 


The solid spheres removed from the probe were analyzed 
for drop size distribution in the Micromerograph—a com- 
mercial sedimentation-in-air analyzer. Two sets of five 
replicate runs were made to determine reproducibility of 
the analysis; one set used a sample with a mass median 
diameter of 22 yu, and the other set used a median of 87 uy. 
For both sets, the rms error was less than 5 per cent of 
diameter, at a given volume percentage, for cumulative 
volumes between 5 and 95 per cent. (Particles larger than 
150 w were always separated from the Micromerograph feed 
and were analyzed by sieving.) 

Absolute calibrations of the analyses were made by com- 
paring the results with manual sizing and counting of images 
from enlarged photomicrographs of 11 different samples 
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Fig. 2. Droplet sampling probe 
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Fig. 3 Path of sampling probe traverse at duct exit 


ranging in mass median diameter from 13 to 95 uw. From 
1700 to 11,000 images were counted for each sample. At 
the mass median diameter (50 per cent cumulative volume), 
there was an rms difference in diameter between the two 
types of analyses of about 8 per cent. Comparisons were 
also made at 15.9 and 84.1 per cent cumulative volumes; 
these percentages represent variations of one geometric 
standard deviation from the mass median for a log-probability 
drop size distribution. At the 84 per cent level, the rms dif- 
ference in diameters was about 6 per cent; but at the 16 per 
cent level, it was about 24 per cent with the Micromerograph 
usually showing larger particles. It was interesting that this 
error—in sizing the smallest part of a sample—was not much 
affected by the absolute level of diameter. A given particle 
diameter was measured more accurately if it occurred at the 
50 per cent point rather than at the 16 per cent point. 

Attributing all of the error to the sedimentation analysis 
probably exaggerates the defects of that analysis. There 
were also variations from one photomicrograph to another 
of a given sample. Therefore, because of the speed and sim- 
plicity of the Micromerograph (and sieve) analyses, they 
were used throughout for the data reported here. Unless 
the smallest drops in a given spray are of critical importance, 
as they might be in some ignition phenomena for example, 
the overall accuracy of the analyses is considered adequate 
for any likely application of the data. 

Development of analytical techniques made possible an 
assessment of the sampling probe. First, particles collected 
from the probe were found to be almost uniformly spherical. 
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There was no consistent distortion of the largest drops, 
which are most difficult to freeze before impaction. On the 
average, less than 5 per cent of the total sample weight 
collected was found as a film on the probe walls; a film would 
prove impact before freezing. 

Air withdrawal rate through the probe was deliberately 
varied (unbalancing the probe nose and stream static pres- 
sures) in a series of runs. Only when the mismatch between 
probe nose and main stream velocities reached about 100 ft 
per sec was any significant shift observed in the size distribu- 
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tion of particles collected. Since a match within 10 to 30 ft 
per sec was easily obtained, sampling velocity was not 
critical. 

The presence of the probe did not greatly distort the ve- 
locity profile upstream of the probe. A small (j-in. long) 
Kiel tube was mounted on the duct axis, and the sampling 
probe was moved across a diameter in a plane just downstream 
of the Kiel tube. For the complete sample probe trave! 
(from directly behind the Kiel tube to completely clearing 
it), the maximum change in air velocity at the axis was less 
than 1 per cent. 

There was no preference for any particular particle siz 
during sampling. This was demonstrated by injecting :. 
“spray”’ of solid spherical particles and collecting a sampl 
in the normal way used for liquid drops. Three analyse: 
of the particles fed are compared in Fig. 4 with analyse: 
from two runs in which the solid particles were collected 
Within analytical error, the results were identical. 

Finally, the adequacy of the traverse path chosen wa: 
assessed. An error in average particle size for the tota 
spray would exist only if both (a) local particle size distribu- 
tion varied appreciably from point to point in the spray, anc 
(b) the sampling probe did not collect for an approximately) 
equal time at all equivalent points in the cross section. 

Local samples were taken from three sprays, each result- 
ing from a ;%-in. tube injecting contrastream at 1 gpm. 
Two sampling distances were tested—23 and 41 in. betweer 
injector tip and sampling plane—and two air velocities— 
350 and 750 ft per sec. The results are plotted in Fig. 5. 
They show that in none of the three sprays did any loca! 
median diameter vary by more than 10 per cent from the 
average between maximum and minimum local median 
diameters. 

The total mass collected by the probe during spray traverses 
was 100 per cent of the theoretical collection with a standard 
deviation of about 20 per cent. Considering that local liquid 
mass velocities varied by as much as 100 to 1 in a given 
spray, this indicates that the traverse used did give a rea- 
sonably weighted sample of the entire spray. 


Experimental Results 


Using the experimental methods outlined above, a study 
was made of the sprays from simple cylindrical tubes. The 
variables studied and the approximate ranges covered were 
as follows: Relative air velocity (V, 200 to 1000 ft per sec), 
liquid injection velocity (v, 4 to 100 ft per sec), injector in- 
side diameter (D, 4% to 7 in.), liquid viscosity (uz, 3.2 to 
11.3 centipoise), air static pressure (P, 1 to 5 atm—corre- 
sponding to an overall air density, p4, change of 0.046 to 0.26 
Ib per ft*) and both costream and contrastream injection. 
The liquid injection rate varied from 0.1 to3gpm. Resulting 
mass median diameters (X) ranged from 19 to 118 yu. The 
experimental data, summarized in the following sections, are 
presented in detail in (5). 


Effect of Relative Air Velocity 


Velocity was the variable with the largest effect, by far, 
on drop diameter. The most useful velocity to consider was 
relative velocity between liquid- and airstreams (i.e., air 
velocity plus liquid injection velocity for contrastream injec- 
tion or minus liquid injection velocity for costream injec- 
tion). At a given relative velocity—with other parameters 
also constant—the same drop sizes were obtained when in- 
jecting either costream or contrastream. 

Most of the data obtained are shown in Fig. 6; the results 
are separated by injector diameter and direction of injection; 
liquid injection velocity is the parameter. If all of the slopes 
of the sixteen lines of Fig. 6 are considered, the average ve- 
locity exponent is 1.33 with a standard deviation of 0.13. 
Two of the lines have markedly low slopes, both resulting 
from the @; tubes at low liquid injection velocity (and both 
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Fig. 6 Effect of relative air velocity on mass median drop 
diameter 


evidently caused by unexpectedly small drops at low relative 
velocity). Excluding these two slopes, the average exponent 
becomes 1.37 with a standard deviation of 0.10. Therefore, 
the difference in precision offered little justification for elimi- 
nating these two slopes, and an average velocity exponent 
was taken as 1.33, i.e., X ~ V~33, This velocity exponent 
compares to the value of 1.0 found by Nukiyama and Tana- 
sawa (4) for small pneumatic atomizers, and to the value 
of 1.11 for a metal alloy and 1.68 for a wax found by Marshall 
(2) for large venturi atomizers. 

There was no consistent change of relative velocity ex- 
ponent with injector diameter, direction of injection or liquid 
injection velocity. 


Effect of Liquid Injection Velocity 


Fig. 6 makes it clear that liquid injection velocity per se 
affects drop size in addition to its contribution to relative 
velocity between liquid and air. This effect is shown in 
Fig. 7 for each injector diameter, combining the costream 
and contrastream injection points. The data are grouped 
at the two levels of relative velocity, 350 and 750 ft per sec, 
at which most of the results were obtained. When a run 
was not made at exactly these relative velocity levels, the 
diameter was corrected by the relative velocity exponent read 
from the appropriate individual line in Fig. 6. 

As liquid injection velocity increases, there is an increase 
in drop diameter under all conditions. The effect is small 
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Fig. 7 Effect of liquid injection velocity on mass median drop 
diameter. Open points are contrastream, solid points costream 


except in one case; again, the ,°; tubes produce drops smaller 
than expected at low relative velocities and low liquid in- 
jection velocities. Excluding this case, the average liquid 
velocity exponent is 0.07 with a standard deviation of 0.02. 
Including it, the average exponent is 0.09 with a standard 
deviation of 0.07. If less weight is given to the single high 
value, a reasonable compromise exponent would be 0.08. 


Effect of Injector Diameter 


The curves of Fig. 7 are crossplotted in Fig. 8 to establish 
the effect of injector diameter. Thus, different injectors 
can be compared at constant relative velocity and at constant 
liquid injection velocity. 

At 750 ft per sec, the results are quite regular; drop size 
increases slowly and consistently as injector diameter in- 
creases at all liquid injection velocities. The line slopes 
range from 0.12 to 0.18, averaging 0.15. However, at 350 
ft per sec, the results for the smallest tube are again spread— 
as would be expected from the high liquid velocity exponent 
observed for this tube in Fig. 7. The lines in Fig. 8 connect 
the data points directly. If an attempt is made to draw 
“best straight lines” through all corresponding points, the 
slopes range from zero to about 0.35, averaging about 0.19. 
Giving a slight weight to this higher average value, a com- 
promise exponent for diameter would be about 0.16. 

The unusual behavior of the ;°; tube at low relative velocity 
and low liquid injection velocity is not unique to one par- 
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ticular geometry of injector. Both the “long” and “short” 
gz tubes shown in Fig. 1 gave about the same results (see 
Fig. 7). Therefore, it is likely that a different mechanism of 
atomization begins to become important for these small 
tubes at low liquid and air velocities. Drops produced by 
such tubes will not be accurately predicted by results from 
the same tubes under other conditions, or from other tubes 
under similar conditions. 
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Fig.8 Effect of injector diameter on mass median drop diameter 
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Effect of Sampling Distance 


The results discussed up to this point were obtained samp- 
ling at distances of 33 to 41 in. from the injector tip. Several 
runs were made to determine the effect of sampling closer to 
the injector. All of these tests used the ;%;-in. tube injecting 
contrastream at a liquid velocity of 12 ft per sec. 

Fig. 9 shows that a reduction in sampling distance from 
41 to 23 in. had no observable effect. In one run there was 
a further reduction, to 11 in. The resulting median diameter 
was slightly higher than the mean of the other points, but th» 
difference, based on this single run, is not regarded as sig- 
nificant. Therefore, one can conclude that no important 
changes in drop size distribution occur at distances great: 
than about | ft for conditions similar to those tested. 


Effect of Air Density 


The effect of air density was studied by varying stati: 
pressure of the air stream while sampling at a fixed location 
All of the tests were run with wax injected contrastream a 
about 4 gpm (16 ft per sec) through the 4%5-in. tube. Stati 
pressure was varied from 1 to 5 atm abs with maximum ai 
velocity limited by compressor capacity. The minimum ai 
velocity was about 200 ft per sec; lower velocities resulte: 
in gravitational classification of particles between injectio: 
and sampling points. 

In the first tests, results at atmospheric pressure for poin' 
sampling in the “pressure duct’’ were compared to results at 
atmospheric pressure in the normal duct with full sampling 
traverses. ‘The data, in Fig. 10, show no evident differences. 

All of the pressure data obtained are plotted in Fig. 11 (in 
two parts to minimize confusion of data points). For each 
set of points at a given pressure level, a best straight line was 
drawn and the slope measured. No slope but one differed 
by more than 6 per cent from the average of 1.33 determined 
for all data at atmospheric pressure. The single exception, 
at 3 atm, could be shifted greatly by change of one point. 
Therefore, the effect of relative velocity on drop diameter 
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eems to be unaffected by pressure in the pressure range 
covered. 

The effect of pressure at given velocity levels was deter- 
ained by erossplotting points from the lines of Fig. 11. 
ig. 12 shows that a fairly uniform pressure effect was found 
when diameter was plotted against (1 + 1/P). The slopes 
of the lines in Fig. 12 vary from about 1.0 to 1.2, averaging 

Because of the scatter of original data points in Fig. 11, 
the exponent of 1.1 was simplified to 1.0. Original pressure 
duct data were then correlated by the group (1 + 1/P)/V**. 
The correlation showed that the equation X = K(1 + P)/V‘* 
‘with K a constant) described the 25 pressure duct median 
diameters with a standard deviation of 6 per cent of diameter. 


Effect of Liquid Viscosity 


Seven tests were conducted in order to determine the 
effect of lowered liquid viscosity on median drop size. This 
was done by spraying wax at temperatures up to 400 F into 
airstreams of the same temperature. The results are shown 
in Fig. 13 for comparison with earlier data at 300 F. There 
is a marked reduction in drop diameter as temperature 
increases. 

Points from the lines of Fig. 13 are crossplotted in Fig. 14 
to determine the viscosity exponent. An average for the 
three lines shown is 0.34. The same exponent was found for 
a pair of runs at 300 and 350 F using a smaller (3%5-in.) injector. 

On increasing temperature, properties other than liquid 
viscosity also change, and these changes can be expected 
to affect drop diameter. However, the effects are small 
compared to the overall effect observed. For example, as- 
suming the influences of surface tension and liquid density to 
be those reported by Nukiyama and Tanasawa (4), and as- 
suming the influence of air density to be the same as reported 
in the previous section, a temperature increase from 300 to 
400 F would cause an overall drop diameter change of less 
than 3 per cent if there were no effect of liquid or air viscosity. 
This compares to the 35 per cent changes observed. 


Correlation of Results 
The results described can be summarized briefly by the 
following proportion 
~ V —1-33y/0.08 70. 161) 70.34 (1 
PA 


where po is the density of air at 300 F and atmospheric pres- 
sure. Three other properties, not studied here, can be ex- 
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pected to affect drop size: Surface tension, air viscosity and 
liquid density. Assuming that drop diameter is proportional! 
to each property raised to some simple exponent, the ex- 
ponents can be determined uniquely by elementary dimen- 
sional analysis (since there are three exponents and threc 
dimensions). The result is 


X ~ pp 


Combining this proportion with the preceding one, the mass 
median drop Weber number (the ratio of air dynamic pressure 
forces to surface tension forces) can be separated to give 


A PAo 
(1 + 


08 70. pP 


XpaV? 19-34 
= K, ——— 
OL o 19-59 


In the experiments, pz was 0.828 (at 300 F) and pa, wa: 
0.000838. It is a convenient approximation to substitut: 
p,/10° for pa, If this is done, and if variables with lik« 
exponents are collected, one reasonable grouping of the 
variables is 


CL PL 


with all variables in consistent units. 

One final substitution can be made which eliminates D 
and v explicitly. Letting W be the mass rate of injection oi 
liquid (e.g., ]b/sec), then 


OL OL PL KL 


This equation, with K; = 0.61, was tested against all of the 
experimental data, including all air and liquid velocities for 
both costream and contrastream injection, all injectors, and 
all air densities and liquid viscosities. The only runs ex- 
cluded were those point samples designed to measure local 
spray variations. 

Fig. 15 shows the results for 101 runs about equally divided 
among the ,8;-, #5- and g,-in. injectors. For the two larger 
tubes, the standard deviation between observed and cor- 
related median drop diameters was 9.3 and 9.6 per cent. 
For the 35 points for the #;-in. injector, the standard devi- 
ation was 20 per cent. However, Fig. 15 shows clearly that 
there are large and consistent deviations for the nine 4; tube 
data points at Weber numbers lower than 25. Excluding 
these nine points, the standard deviation drops to 12 per cent. 

All nine of the excluded points resulted from 4y-in. tubes 
operating at relative velocities below 350 ft per sec. Earlier 
data showed repeatedly that the smallest injector did not 
fall into the general pattern at low relative velocities. There- 
fore, the failure of these points to fit the correlation is not 
surprising. 
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The drop diameters observed in this study were 35 to 100 
per cent larger than predicted by Nukiyama and Tanasawa 
(4) from their work in small pneumatic atomizers. How- 
ever, there was much better agreement (less than about 30 
per cent difference) with results cited by Marshall (2) for 
larger venturi atomizers. 


Other Spray Properties 


Because the entire drop size distribution was measured 
in each run, properties of the spray other than mass median 
diameter were also determined. 

The distribution of dispersions is shown in Fig. 16. Dis- 
persion is defined here as where and X45.9 
are the drop diameters at cumulative volume percentages 
of 84.1 and 15.9, respectively. No significant correlation 
could be found between dispersion and any of the variables or 
their simple interactions. The median dispersion was 
1.66, and the average was 1.71. About two-thirds of the 
data points lie within 10 per cent of 1.7, and 97 per cent lie 
within 25 per cent of 1.7. 

As would be expected, the ratio of Sauter mean diameter 
(the diameter of a drop with the same surface-to-volume 
ratio as the total spray) to mass median diameter decreased 
as dispersion increased. An empirical value for this ratio 
was (1.43 — 0.43 q) where q is the dispersion. The average 
Sauter mean diameter was about 70 per cent of the mass 
median, but the range extended from 36 to 87 per cent. 

Each drop size distribution was matched by a best fitting 
upper limit equation. The latter was proposed by Mugele 
and Evans (3) and is a modification of the log-probability 
equation that assumes a maximum drop size. By adjusting 
two constants of the upper limit equation (the third constant 
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was the mass median diameter, and the equation was con- 
strained to pass through the observed value), very good data 
matching was obtained. In only 2 per cent of the runs did 
the rms difference between predicted and observed volume 
percentage (at any given drop diameter) exceed 4 per cent 
absolute. Two typical matches between the upper limit 
equation and the experimental data are shown in Fig. 17. 
One match illustrates a very good fit with an rms error of 
0.4 per cent. The other illustrates a relatively poor fit with 
an rms error of 3.0 per cent. However, even the ‘‘poor’’ fit 
does a reasonably good job of describing the data. The 
experimental data were obtained as continuous traces on the 
Micromerograph; discrete data points from these traces 
are shown in Fig. 17 to distinguish more clearly between the 
data and the matching curves. 

A “maximum diameter” was one of the two constants 
selected to give the best fitting equation. In some cases, 
this maximum diameter was larger and in some cases smaller 
than the largest drops actually observed; however, the 
analytical data were quite erratic in measuring the last few 
per cent of the spray containing the largest drops. In gen- 
eral, the maximum diameter varied in the same way that the 
median diameter did, but the ratio of maximum to median 
showed a great deal of scatter. 

Since only three constants are required to define the upper 
limit equation, the maximum diameter is a unique function 
of the dispersion if the Sauter mean diameter is also specified 
as a function of dispersion. Assuming the Sauter mean to 
vary with dispersion as defined above, the ratio of maximum 
diameter to median diameter is approximately equal to 
(2.66q—0.56q?—1.10). This is an empirical expression useful 
for q up to about 2.5. 


Conclusions 


The atomization of liquids by large high velocity air- 
streams occurs by direct action of the airstream on the ex- 
posed liquid surface. Therefore, relative velocity between 
liquid- and airstreams is of primary importance. Physical 
properties of the fluids do affect spray fineness, but their net 
influence is less critical. The exact way in which the liquid 
is introduced to the air, i.e., the geometry and operation of 
the injector, is of least importance, particularly at very high 
air velocities and for customary ranges of the variables. 
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Some Thermal Stress Design Data 
for Rocket Grains 


M. L. WILLIAMS! 


California Institute of Technology 
Pasadena, Calif. 


Supplementing a previous paper (JET PROPULSION, June 1957) giving preliminary elastic stress 
concentration factors in internal burning star rocket propellant grains due to pressure differences 
across the web, similar design data for thermal stresses are presented for the particular case of a 
steady-state temperature difference across the grain, or for a specified heat flux. After determining 
the internal temperature using a voltage analogy, the Biot-Muskhelishvili analogy is employed to 
transform the thermal stress field into an equivalent temperature independent edge dislocation 
problem. The appropriate displacements are then applied to a series of photoelastic specimens of 
slotted grains, and the stress rise over that existing in a thick walled cylinder of the same web is 
evaluated experimentally. Following the previous procedure, the data are presented in design chart 


form. 


N A PREVIOUS paper (1),? some photoelastic data were 
presented for a series of internal burning star-shaped rocket 
grains subjected to pressure loadings. The design curves 
given were arranged in the form of stress concentration fac- 
tors, normalized upon the stress which would exist in a thick 
walled cylinder in the same web fraction, as a function of the 
geometric proportions and number of star points. Inasmuch 
as thermal loadings contribute another major class of design 
problems in rocket grains, some exploratory investigations 
were conducted upon geometries similar to the pressure loaded 
case, for a steady-state temperature difference across the web 
of the elastic specimen. Although the simulated loading con- 
dition is, of course, restricted, it is believed that the quantita- 
tive information obtained will provide a firmer basis for the 
necessary design extrapolations. 

In considering a typical model for elastic thermal stress 
analysis, it is found that the thick walled cylinder is the 
simplest configuration with which to work. With particular 
reference to rocket grains, Geckler’s work (2) supplies the 
basis for a fairly comprehensive study of this problem. On 
the other hand, and similar to the pressure loading case, when 
the wall thickness is nonuniform, as in the case of an internal 
burning star, analytical solutions are not readily obtained. 
Furthermore, numerical methods are not easily applied in 
areas of high stress concentration, and one therefore naturally 
turns to experimental methods. 

Biot (3) and Muskhelishvili (4) independently pointed out 
that certain classes of thermal stress problems, to which the 
usual grain configurations belong, are analogous to superim- 
posing appropriate displacements on the boundary which are 
uniquely determinable by the prescribed steady-state tem- 
perature distribution or heat flux into the grain. In 1939 
Weibel (5) exploited this analogy photoelastically and, after 
experimentally verifying the results by comparison with the 
analytical solution for the thick walled cylinder, evaluated the 
thermal stresses in a square thick walled tube. It required 
therefore only minor modifications to Weibel’s experimental 
setup to test a rocket grain configuration. Preliminary tests 
on several slotted grains have been completed wherein the 
maximum thermal stress at the star point was found and then 
divided, for normalization, by the thermal stress which would 
exist in a thick walled cylinder for the same temperature dif- 
ference. The resulting thermal stress concentration factors 
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can then be presented as functions of the design parameters in 
charts for easy use. 

More specifically, there are actually two concentrations 
present: (a) a geometric factor representing the difference in 
planar rigidity of the star from a thick walled cylinder of the 
same web fraction, and (b) a temperature factor due to the 
concentration of heat flux at the star points, which is of course 
not present in the cylinder. While the two effects can be 
combined, a comparison of the geometric factor with the 
(geometric) factor for pressure loading indicates they are 
quantitatively very similar. In order to explain this associa- 
tion, consider for example the hoop stress distribution through 
the web when the grain is loaded by internal pressure con- 
trasted with a temperature loading, say hot outside and cold 
inside. The first case tends to produce a more or less uniform 
tensile stress across the web, whereas the latter will lead 
toward internal tensile and external compressive stresses, 
approximately linearly distributed across the wall. 

If one next turns to Neuber’s curves for stress concentration 
factors (6) for a flat rectangular bar notched on one side (Fig. 
6)—with the notch associated with the star point at the 
minimum web section, it is found that when the bar is pulled 
in tension the factors are qualitatively similar to, but quanti- 
tatively lower than, those when the bar undergoes bending. 
As these loadings are characteristically those produced by 
pressure and temperature, respectively, in the web of the 
grain, one would expect the relative increase of stress at the 
star point to be in approximately the same relation. The ex- 
perimental results substantiate this conjecture (1)—higher 
values for the geometric thermal factor. 

It is perhaps worth emphasizing at this point that the tem- 
perature distribution is usually a requirement for the subse- 
quent determination of thermal stress* with one possible 


3 Generally speaking, the three-dimensional classical thermal 
stress boundary value problem requires a simultaneous solution 
of four coupled linear partial differential equations in the three 
displacements and the temperature. It is only a consequence of 
the fact that thermal inputs do not, in most engineering prob- 
lems, accelerate elements of the elastic body sufficiently for those 
displacements to affect the temperature distribution. Thus it is 
permissible to uncouple the four equations and solve the fourth 
(heat conduction) equation independently of the displacements, 
by hypothesis, and with the temperature field known, solve the 
remaining three equations for the displacements in three direc- 
tions corresponding to the temperature. See the full equations 
for example in the Parkus paper (7). The only area in rocket 
grains where this separation might have to be further explored is 
under deflagration-detonation conditions, where the chemical 
energy terms would also have to be included in the above set. 
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photothermoelastic exception (8). Even presuming that the 
temperature or heat flux boundary conditions are known at 
the surface of the solid, it is not always an easy thing to de- 
termine the internal temperatures. Aside from analytical 
solutions which exist in certain cases either by direct quadra- 
ture or conformal mapping techniques for steady state, one 
again thinks of experimental determination or analogs. 
Thermocouples have frequently been used, particularly where 
transients are required, but for steady-state and relatively 
simple boundary conditions a voltage analogy works well. 
This technique, using Teledeltos conducting paper‘ rather 
than a fluid medium, was actually the one used in preparing 
to find the thermal stress concentration factors. 

An alternate semi-analytical technique may also be used for 
determining approximate temperature distributions. Mini- 
:aum principles have been determined for heat conduction by 
Biot (9), Chambers (10) and Parkus (7), that is, a functional 
representation of the temperature distribution may be as- 
-umed, say a power series with unknown coefficients for which 
‘he values are determined by setting the first variation of the 
1eat conduction integral with respect to these unknowns to 
zero and solving the resulting algebraic equations. For those 
icquainted with structural analysis, the method is identical 
‘n principle with the Rayleigh-Ritz procedure. 


Analogy 


Consider the case of a hollow cylinder in plane strain, that is, 
the axial displacement is restrained to zero. Then, if steady- 
state temperature conditions exist and the temperature is 
symmetric and independent of the axial length, the problem 
can be formulated as one of plane stress. Specifically applied 
to rocket grains, one would have a thin slab of uniform thick- 
ness with outer and inner boundaries C, and C;, respectively 
(Fig. 1). The steady-state temperature distribution can then 
be calculated from a solution of V?7'(z, y) = 0 for the arbitrary 
boundary conditions on the temperature or heat flux upon C, 
and C;. Presuming a thermal field is the only loading, the 
stresses can then be calculated from the applicable equations 
of elasticity and the associated boundary conditions. 

The principle of using a photoelastic technique to find 
thermal stress was first employed by Weibel (5), based upon 
the independent theoretical work of Muskhelishvili (4) and 
Biot (3). The essence of the matter has also been concisely 
summarized by Timoshenko and Goodier (11), whose treat- 
ment wil] be outlined for ready reference. 

First, it can be shown that if the cross section of a cylinder is 
simply connected, i.e., without holes, then the steady-state 
thermal stresses in the plane are zero. Hence, it is concluded 
that in a cross section which is multiply connected, i.e., having 
one or more holes, thermal stresses present in this case can be 
completely relieved by making the cross section simply con- 
nected through the simple expedient of making cuts connect- 
ing the various inclusions. Reversing the argument then, the 
thermal stress distribution can be restored to the body by 
imposing the proper displacements and rotations in closing the 
various slits. 

For the thermal stress in the singly connected cross section 
to be zero, or = o, = Tz, = 0, the normal strain must be 


= = (1+ v)aAT =e [la] 
Ox 
(1+ v)aAT =e [1b] 
= 
and the shear strain 
ov ou 
zy 2 
(2) 


‘Such conducting paper is supplied by the General Electric 
Co., Schenectady, N. Y., under the trade name Teledeltos. The 
silver conducting paint was supplied by Micro-Circuits Co., New 
Buffalo, Mich. 
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Fig. 1 Geometry of the grain configuration 


The preliminary steps are completed upon recalling the defini- 
tion of rotation at a point 


1 f ov ou 
3 
i? 


From Equations [1 and 2] we have 


ou a ov ov ou 
or oy or oy 


which Cauchy-Riemann equations show that the u and v dis- 
placements in the x and y directions, respectively, are har- 
monic functions. Moreover the dilatation ¢ and rotation w 
are themselves harmonic, which is easy to verify, e.g. 


dy 2dy 2 or Lay oda 
Hence, one can form the complex function 


€(z, y) + y) = F(x + ty) = Fe) [4] 


From the property of line integrals along a path in the body, 
beginning at a point uo, %, one can show that 


z (de 


Now, whereas in a singly connected region these integrals 
must vanish when the contour is closed on itself, z —> 20, in 
order that the thermal stress vanish, in a multiply connected 
region this is not necessarily true. The values of the inte- 
grals—which depend only upon the given temperature dis- 
tribution—give the displacements Au and Av and the rotation 
Aw, which must be released at the cuts in order that the 
thermal stress vanish. Or, conversely, if these same dis- 
locations were imposed on a cross section (at room tem- 
perature), they would produce stresses identical to those of 
the temperature distribution with which they are associated. 

This means that if the original region is cut (say across AA, 
Fig. 1), then the thermal stresses may be simulated by im- 
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Fig. 2 Sectional thermal contour showing lines of constant 
voltage or temperature 


posing the resulting translations and rotation as boundary 
displacements at the cut. Using a deformeter similar to that 
of Weibel’s, a birefringent specimen in polarized light can be 
subjected to the dislocations, and resulting fringes interpreted 
in terms of principal stress by the Stress Optic law® 


2 
[7] 


where on the internal boundary, one of the principal stresses 
is zero. Thus the steady-state thermal stress variation along 
the boundary may be obtained relatively simply. 

The model stresses obtained in this plane stress experiment 
are now transferred to an equivalent plane strain case for a 
long restrained cylinder of the model material by increasing 
the stresses by the factor 1/(1 — v*), where v is Poisson’s 
ratio. 

Finally, the stresses in a prototype grain material are ob- 
tained by multiplying the plane strain model stresses by the 
ratio (E,/E,,) where E is Young’s modulus, and the subscripts 
p and m indicate the prototype and model, respectively. 
The axial stress when the ends are restrained becomes 


Czy = —E,a,T (2, y) + V,(Orp Fup) [8] 


Determination of the Temperature Distribution 


As most grain shapes have more or less irregular cross-sec- 
tional shapes, it is somewhat difficult to solve for the tem- 
perature distribution analytically. For the same reason, a 
numerical solution using relaxation methods is also somewhat 
cumbersome, because of the large number of irregular stars 


5The model constant 2f/t equaled 376 psi per fringe for the 
tests to be discussed presently: f is a property of the model ma- 
terial and ¢ is the thickness of the specimen. 
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needed in the network. The method used in this particular 
case was to take advantage of the similarity between the 
steady-state heat and potential equations. By cutting out a 
geometrically similar contour from a piece of conducting 
paper and imposing arbitrary boundary voltages, the lines of 
constant voltage (or temperature) were easily found by run- 
ning a probe, which was in series with a voltmeter, across the 
paper. 

In applying this technique to a geometry of interest in this 
investigation, the shape shown in Fig. 1 was chosen. First, 
the temperature distribution was measured experimentally 
using the Teledeltos conducting paper. A nominal 3 v dif- 
ferential was imposed by a standard dry cell after painting the 
boundaries C, and C; with conducting silver paint to insure 
the same level of potential along each boundary. A typical 
set of experimental temperatures or voltage contours is shown 
in Fig. 2. 

For grains which have a circumferentially periodic sym- 
metry, it turns out that along any radial cut, the radial dis- 
placement is identically zero. The displacement therefore is 
only tangential and proportional to the radius from the center 
of the grain, or alternately, a pure rotation of the slit. Now 
in Equation [5], F(z) = e+ iw = (1+ v)aAT+ iw. How- 
ever, the temperature must be an even function of 6 with re- 
spect to AA, whereas the rotation must be odd. Thus, upon 
writing the line element dz along a constant radial path ro as 
—ro sin@ dé + iro cosé dé the integral becomes in polar form 


Auy + = [a+v)a AT’, cosn6 +7 > wy sin nO} X 
(—ro sin 6 dO + iro cos 6 dO) 


which reduces by the orthogonality relations to 


Qn 
Au, + = iro [1+ v)a AT», cos 8 cos nO — 
wn sin sin nO 


0 + iro (constant) 


so that 
Au, = 0 [9] 


Ave = Bro = constant {10} 


Turning to Equation [6], one has 


ow ow 
Aw = f (2 dr + a0 


using the Cauchy-Riemann relation. The rotation then, 
B =} Aa, is 


B= (1+ v)ar J, 


and the physical significance of the rotation is immediately 
evident. It is the amount of heat flow across a constant 
radius boundary. Specifically, if the arbitrary radius 7 is 
chosen as the external surface of the grain, the analogy re- 
quires a rotational dislocation of the slit, proportional to the 
total heat flux into the grain. 

For cases where the flux is prescribed or known on the 
boundary through fluid mechanics considerations or other- 
wise, evaluation of Equation [11] is straightforward. On the 
other hand, if a constant temperature condition is prescribed 
the boundary flux will, of course, vary if the cross-section 
geometry is other than constant wall thickness. 


[11] 
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Fig. 3a Thermal concentration factor—four-slotted grain 
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Fig. 3c Thermal concentration factor—eight-slotted grain 


As a specific example useful in later normalization, the 
steady-state temperature distribution in a thick wall grain 
with web differential AT’, such that, T(a) = 0, T(b) = To, is 

AT _ tog (r/a) 
log (b/a) 


giving the heat flux as 


oAT To ! 


or log (b/a) 


so that the associated reference rotation for the cylinder be- 
comes, using Equation [11] 
2m(1 + vials 


For the grain of Fig. 1, Equation [11] was integrated nu- 
merically on the basis of the voltages measured and shown 
in Fig. 2 


2x 0(AT7'/7>) | 
nary 
(1 + vial d(r/re) dé 
o0(AV/Vo) 
| 


by the similarity of the temperature and voltage distributions. 
In the tests, the flux was measured across several different 
control boundaries rp and averaged. It was convenient to 
measure the voltage drop across a fixed Ary interval at V equal 
spaced intervals A@.6 The finite difference form of Equa- 
tion [14] therefore becomes 


B 
[14] 


roN 


Fae 


N 
B= (1 + > AV, {15] 


In such a fashion, the thermal concentration factor K; was 


® If one is interested only in total flux, as opposed to the com- 
plete determination of internal temperature distribution, it 
proves adequate and simpler to measure the integrated flow, be- 
cause the resistance of the conducting paper is constant over the 
interval Ar. 
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Fig. 3d Thermal concentration factor—ten-slotted grain 


obtained as the ratio of Equations [15] to [13] 


B ro Ab log (b/a) N 
K,= =“ av 16 
The test results for a series of four-, six-, eight- and ten-slotted 
grains are given in Fig. 3.7 


Alternate Methods 

Before leaving the subject of temperature distributions, it 
may be remarked that analytical expressions for temperature 
distributions may sometimes be obtained. One case, for 
example, is a wedge of included angle ¢ with, say, zero tem- 
perature on the radial edges and a temperature 7 along the 
unit radius. Using the even cosine expansion in @ about the 
bisector ot the wedge angle, one can easily find 


Indeed a closed form solution can be obtained if @ = 7/4 


(12) 
2r4 sin 40 


Such formulas may be useful in analyzing limit cases, although 
for the example chosen in Equation [17] the input heat flux is 
infinite because of the temperature discontinuity on the 
boundary. 

Sometimes boundaries of interest can be obtained by con- 
formal mapping techniques. In this case the mapping func- 
tion z = f(¢) = Cg (1 + még") deforms a circle into a wavy 
corrugated form with mean radius c, containing n waves or 
star points. The depth of the star point is mc measured from 
the mean radius (13). Thus this contour resembles an inter- 
nal burning surface and, if the web fraction is not too small, 
can be used to represent the grain cross section for tempera- 
ture or even for stress calculations. : 

As a final simple example, one can state a minimum princi- 
ple for heat conduction by analogy with the Rayleigh-Ritz 
(14) procedure used in determining approximate deflections of 
a membrane whose governing equation is also of the Laplace 


7 In obtaining the experimental data for this investigation, it 
is a pleasure to acknowledge the assistance of R. R. Parmerter. 
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Fig. 4a Isochromatic fringe pattern simulating thermal stress 
in a rocket grain 


Fig. 4b Isochromatic fringe pattern for a pressure loaded grain 


type. Theorems in a broader form have already been dis- 
cussed by various authors (7, 9, 10). 


The Dislocation Stresses 


Having determined the temperature, or more precisely the 
required rotation 8, the problem of imposing the rotation 
upon the plane stress photoelastic specimen remains. This, 
however, is a straightforward procedure and can be ac- 
complished without difficulty in the photoelastic laboratory. 
Using Weibel’s deformeter (5) as the loading mechanism, 
fringe patterns, and specifically the maximum fringes at the 
points of maximum stress, were recorded for several values of 
8. A typical fringe pattern is shown in Fig. 4 along with a 
similar pattern for pressure loading from (1). Inasmuch as 
the stress is directly proportional to deformation, this per- 
mitted a (linear) curve to be drawn through the test data in- 
stead of relying on measurements at only one loading. The 
test data then became a plot of maximum fringe order at the 
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point of the slot vs. 8B, from which the slope of the plotted 
curve could be readily obtained. 

Again, the plane stress slice with a constant wall thickness 
was used for normalization. It can be shown (15) that wedge 
opening will cause a tangential stress o¢ on the inside boundary 
of the model 


BE 2 log (b/a) 
Oom*(a) = E (a/b)? [18 


Recall, however (17), that when the stress in the cylinder ot 
the model material is desired, the plane strain solution is ob- 


tained from the plane stress solution Equation [18] by re- 
placing by — v?), so for the cylinder 


— v?) 1 — (a/b)? 


[18a 


Inasmuch as the radial stress is zero, i.e., ¢-(a) = 0, by the 
boundary condition, and as o, and gg are the principal 
stresses, the Stress Optic law (7) can be used to obtain 


BE [ _ 2 log 


and 
dn* E | 2 log oe 20) 
4x(2f/t) 1 — (a/b)? = 
From experimental data, however, we have a relation 
similar to Equation [19 or 18] for the actual contour; a 
geometric factor Ky, for identical 8, can now be defined 


Nmax dn /dB 


dn/dB 
= E 2 log 
4 w(2f/t) 1 — (a/b)? 


The results for a series of slotted grains are presented in 
Fig. 5. 


Thermal Stress—Constant Temperature 
Across Web 


Returning to the thermal stress problem in the grains, the 
final stress concentration factor will depend upon both K; and 
K,. The data can be arranged for convenient design use in 
the form shown in Fig. 7, which illustrates the case of a con- 
stant temperature difference across the wall. The actual 
maximum thermal stress at the star point in the cylinder 
under plane strain is, using Equation [18a] and recalling that 
the negative wedge angle must be used to obtain the desired 
doubly connected region 


E 2 log (b/a) 
= Ky E | B) (22) 


but from Equation [16], 8 = K,8*. Thereupon using the 
normalizing value of B* from Equation [13] 

EaTy | 
(1 — v) log (b/a) 


2 log (b/a) 
1 — (a/b)? 


oo. = KgK;y 3 [23] 
Upon setting Kz = K; = 1, Equation [23] is found to yield 
the proper stress for a hollow thick walled cylinder (16), 
which normalizing value a¢,* is shown in Fig. 8a as a function 
of b/a for convenient reference. 


Thermal Stress—Constant Heat Input 


Suppose now an alternate problem is considered, namely 
that various slotted configurations of the same outside diame- 
ter are compared under the condition that the total heat flux 
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1e into the grain is the same. Thus if 
or | 
K—| = [24] 
3] or 
then 
ld 
), fr? = - Qo 
ro Jo dé > 1h [25] 
so that Equation [11] gives the rotation to be applied as 
B = (1 + v)a2rbQo/K [26] 

ly But, as Qo is constant for all shapes, it will specifically be ~ 
e- the same for a thick walled cylinder, or 8* = 8, which fol- Fig 6. Stress concentration factors for notched bars, from 
1x lowing the previous notation implies K; = 1. In this fashion Neuber (6) 
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Fig. 7a Combined stress concentration factor—four-slotted grain 
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Fig. 7c Combined stress concentration factor—eight-slotted grain 


Equation [22] may be used as before, but now substituting 
Equation [26] there results for the maximum stress in the 
slot 


(27] 


= 


K — EaQs/k E 2 log(b/a) | 
21 ») 1 — (a/b)? 


where, as before, K, is taken from Fig. 5, and the normalized 
stress o¢,* = a¢,./K, is plotted in Fig. 8b. 

As a final remark, the heat flux is often estimated in terms 
of film coefficients, Q = h[T(b) —T,], where T,, is the 
(known) temperature at the medium surrounding the grain. 
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Fig. 7b Combined stress concentration factor—six-slotted grain 
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Fig. 7d Combined stress concentration factor—ten-slotted grain 


Conclusion 

The question always arises as to the validity of applying 
elastic stress concentration factors to viscoelastic propellant 
grain materials. Although it is not the intent of this paper 
to dwell upon this important fact, it may be pointed out that 
there are many practical cases in which a direct association for 
an elastic and viscoelastic state can be made (18); hence one 
has an important reason for devoting time to elastic 
considerations which at first sight appear inappropriate. 

Further development and engineering applications along 
these lines will be reported subsequently. 
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Fig. 8 Thermal stress in a thick walled cylinder, used as a nor- 
malizing quantity: 8a Constant temperature difference across 
web; 8b constant heat input to the web 
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Interplanetary Trajectories Under 


Low Thrust Radical Acceleration 


JACK COPELAND' 
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Interplanetary trajectories resulting from constant radial acceleration (toward or away from the 
sun) are obtained in terms of tabulated functions. For sustained acceleration greater than a critical 
value, the vehicle spirals to escape. For accelerations less than critical, the vehicle moves in pre- 
cessing orbits with perihelion or aphelion on the initial orbit using outward or inward thrust, re- 
spectively. Some examples of trajectories are given and their implications discussed. 


N THE many discussions of space travel that have ap- 
peared in recent literature, it is generally agreed that 
interplanetary corrective maneuvers will have to be made 
if the space vehicle is to intercept and orbit the destination 
planet. Whether the corrections are made continuously or 
by impulses, the guidance and control equipment must per- 
form their functions properly during months of unsupervised 
Received May 21, 1958. 
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operation. The question of system reliability is clearly of 
paramount importance, even if the guidance equipment is 
aided by human intelligence. 

One way to increase reliability is to make the guidance 
program so simple that a mistake is unlikely, or if a minor 
miscalculation is made, it will little affect the overall success 
of the mission. Since maneuvers must be made, the thrust 
must be changed at least in direction; however, an easily 
mechanized thrust program requiring a minimum of com- 
puter decisions for the major portion of the trip would allevi- 
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ate the problem. It should be emphasized that we are con- 
cerned here with the interplanetary phase of the guidance 
program; that is, escape from the Earth has been effected by 
some means and we are now faced with the problem of guid- 
ing the vehicle so as to intercept the destination planet. 

With simplicity a prime consideration, a constant radial 
acceleration, toward or away from the sun, looks quite 
attractive. From a mechanization standpoint, the thrust 
direction relative to the vehicle is easily determined. Mo- 
mentum conservation will reduce the vehicle velocity relative 
to the destination planet, thereby facilitating the corrective 
maneuvers necessary to effect final planetary capture of the 
vehicle. Although perhaps not optimum, when such factors 
as fuel economy or transit time are considered, such a thrust 
program because of its simplicity deserves more than cursory 
consideration. 

Fortunately, the equations of motion for this special case 
are not so intractable as to discourage further investigation. 
In the analysis that follows certain simplifying assumptions 
are made; consequently, the results are not exact, particularly 
in the region of the Earth and the destination planet. The 
general shape of the trajectories are given, however, along 
with the transit times, which are not prohibitive. 


Analysis 


A complete analysis of the problem would involve many 
variables among which are the effects of Earth and moon 
at departure, and of the destination planet at the end of the 
trip. For study purposes, however, it is sufficient to con- 
sider the problem reduced by the following assumptions: 

1 The perturbations of Earth, moon and destination 
planet are ignored; that is, the general n-body problem has 
been reduced to a simple two-body problem involving the sun 
and the vehicle. 

2 The vehicle is initially moving in what is essentially 
the Earth’s orbit with approximately the Earth’s orbital 
speed. 

3 The orbits of the Earth and destination planet are con- 
sidered circular and coplanar. 

Although these assumptions cannot be used for actual tra- 
jectory calculations, they do not invalidate the general con- 
clusions which result. 

Let us locate the vehicle in a heliocentric coordinate system 
by a radius vector R and a polar angle ¢. If the vehicle 
acceleration is constant and directed radially, the equations 
of motion are 


@R do \? Ro \? 
dt? -n(%) +0(*) 
d 
R? “e = constant [2] 


where 


Ry = one astronomical unit 
go = sun’s gravitational acceleration at Ro 
A the vehicle’s acceleration 


For convenience, we make the following changes of variable 


(1)? 


In other words, the unit of distance is the astronomical unit, 
the unit of time 1/2 m years, and the unit of accleration go. 
Velocity is measured in units of the Earth’s orbital velocity, 


* Numbers in parentheses indicate References at end of paper. 


268 


30 km/sec. With these changes, Equations [1 and 2] be- 
come 


1 
[4] 


rd = r(0)¢(0) = 1 [5] 


where the dots indicate differentiation with respect to r. 
The use of Equation [5] to eliminate ¢ from Equation [4] 
gives 

1 1 


=a [6 | 


Integration of Equation [6] assuming an initially circular 
orbit at r = 1, gives 
2 ] 
— —~—42a(r—-1)-1 [7] 
r? 


Except for one special case, the second integration to find 
r(r) or r(@) results in elliptic integrals; however, a great 
deal of information about the trajectory can be deduced from 
an examination of 7(r). 

It is apparent from Equation [7] that the radial velocity 
will be zero at those values of r that satisfy 


2 ar? — (2a + 1)r? + 27 —1 =0 [8] 
Equation [8] is satisfied by r = 1 and 


1+ V1-8a 
r= 9] 
4a 
It is clear that r is bounded for a < 4; however, if a > } no 
real roots exist, which means that escape from the solar system 
will be accomplished if the thrust can be maintained long 
enough. This behavior is well-known and has been men- 
tioned in the literature (1, 4). A dimensionless acceleration 
of 4 corresponds to an actual acceleration of about 0.8 X 
10-4g. Since the acceleration to be expected from the use of 
ion rockets may be of the order of 10~4 g (2), the trajectories 
for a < } are of more than academic interest. 
To obtain the trajectory ¢(r), Equation [7] must be inte- 
grated. Wehave 
— =r V(r 1) (2 ar? $1) [10] 
dp 


or) = [11] 
r = 
1 V(z —1)(2 ax? — 2 +1) 


The time tor may be found as follows 


(r) dr [12] 
1? 1 V(r — 1) (202? — x + 1) 


The radial speed is given by Equation [7] 


r 
and the @ component by Equation [5] 
1 
= — 


Therefore the velocity at r is given by 


r 


There will be four different types? of trajectories (4) corre- 


® Reference (4) examines the complete elliptic integrals for our 
case 2, and considers a< }. 
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“spirals” in toward the sun, reaching r = ? after about 
1.7 months. At this point 7 becomes positive and the vehicle 
returns to r = 1 after an additional 1.7 months at which time 
it leads the Earth by about 20 deg. If the vehicle were 


a <0; 0<a<}; a=}; a> 


om 


sponding to: 
Casel: a <0 


For a < 0,7 = 0 has three real roots, two positive and one 


negative. The negative root has no physical significance. It unperturbed by either Venus or the Earth, it would continue, 
may be noted that as a —> ©, the roots of * = 0 are at zero oscillating between r = 1 and r = 2 with a period of 3.4 
and unity; that is, if the sun were a point mass, collision with months. 
the sun by means of a constant finite radical acceleration 
would he impossible. Case 2: 0<a<} 

Here, the three roots of * = 0 are all real and positive. 

_ asin? 6+ b From (3) we find 
sin? 6+d 
2 6? a 

reduces Equations [11 and 12] to o= a II (« 5 r) + constant [19] 


2 
d= I1(6, a/b, k) + cF (6, »| + constant [l4 
yal b 


2 
T= [see k) + (11 — a3) E (6, k) + 


here 
a = —23k? 
b = 
— Ze 
— 


asin? 
esin? 6+ 1 


6 = (22 — 23)k? 


= (nm —2 — Xs) 


and a; > a» > 23 are roots of # = 0. The constants of integra- 
tion should be chosen so that @ and 7 vanish at r = 1; how- 
ever, this is not essential for plotting the trajectory. 
The functions F(0, k), E(9, k) and II(@, a/b, k) are elliptic 
integrals of the first, second and third kind, respectively (3). 
As an example let us take a = — 3. The three roots are 


= 3/4 = -3 [16] 


Fortuitously, a = — % results in perfect squares under the 
radicals and produces an 22 near Venus’ orbital distance. 
Equations [14 and 15] become (3) 


6 


tan 7/6 


l 5 
= F 6, + aretan 
( 4 4/1 — (1/16) sin?x/6 


where sin? 6 = 4(r — 3)/(r + 3). 
The trajectory $(r) is shown‘ in Fig. 1. 


| 


The vehicle 


4 The author is indebted to one of the reviewers for calling his 
attention to the existence of errors in evaluating the complete 
integrals of case 2. This prompted a complete recomputation 


of the trajectories for a = —% and a = }. Both trajectories 


differ significantly from those in the original paper. 
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(x — a2) sin 20 


— k? sin? 6 


4 + constant [15] 


k) — (21 — 25)B(8, k)] + constant [20] 


where 
a = — 
= 


Y= — 23) V/2 a(x; — 
asin? + b 


6 = 22-23 
= 


and 2; > x2 > x3 are roots of *(r) = 0. 
In this case, a = } is convenient for an iiiustrative example. 
Equations [19 and 20] become 


3 | 3 tan 6 


where r = 3 sin? 6+ 1. 


[17] 
tan 6 


aretan 


4 4/1 —(1/16)sin? 8 


sin 2 6 


18] 
8V1 


sin 7/3 


8 VW1—(1/16)sin? 


The trajectory ¢(r) is shown in Fig. 2. The vehicle spirals 
outward reaching r = 1.5 after 12 months. At this point, 
* becomes negative and the vehicle spirals inward, reaching 
r = 1 after an additional 12 months at which time it leads 
the Earth by about 108 deg. As can be seen from the 
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figure, the unperturbed trajectory is a quasi-elliptical orbit 
whose perihelion rotates about 108 deg per period. 


Case 3: a =} 


For this case the elliptic integrals degenerate into elemen- 
tary functions. We get 


r dz 
1 2(2—2)V2—-1 


The change of variable tan? £ = x — 1 permits easy integra- 
tion. The result is 


1 
af tan (= + :) | [24] 


[23] 


G(r) = 2 


Fig. 1 Trajectory for a = — 29, showing the relative positions of 
Earth and vehicle. The heavy dashed curve is the path the 
vehicle follows if unperturbed by either Venus or Earth 


Fig. 2. Trajectory for a = 14, showing the relative positions of 

Earth and vehicle. The vehicle reaches r = 1.5 after 1 year 

and, if unperturbed by Mars, continues on a path shown by the 
heavy dashed line 
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The time tor is 


x dx 


)=2 
1(2—2)Vz—-1 


[25) 


Here the change of variable sin? ¢ = x — 1 is convenient. 


We get 
= 4| 2 log, tan sin [26 | 


In this case 7 = 0 has three positive real roots, two of which 
are equal; that is, one root is unity, and r = 2 is a double root. 
Inspection of Equation [4] reveals that 7 (2) = 0; therefore, 
r approaches 2 asymptotically. The time to r = 1.5 is 
8.2 months. The trajectory $(r) is shown in Fig. 3. 


Case 4: a> } 


Here * = 0 has only one real root, unity, and a pair of 
imaginary roots. 
From (3) we get 


| 


/1 — k*sin? 0 — F(0, k) + E(6, k) 
sin 6 


1 
= are sin (k sin 8) — F(8, «| [28] 


Qa 
where 
9) 1 


— 
i 1 — cos 0 8a 


Fig. (4) shows $(r) for a = 3. 

The time to reach escape is perhaps of interest. The 
potential that must be overcome is 1/r._ When 1/r is equated 
to 3 (7 + r*d?), it is found that escape velocity will be 
reached when r = 1+ 4a. For a = 3, escape velocity is 
reached atr = 4. The time to escape is 16 months. 


t= 8.2 Mo. 
Ve 0.85 Vimars 


6.6 Mo, 
14 
/ t25.2Mo, 
r= 


Fig. 3 Trajectory for a = }¢ (although not apparent from the 
figure, the approach to r = 2 is asymptotic) 
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Fig. 4 Trajectory for a = %. After 16 months the vehicle 
reaches escape velocity atr = 4 


Discussion 


The question of what terminal corrective maneuvers must 
be made remains unanswered, since the planets were ignored 
to facilitate analysis; however, the destination planet may be 
replaced for discussion purposes without substantially alter- 
ing most of the vehicle trajectory. A qualitative discussion 
of the capture problem may then be made by assuming the 
interplanetary phase of the mission has been successful, so 
that near the nominal intercept point the vehicle is approach- 
ing the planet on what is approximately a collision course. 
Consider case 4, a = $, asa first example. Under the assump- 
tion of circular planetary orbits, the velocity of the vehicle 
relative to Mars in the region of r = 1.5 would be about 8 
km/sec. This is 3 km/sec in excess of escape velocity from 
the surface of Mars. Obviously, the vehicle in this case 
will either collide with Mars or orbit hyperbolically unless 
sufficient corrective impulse can be applied in the proper 
direction at the proper time. 

For the case of a = 3, the situation is not quite so critical. 
The vehicle velocity relative to Mars at r = 1.5 would be 
only slightly greater than escape (5.5 km/sec). Corrective 
maneuvers must still be made; however, less impulse is needed, 
and, since the vehicle approaches Mars more slowly, there is 
more time in which to make the necessary corrections. 
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For a = 3, the velocity of the vehicle at r = 1.5 is only 
4 km/sec relative to Mars at the same point. This is four- 
fifths of escape velocity from the surface of Mars. <A very 
attractive approach trajectory is one with aphelion somewhat 
beyond 1.5, say 1.55, and crossing r = 1.5 perhaps 2 months in 
advance of Mars. Mars, moving more rapidly than the vehicle, 
would then come between the vehicle and the sun as the vehicle 
approached aphelion. The vehicle would spend about 3 
months within some 5 million miles of Mars with a relative 
velocity of the order of 4 km/sec. Even at 5 million miles, 
however, the presence of Mars enters the problem only as 
a minor perturbation, since at that distance the gravitational 
acceleration of the sun exceeds that of Mars by four orders 
of magnitude. Whatever the approach trajectory, the 
terminal maneuver problem is a three-body one and can be 
discussed here only qualitatively. 

The situation for capture by Venus (a ~ %) is not so favor- 
able. The velocity of the vehicle relative to Venus is about 
12 km/sec compared to 10 km/see for escape from Venus, 
and only about 13 months are spent within 5 million miles 
ofr = 0.75. 

The bounded trajectories for a < } are interesting. In 
particular, a trajectory for a ~ } with its orbital period of 
2 years is practicable for a reconnaissance of Mars with sub- 
sequent recovery of the information from the vehicle upon its 
return near perihelion. Successful recovery would be a 
difficult problem but not manifestly more so than that, for 
example, of establishing the vehicle in a satellite orbit about 
Mars and transmitting the information there obtained to 
Earth, some 50 million miles distant. Reliable operation for 
a period of 2 years is indeed questionable; however, a mini- 
mum of guidance equipment would be required. Of course, 
the final feasibility of such a mission can only be determined 
after a more rigorous solution of the problem, including the 
effects of both Earth and Mars, has been made. Here, the 
assumption of coplanar orbits is not necessary or even desir- 
able; advantage should be taken of the 4 million mile sepa- 
ration that would result to reduce the perturbation of the 
vehicle by Mars. A similar mission for reconnaissance of 
Venus (case 1) is of equal interest. It has two attractive 
features: Duration of only about 33 months, and an un- 
obstructed communication range of less than some 75 million 
miles. 

The results given here for heliocentric motion may be easily 
extended to geocentric motion by a simple change of scale. 
For example, a vehicle leaving a circular orbit above the 
Earth under a constant } g radial acceleration would spiral 
outward approaching about 2 Earth radii asymptotically. 
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Reflection of Weak Shock Waves From 
Nozzles With No Flow and Critical Flow 


CHARLES E. FEILER! 


National Aeronautics and 
Space Administration 
Cleveland, Ohio 


The reflection of weak shock waves from nozzles having contraction ratios of 16, 4 and 2, and con- 
vergence half-angles of 60, 30 and 15 deg was studied in a shock tube. The experiments were per- 
formed with no initial flow through the nozzle and also with critical flow, and the results compared 
with theory. Contraction ratio, as predicted by theory, had much less effect on reflected shock 
strength with critical flow through the nozzles than with no initial flow. Convergence angle had 
little effect on reflected shock strength for either flow condition. 


IGH-frequency combustion instability poses a serious 
problem in the development of liquid propellant rocket 
engines. The mechanism that sustains these oscillations is 
obscure. For the case of longitudinal oscillations, the pres- 
sure waves reflect directly from the engine nozzle, and the 


have also been determined (6, 7). This paper describes the 
results of an experimental study of the reflection of weak 
shock waves from nozzles having contraction ratios of 2, 4 
and 16, and convergence half-angles of 15, 30 and 60 deg. 
Experiments were conducted for both no flow and critical 


flow initially through the nozzle, and the results compared 
with theory. The range of incident shock strengths (pres- 
sure ratio across the shock) varied from about 1.8 to 3.5. 


engine behaves acoustically as though it were a closed organ 
pipe (1, 2).2. In (3), it has been reported that the severity 
of combustion instability is related to the nozzle convergence 


angle. It was also postulated that the nozzle might have 
My an amplifying effect on the pressure wave associated with this Apparatus 
instability. 
3 The reflection of shock waves from the closed or open end The dimensions of the shock tube and locations of the pres- 


sure transducers, nozzle and injector are shown in Fig. 1. 
A pneumatically operated plunger was used to burst the 
diaphragm separating the high and low pressure sections. 
Driver pressure was measured with a Bourdon tube gage. 
The driver gas used was nitrogen; the gas in the low pressure 
or test section was air. Both gases were at room temper- 
ature. Diaphragm materials varied with the shock strength 
desired and were selected by trial. 


of a shock tube has been studied both experimentally and 
analytically (4, 5). Methods of analysis applying to nozzles 
for both no initial flow and critical flow through the nozzle 
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DRIVE STATION TEST 
SECTION SECTION 2 3 
DRIVE| | 
PRESSURE 7 22" 12" 12"— 12" 
CYLINDER” '\— DIAPHRAGM we PRESSURE 
PLUNGER REMOVABLE LUSH-MOUNTED TAP 
INJECTOR PRESSURE 
(SEE DETAIL TRANSDUCERS 
NOZZLE 
SHOCK-TUBE ASSEMBLY DETAIL 
BELOW) 
27] NOZZLE CONFIGURATION 
NOZZLE|CONTRAC-] DIAM,| a, [DIAM.| 2z, 
TION x, NI |DEG] y,IN IN. 
RATIO 
16 0.517] 30] 0.571.635 
2 4 1.033] 15] 1.132.210 
3 4 1.033] 30 | 1.13 |1.250 
4 4 1.033] 1.13] 810 
CROSS SECTION 5 2 1.462] 30 | 1.60] .913 


OF INJECTOR 
Fig. 1 Shock-tube assembly, injector and nozzle configuration 
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The configuration and dimensions of the nozzles are shown 
in Fig. 1. Contraction ratios (area shock tube/area nozzle 
throat) were 2, 4 and 16. 

A sketch of the injector used to provide airflow into the 
test section is shown in Fig. 1. For the air supply pressure 
used, the pressure ratio across the entry holes was subcritical, 
and thus the flow velocity into the shock tube was subsonic. 


Instrumentation 


High-frequency-response pressure transducers were mounted 
flush at four stations along the test section, as 
shown in Fig. 1. The outputs of these transducers were 
amplified and fed into a four-channel oscilloscope. Pressure 
at these four stations was recorded continuously on film. 
Accuracy of the transducers was within +1 per cent of full- 
scale deflection or about +2 lb per in.?. At the lower pres- 
sures, 25 to 30 lb per in.?, this error could be of the order 6 to 
8 per cent. Impact of the shock waves excited the natural 
frequency of the pick-up diaphragms. This effect, shown 
in the oscillogram of Fig. 2, increased the difficulty of meas- 
uring pressures accurately. 

Shock velocity was determined from the transit time of the 
shock between stations. Accuracy of velocity measurements 
was limited largely by time resolution on the film. Velocities 
appeared to be accurate to within about +3 per cent over the 
entire range (1500 to 1900 ft per sec). 


Procedure 


Data Recording 


The oscillogram of Fig. 2 is typical of the data obtained. 
Deflection of the traces occurs as the shock wave passes each 
station. Incident shock strength P; was obtained from the 
average value of pressure measured at all four stations. For 
no-flow tests, the initial pressure in the test section was 
atmospheric; for the flow tests, it was about 35 lb per in.? 
abs. A typical 2, é-plane diagram, as obtained from the 
oscillogram, is also shown in Fig. 2. The position or time of 
the interaction of the reflected shock with the contact surface 
varied with initial shock strength and usually occurred near 
station 3. As a result, reflected shock strength could be 
obtained from station 4 only. 


Calculation of Reflected Shock Strengths* 


The wave system set up by the shock waves is shown in 
Fig. 3. For the reflected wave (Fig. 3B) the following equa- 
tions may be obtained for a gas of heat capacity ratio equal 
to7/5 


5 1 
— x 
+1 


(1] 


U3 = M, 


Cy = +6) +6) 2] 
6P; +1 +1 

=A [3] 

+ Un? = 5Ca? + Us? [4] 


Equations [1 and 2] relate properties behind the reflected 
shock to the pressure ratios across the incident and reflected 
shock waves. Equations [3 and 4] relate flow properties 
(assumed to be steady-state) behind the reflected shock to 
those at the nozzle throat, that is, across the converging sec- 
tion of the nozzle. For no initial flow through the nozzle, 
M, = 0, calculation showed that for the weakest incident 


*A list of definitions of symbols appears at end of paper. 
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Fig. 2 Typical oscillogram and x, (/-plane diagram. Nozzle 1: 
P, = 3.19; Pr = 250% Uy = 0.93 
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(A) WAVE SYSTEM BEFORE REFLECTION 
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(B) WAVE SYSTEM AFTER REFLECTION 
Fig. 3 Wave system in shock tube 


shock strength (P; = 2.0) sonic flow occurs at the nozzle 
throat after shock reflection for contraction ratios of about 
1.4 or larger. Thus for the range of contraction ratios 
and incident shock strengths studied, the nozzle was choked 
after reflection of the incident shock. The boundary condi- 
tion was then Us; = Cg. For this boundary condition there 
results a reflected shock wave; however, for weaker incident 
shock strength or smaller contraction ratios the flow through 
the nozzle would be subsonic, and the reflected wave would 
be an expansion wave due to the different boundary condition. 

For critical flow through the nozzle initially, the boundary 
condition was also Us, = Cm. In this case the value of M, 
is a function of the contraction ratio and can be calculated 
from stream-tube area relations. M;, was also equal to M;, 
where M; is the Mach number of the flow behind the reflected 
shock. 
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Fig. 4 Variation of reflected shock strength with incident shock 
strength for various nozzle configurations. No initial flow 


The system of equations was solved by iteration. Values 
of P, calculated for the two cases are shown in Table 1. 


Experimental Results 


Reflected shock strength is shown as a function of incident 
shock strength in Figs. 4 and 5 for the initial conditions of no 
flow and critical flow, respectively. No data were obtained 
for nozzle 5 with flow, since its critical flow rate consider- 
ably exceeded the capacity of the available system. Con- 
traction ratio had a much greater effect on reflected shock 
strength with no initial flow than with critical flow initially, 
as predicted by the theoretical calculations. Convergence 
angle gave no observable effect for either flow condition. 

The scatter in the data is of the same order as the 
experimental error previously discussed. Experimental data 
tended, in general, to fall below the corresponding theoretical 
curves. A recent study of the interaction of a reflected 
shock with the boundary layer created by the incident shock 
shows that for a closed-end tube the reflected shock is appre- 
ciably attenuated (8). This effect may partly account for 
the observed deviation. 

The data of Fig. 4 are crossplotted in Fig. 6 to show the 
variation of reflected shock strength with contraction ratio. 
The curves terminate on the abscissa at the contraction ratio 
that just gives sonic flow for that particular incident shock 
strength. 
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Fig. 5 Variation of reflected shock strength with incident shock 
strength for various nozzle configurations. Initial flow critical 


274 


Table 1 Calculated values of reflected shock strength 


No flow, Sonic flow, 
A,/A, M; P, M, = 0, Pr M, = M;, Pp 
2 0.305 1.5 are 1.40 
2.0 1.25 1.75 
2.5 1.49 2.04 
3.0 1.70 2.29 
+ 0. 1466 1.5 1.20 1.44 
2.0 1.53 1.81 
2.5 1.83 2.14 
3.0 2.10 2.43 
16 0.0363 1.5 1.39 1.46 
2.0 1.78 1.86 
2.5 2.13 2.21 
3.0 2.43 2.52 


Conclusions 


The results show, as predicted by theory, that contraction 
ratio has much less effect on reflected shock strength when 
there is critical flow initially through the nozzle than when 
there is no initial flow. Also, convergence angle had no 
observable effect for either flow condition. The insensitivity 
of reflected shock strength to contraction ratio with critical 
flow also supports the observation that a rocket engine behaves 
acoustically as a closed organ pipe during longitudinal insta- 
bility. 
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Nomenclature 
A, = area of shock tube 
A, = area of nozzle throat 
c = speed of sound 
M = Mach number 
P, = pressure ratio across incident shock 
Pp = pressure ratio across reflected shock 
u = particle velocity or contact surface velocity 
Ci; = c;/c; (e.g., Ca = ¢5/e1), dimensionless quantity 
M; = u;,/c;, dimensionless quantity 
U;; = u;/c;, dimensionless quantity 
INCIDENT 
26- _ SHOCK 
STRENGTH, 
——— THEORY 
2.2 EXPERIMENTAL 
REFLECTED 
SHOCK 
STRENGTH, 
Pr 


4 .4 8 1.0 
RECIPROCAL OF CONTRACTION RATIO 


Fig. 6 Variation of reflected shock strength with contraction 
ratio. No initial flow 
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Subscripts 

i,j = regions of constant flow properties (Fig. 3) 
1 = low pressure or test section initially 

2 = behind incident shock 

3 = behind rarefaction wave 

4 = high pressure or driver section initially 

5 = behind reflected shock 

6 = nozzle throat 
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Turbulent Flame Studies in a 


Two-Dimensional Open Burner 


J. H. GROVER,’ E. N. FALES® 
and A. C. SCURLOCK* 


Atlantic Research Corp. 
Alexandria, Va. 


A two-dimensional open burner has been developed for studying turbulent flames. Stoichiomet- 


ric natural gas-air flames with inlet flow velocities of 15, 20, 30 and 40 ft per sec were investigated. 
Densitometer traverses of time exposures were used to locate the position of the mean flame front 
and its root-mean-square displacement. Flash tube exposures of the flame with one-micron mag- 
nesium oxide particles in the flow were taken to determine the flow lines into and through the flame 
front. The large reduction of brightness that occurred when the flow passed through the flame 
front, due to the reduction of particle concentration, shows clearly the outline of the instantane- 
ously wrinkled flame front and the islands of unburned gas often present in the burned gas. From 
the intercepts of the flow lines and the mean flame front, the local turbulent flame velocity along 
the flame front was calculated. The turbulent flame velocity increased rapidly with distance above 
the burner rim. The flame velocity increase can be entirely accounted for by flame wrinkling. 
Both this result and the appearance of the flash tube photographs, showing the instantaneous flame 
front, support the wrinkled flame concept of turbulent flame propagation. 


pipe turbulence. No account was taken by Karlovitz, how- 


N A REVIEW of the status of experimental turbulent 
flame studies prepared by two of the present authors and 
published in 1954 (12)> attention was directed to the inade- 
quacy of the then existing experimental results for turbulent 
flames propagating in homogeneous mixtures (as distin- 
guished from turbulent diffusion flames). It was pointed out 
that most of the earlier studies (1,3,13,16) with such flames 
were concerned with the overall effects of turbulence on the 
entire flame brush and had not attempted to follow the local 
development of the flame at various distances from the stabi- 
lizer. One notable exception was the work of Karlovitz and 
his associates (8) who measured turbulent flame velocity vs. 
distance from the burner rim for bunsen flames disturbed by 
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ever, of the flow divergence in the unburned gases; and it 
now appears that this divergence is of considerable impor- 
tance. Since the present authors in 1953 (11) outlined their 
concept of the role of eddy diffusion, and thus of flame ex- 
posure time or distance downstream from the stabilizer, on 
the development of the turbulent flame, it has been recog- 
nized that a better understanding of turbulent flames could 
be achieved from analysis of measurements of local turbulent 
flame velocity vs. distance from the stabilizer under various 
conditions of turbulence, mixture composition and _ flow. 
The recent work of several investigators has been directed 
toward this objective (2,10). 

The purpose of the present investigation of turbulent flames 
propagating in homogeneous mixtures has been to determine 
the location of both the time-average and instantaneous 
flame fronts, to obtain measurements of the flow in the vicin- 
ity of these flame fronts, and from these data to determine 
accurately the local turbulent flame velocity point-by-point 
along the mean flame front. This has been accomplished 
through the use of a “two-dimensional” open burner. 
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Fig. 1 Schematic of two-dimensional open burner 
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Fig. 2 Photographs of flame and flow in two-dimensional open 

burner (stoichiometric natural gas-air mixture; approach veloc- 

ity, 20 ft per sec; four-mesh screen 2 in. upstream of burner 
rim) 
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Description of Apparatus 

The 1 X 3-in. rectangular burner that was employed in 
these studies is shown schematically in Fig. 1. The 3-in. 
sides of the burner are open. Transparent confining walls 
extend along the 1-in. sides (the front and back of the burner). 
These serve to constrain the flow, causing it to be essentially 
two-dimensional] and making it possible to view and photo- 
graph the flame. 

Air is supplied from an oil-free compressor at a pressure of 
approximately 60 psi. The compressed air is dried by passage 
through a separator for removal of entrained moisture and 
through a silica gel bed for absorption of water vapor. The 
air is then metered by rotameter and passes through a contro! 
valve which regulates the flow rate and reduces the pressure. 
The air then enters an inspirator, as the primary stream, where 
the air pumps in low pressure natural gas (about 88 per cent 
methane), and the two streams are thoroughly mixed. 

The natural gas, obtained from the Alexandria city main, 
is also metered by rotameter. 

The air-fuel mixture proceeds from the inspirator into the 
‘alming section. The gases then pass through a nozzle with 
a contraction ratio of 25:1 where the turbulence intensity is 
reduced to less than 1 per cent. The 1 X 3-in. exit of this 
nozzle serves as the entrance to the burner. 

Turbulence in the gases entering the burner was generated 
by placing a screen in the nozzle throat. Screens employed 
in this study were of stainless steel with a ratio of mesh width 
to wire diameter of 5:1. 

As is shown in Fig. 1, removable spacers extend downstream 
for a short distance along the 3-in. sides of the burner, the 
purpose of the spacers being to make possible a variation of 
the distance between the flame and the turbulence generating 
screen. Immediately downstream of the spacers, propane- 
air pilots extend along the 3-in. sides of the burner. 

To prevent the flame from stabilizing on the glass walls, it 
was necessary to maintain them at a relatively low temper- 
ature. This was accomplished by employing double plates 
with circulation of cooling water between the inner Vycor 
glass and the outer window glass. 

For visualization of the flow into and through the flame, 
1-u particles of magnesium oxide were injected into the com- 
bustible mixture just prior to its entry into the calming sec- 
tion. These particles, which upon entering the burner are 
evenly distributed throughout the flow, are illuminated by a 
short-duration flash and photographed. The flash duration 
is regulated so that the individual particles appear in the 
photographs as short white streaks, which indicate the local 
instantaneous direction of flow. 

Gilbert, Davis and Altman (6) concluded that 1-y particles 
would correctly trace the streamline direction through a 
propane-air flame. It is thus believed safe to assume that the 
magnesium oxide particles follow the flow in these flames with 
at least the precision to which the direction of the particles 
are measured. 

To obtain a photograph of these very small particles with- 
out obscuration by radiation from the flame, the particles 
were illuminated by the light of a cloud-chamber flash tube. 
The flash tube is powered by a bank of eight 25-uf capacitors. 
These capacitors are charged to a potential of 2000 v, and their 
discharge through the flash tube is synchronized with the 
camera shutter. The flash duration is about 0.8 millisec. 

The light flux from the flash tube is intensified by focusing 
it at the center of the burner with a pair of cylindrical lenses. 
The lenses are focused so that a flat beam of light about 2-in. 
thick slices through the center of the 3-in. wide flame. 

The camera was constructed around an f8, 36-in.-focal- 
length, aerial-survey camera lens. The lens was equipped 
with a high speed between-the-lens shutter with a speed of 
about 7 millisee which was synchronized with the flash tube 
and used in obtaining the particle tracking exposures. In 
addition, a pneumatic shutter was installed back of the lens 
for obtaining time exposures. 
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Fig. 3 Brightness distribution for the densitometer traverse 

2.42 in. above the nozzle (stoichiometric natural gas-air flame 

with approach velocity of 20 ft per sec and four-mesh screen 
2 in. upstream of burner rim) 


Experimental Results 


All of the present work has been conducted with stoichio- 
metric natural gas-air flames. The principal results were ob- 
tained with flow velocity at the entrance to the two-dimen- 
sional burner maintained at 20 ft per sec and with a four-mesh 
screen situated in the nozzle throat 2 in. upstream of the 
burner rim. Flame photographs and some limited results 
were also obtained at other approach flow velocities and with 
an eight-mesh screen in the nozzle throat. 

Shown in Fig. 2 are examples of the two types of photc- 
graphs employed for analysis in this study. The first is a 
time exposure of the turbulent flame with the flame itself as 
the sole source of light. The diffuse nature of the flame 
bush and its widening with distance from the burner rim 
are clearly observed. The second is a flash tube exposure 
with magnesium oxide particles in the flow. Duration of the 
flash was chosen so that each particle would appear on the 
photograph as a short dash indicating the local instantaneous 
direction of flow. The instantaneous boundary between 
burned and unburned gases, assumed to be the flame front, 
is clearly shown. This results from the approximately eight- 
fold (7.62) increase in specific volume of the gases associated 
with combustion and the attendant eightfold decrease in 
volumetric concentration of magnesium oxide particles upon 
passage through the flame front. 

The flame front is observed to have a wrinkled appearance. 
At its top, a small island of unburned mixture has been 
pinched off by the flame propagation and is now surrounded 
by burned gases. The presence of a continuous sharply de- 
fined flame front in these and similar photographs (see also 
Fig. 4) lends support to the wrinkled flame concept of turbu- 
lent flame propagation as distinguished from the picture ad- 
vanced by some (14) of a diffuse reaction zone. Since these 
photographs are essentially two-dimensional, depending for 


¢ Here the speed of the camera shutter is sufficient to prevent 
radiation from the flame itself from exposing the film and mask- 
ing the light reflected from the particles. 
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Fig. 4 Typical flash photographs of magnesium oxide particles 
in a stoichiometric natural gas-air flame with approach velocity 
20 ft per sec (four-mesh screen 2 in. upstream of burner rim) 


exposure only on the reflected light from small particles mov- 
ing with the flow in a thin flame slice equidistant from the 
burner walls, they are less ambiguous than photographs ob- 
tained by shadow or Schlieren techniques which are dependent 
on refraction of light rays during their passage through the 
entire flame zone. 


Analysis of Time Exposure for Position of Mean Flame 
Front and Its rms Displacement 


From the negative of the time-expesure photograph—the 
first photograph of Fig. 2—both the mean position of the 
flame front and the rms displacement of the flame front from 
its mean position have been determined from densitometer 
traverses across the flame. The mean flame front was taken 
as the locus of points of maximum density. 

In order to determine the rms displacement of the flame 
front, it was necessary to relate the film density to the rela- 
tive luminosity or brightness of the flame. Film density was 
found to be proportional to 0.64 times the logarithm of the 
luminosity over most of the range of density. Using this rela- 
tionship and the density traverses, the distribution of flame 
luminosity was obtained. In Fig. 3, this distribution is 
shown for one of the density traverses. The logarithm of the 
relative brightness is plotted vs. the square of the distance 
from the mean flame front. Except for the dim regions of 
the flame where the background and reflected light are not 
negligible compared with the light of the flame, the data fall 
on a straight line. Thus the distribution of brightness is 
essentially Gaussian, and the rms displacement occurs at the 
point where the luminosity is 61 per cent of the maximum.? 


7 Since the distribution is Gaussian 


—a\? 
b => 
oV 2a 
where b is relative brightness, o is a constant equal to the rms 


displacement; x is distance, and a is distance at the mean; and 
=e-'; = 0.6065. 
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Fig.5 Mean flow lines, mean flame front and rms displacement 

of the mean flame front (stoichiometric natural gas-air flame 

with approach velocity of 20 ft per sec and four-mesh screen 
2 in. upstream of burner rim) 


Since the density is proportional to 0.64 times the logarithm 
of the luminosity, the rms displacement occurs at the point 
where the film density is 0.139 below its maximum value;* and 
the rms displacement can be obtained directly from the 
densitometer traverses. 


Analysis of Flash Tube Photographs for the Mean Flow 


From a set of 20 flash tube exposures taken at different 
times under conditions corresponding to those for the second 
photograph of Fig. 2, the mean flow profile in the neighbor- 
hood of the flame has been obtained. A set of eight of these 
flash tube photographs used in the analysis are shown in Fig. 
4. The photographs were projected one-by-one on a large 
grid, and the angle made by the particle tracks with the 
vertical was measured at 165 points in the flow. The read- 
ings for the right half of the flame were averaged with those 
for the left half, so that the average angle for each point was 
the average of 40 readings. From these data, curves of 
average angle vs. lateral distance from the center of the flame 
were drawn for each 0.5-in. increment of distance downstream 
from the burner inlet. From these results, mean flow lines 
have been constructed. 

These mean flow lines, with the mean flame front and its 
rms displacement superimposed, are shown in Fig. 5. The 
rms displacement is observed to increase rapidly with dis- 


8 Da — Da = y log (ba/ba) = 0.64 log 0.6065 = —0.139, where 
De and Da are the film densities at the rms displacement and 
the mean, respectively, and y is the measured proportionality 
factor between the film density and the logarithm of the relative 
brightness. 
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tance from the burner rim. The flow divergence in the 
gases approaching the mean flame front is found to be 
substantial. 


Determination of Local Turbulent Flame Velocity 


From the mean flow lines and position of the mean flame 
front, and assuming conservation of mass within any one mean 
stream tube, the local turbulent flame velocity S7 has been 
determined as a function of distance measured along the 
flame front from the burner rim.’ The results are shown 
in Fig. 6. Sr is found to increase rapidly from 1.5 ft per sec 
at 1.5 in. above the burner rim to about 7 ft per sec at ap- 
proximately 4 in. above the burner rim at the flame tip. 
These values may be compared with a laminar flame velocity 
Sx for the mixture of 1.2 ft per sec (9). 

The time of exposure of a mean flame element or its age ¢ has 
also been estimated from flow velocities obtained from the 
known divergence of the central stream tube and the assump- 
tion that, in the approach flow, the component of flow ve- 
locity parallel to the burner axis at any given distance 
downstream of the chamber entrance remains uniform. 
Flame element age is obtained by graphic integration of the 


relation 
where 
h = the distance downstream of the chamber inlet 


the local velocity in the central stream tube 


S 
II 


A plot of V vs. h is shown in Fig. 7, and the Fig. 6 results re- 
plotted as S7/S_, vs. t are shown in Fig. 8. From Fig. 7 the 
flow divergence in the gases above the burner inlet is found 
to result in a reduction in flow velocity parallel to the mean 
flame front from approximately 20 ft per sec at the burner 
entrance to about 8 ft per sec at the flame tip. 


® For any one stream tube, the assumption of conservation of 
mass yields puAsU = puAsS7, where px is the density of the un- 
burned gases; A, is the cross-sectional area of a stream tube at 
the chamber entrance (assumed to be uniform); A; is area of the 
mean flame front through which the stream tube passes, and 
Sr is the local turbulent flame velocity (actually the local turbu- 
lent flame velocity averaged over the area A,). 


Laminar flame velocity 


TURBULENT FLAME VELOCITY — ft/sec 


2 3 4 
LENGTH ALONG FLAME FRONT—inch 


Fig. 6 Turbulent flame velocity vs. flame length (stoichiometric 
natural gas-air flame with approach velocity of 20 ft per sec and 
four-mesh screen 2 in. upstream of burner rim) 
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Fig. 7 Velocity along central stream tube (stoichiometric 
natural gas-air flame with inlet flow velocity of 20 ft per sec 
and four-mesh screen 2 in. upstream of burner rim) 


Position of Mean Flame Front and Its rms Displacement 
From Flash Tube Photographs 


As a check on the results obtained from the densitometer 
traverses of the time exposure photograph, the position of the 
mean flame front and its rms displacement were determined 
independently by analysis of a set of flash tube photographs. 
The outline of the instantaneous wrinkled flame front, as 
denoted by the rapid change in particle concentration, was 
traced for 18 flash exposures of this flame. At short in- 
tervals downstream from the inlet, measurements were ob- 
tained from each tracing of distance of the right and left 
flame fronts from the chamber center line. When flame 
wrinkling resulted in the presence of more than one flame 
front at any chosen height, a measurement for each individual 
flame front was taken. Flame symmetry was again assumed, 
and measurements from right and left flame fronts .were 
grouped together for analysis. 

The position of the mean flame front and its rms displace- 
ment determined for each interval from these measurements 
are plotted as points in Fig. 9. The results of the densi- 
tometer method are shown as the solid lines for comparison. 
Agreement is excellent on the position of the mean flame front 
and reasonably good on rms displacement. Since the densi- 
tometer method is dependent for flame front location on light 
from the luminous flame zone, whereas the flash tube photo- 
graph method is dependent for flame front location on the 
rapid change in concentration of suspended particles, the 
excellent agreement of the two methods on location of the 
mean flame front gives strong support to the assumption that 
the steep gradient in suspended particle concentration does, in 
fact, correspond to the luminous flame front. Since the flash 
tube photographs are of a thin slice of flame in the center of 
the burner, whereas the time exposure photographs are of the 
entire flame including that adjacent to the front and rear 
walls, this excellent agreement is also interpreted as evidence 
that, under the conditions studied, the boundary effects of 
the confining walls are relatively small. 
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Fig. 8 S7/S, vs. time (stoichiometric natural gas-air flame 
with approach velocity of 20 ft per sec and four-mesh screen 
2 in. upstream of burner rim) 
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Fig. 9 Mean flame front and rms displacement as obtained by 

two methods (stoichiometric natural gas-air flame with approach 

velocity of 20 ft per sec and four-mesh screen 2 in. upstream of 
burner rim) 
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Fig. 10 Flash photographs of magnesium oxide particles in two- 

dimensional open burner flames at various flow velocities (stoich- 

iometric natural gas-air mixture; screen 2 in. upstream of burner 
rim) 
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Instantaneous Flame Surface From Flash Tube 
Photographs 


By further analysis of the flash tube photographs, an addi- 
tional test of the wrinkled flame concept has been conducted. 
At flow velocities of both 20 and 40 ft per sec, and employing 
the same four-mesh turbulence generating screen as before, 
the mean instantaneous flame area has been approximated 
from flame perimeter measurements obtained from flash tube 
photographs taken at these two flow velocities. From thes: 
areas and the known mass flow at the inlet, laminar flame 
velocities have been calculated for comparison with th« 
known value of 1.2 ft per sec. In order to convert the mean 
instantaneous flame perimeter, obtained by averaging thi 
individual perimeter measurements from a set of flash photo- 
graphs, into a mean instantaneous flame area, an assumption 
must be made in regard to flame wrinkling in the direction 
normal to the flame slice depicted in the photographs. It 
would seem reasonable, in the absence of other information, 
to assume that flame wrinkling in this direction is, on th« 
average, the same as that in the plane of the photographs. 
This was the assumption chosen. Thus, to determine the 
mean instantaneous surface area per unit of flame depth, the 
mean instantaneous flame perimeter is multiplied by the ratio 
of mean instantaneous flame perimeter to perimeter of the 
mean flame front. The results of the analyses are given in 
Table 1. At an inlet velocity of 20 ft per sec, a flame ve- 
locity of 1.35 ft per sec based on the mean instantaneous 
flame area is obtained. At an inlet velocity of 40 ft per sec, 
the flame velocity similarly determined is 1.30 ft per sec. 
These compare well with the known laminar flame velocity 
of 1.2 ft per see and indicate that, at least under the condi- 
tions of these experiments, the turbulent flame is nothing 
more than a randomly fluctuating wrinkled laminar flame. 


Table 1 Flame velocity from mean instantaneous flame | 


area 
Inlet velocity, ft/sec 20 40 | 
no. flash tube photographs analyzed 18 19 
mean instantaneous flame perimeter, in. 10.8 19.13 
standard deviation, in. Le 1.3 
perimeter of mean flame front, in. 19 HS 


mean instantaneous 
flame perimeter 
perimeter of mean 
flame front 

mean instantaneous flame area per unit 
of flame depth, in.?/in. 14.8 30.8 

flame velocity based on mean instantane- 
ous flame area, ft/sec 

laminar flame velocity (9) Sz, ft/sec 


wrinkling factor = 


Flame Photographs Over a Range of Flow Velocities 


Shown in Fig. 10 are flash tube photographs taken at inlet 
velocities of 15, 20 and 30 ft per sec. In the top set of three 
photographs the four-mesh screen was employed for turbu- 
lence generation. An eight-mesh screen was used in the lower 
photographs. Flame height and turbulent flame velocity 
are found to increase with flow velocity. For a given flow 
velocity, flame height is always greater with the eight-mesh 
screen, and turbulent flame velocity is always found to be 
greater with the four-mesh screen. Qualitatively, these 
results are as would be expected and are consistent with those 
obtained previously by Wohl et al. (16). 


Comparison With Theory 


Scurlock and Grover (11, 12), in the development of their 
theory, consider an initially flat flame of infinite extent 
propagating into a combustible mixture agitated by isotropic 
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Table 2 Calculated values of J, and C, 


Distance above burner rim, in. 1.2 if 
measured Sr, ft/see 1.4 TR 
measured (Y?2)'/2, in. 0.05 0.083 
Dryden ls, in. 0.0486 0.0493 
calculated /, assuming C; = 4, in. 0.165 0.250 
ealculated C; using Dryden 0.348 0.155 


2.0 2.5 3.0 3.5 
1.8 2.2 2.9 42 
0.123 0.165 0.211 0 275 
0.0511 0.0526 0.0541 0.0556 
0.220 0.221 0.193 0. 164 
0.216 0.228 0.314 0.459 


turbulence. Only turbulence with a scale very large relative 
to the thickness of the laminar flame front is considered, the 
turbulence assumed only to wrinkle and extend the flame 
surface, the laminar flame velocity remaining constant and 
unaffected either by any increase in transport properties or 
by flame curvature. It is proposed that as time passes the 
initially flat flame becomes more and more wrinkled by eddy 
diffusion of flame elements away from the mean position of the 
flame, and the flame area at first progressively increases. 
Assuming that the depth of the wrinkles is proportional to 
the rms displacement of a flame element from its mean posi- 
tion (Y?)'/? and that the average base width of these wrinkles 
is proportional to the Eulerian scale of turbulence (2, the 
following relation for the ratio of wrinkled to initially flat 
flame area A7/A,, was obtained from simple geometric con- 


siderations 
A T Sr y2 
St 


where C, is a combined proportionality constant somewhat 
arbitrarily assigned a value of 4, based on the assumption that 
mean wrinkle depth equals 2(¥2)'/?, and that the average 
base width of the turbulent flame wrinkles is about four 
times the scale of turbulence. Although neither of these as- 
sumptions can be considered quantitatively exact, geometric 
considerations and previous experimental results indicate 
that they should be fair approximations. 


Indication That Scale of Flame Turbulence 
Is Substantially Larger Than Anticipated 


Since experimental values of S; and (¥?)'” are available 
from the present investigation, Equation [1] may be em- 
ployed either: 

1 With the constant C; assigned the value of 4, after Scur- 
lock and Grover, to calculate the scale of turbulence (, at 
various points along the flame front for comparison with that 
predicted by the relation of Dryden et al. (5) for turbulence 
behind screens in cold fiow. 

2 With the values of /2 predicted by the Dryden relation 

to calculate values of C}. 
The results of such calculations are shown in Table 2. For 
C, = 4, the calculated values of scale are approximately four 
times those predicted by the Dryden relation. Assuming 
values of scale from the Dryden relation, C,; would have a 
value of approximately 0.25. Further work is necessary to 
determine where the truth lies, but one is left with the sug- 
gestion that the flame may be under the influence of turbu- 
lence with an effective scale larger than would be expected 
from the presence of the four-mesh screen. 

Preliminary measurements of wrinkle base width from the 
Fig. 4 photographs yield values in the range of 0.4 to 0.8 in. 
compared with a wrinkle base width of 0.2 in. required for a 
C, = 4. This may be interpreted either as a further argu- 
ment for the presence of a turbulence scale larger than that 
predicted by the Dryden relation, or as evidence, contra- 
dictory to (7), that average base width of the flame wrinkles 
is substantially greater than four times the turbulence scale. 
An increase in effective scale might result either from the 
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divergence of the approach flow associated with the presence 
of the turbulent flame, or it might be attributed to the exist- 
ence of flame-generated disturbances, the scale of which 
would not necessarily be the same as that generated by the 
screen. 


High Effective Turbulence Intensity in Flame Suggests 
Presence of Flame Generated Turbulence 


By modification of Taylor’s (15) theory of eddy diffusion 
by continuous movements to take into account flame propa- 
gation, Scurlock and Grover arrived at the following relation 
for the effect of eddy diffusion on (¥?) 


d(¥?) t 
= ff R.R,dt [2] 
dt 0 
where 
v’ = turbulence intensity (assumed to be isotropic 
and thus the same in all three directions) 
R,and R, = the Lagrangian and Eulerian correlation coeffi- 
cients of fluctuating velocity 
l = the exposure time or age of the flame element 


under consideration 


For short-time intervals, R, and Ry may be assumed to equal 
1, and Equation [2] may be integrated to obtain the following 
simple relationship 


(2)? = v't [3] 


For long-time intervals, R, and R, fall to zero, and the follow- 
ing relation may be obtained 


(¥2)' = [4] 


where isa mixed scale.'° 

Thus assuming a constant v’ in the vicinity of the flame, it 
would be anticipated from Equations [3 and 4] that (¥?)'”? 
would rise linearly with exposure time for a short interval 
and would then increase at a lower rate proportional to the 
square root of exposure time. A more rapid increase in 
(Y?)'/: would indicate that effective turbulence intensity was 
increasing with exposure time. 

Using experimental values of (¥?)'/? and t, Equation [3] 


10 Using Taylor’s definition of Lagrangian scale 


R,ct 
0 
f, Rydu 


and, accepting after Dryden (4) the approximations Rk; = exp — 
v't/l and Ry = exp — Srt/le, the following relation may be ob- 
tained from the mixed scale /,’ 


Equation [3] would apply when t< 7)’. Equation [4] applies 
when t >> 7)’. 


and Eulerian scale 
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Distance above burner rim, in. 
measured (¥?)'/2, in. 
time (corrected for velocity reduction from flow divergence), sec 
v’ (cold flow from Dryden (5)), ft/sec 
v’ (from Eq. [3]), ft/sec 
v’ (from Eq. [4]; Dryden from Table 1), ft/sec 
v’ (from Eq. [4]; /2 from Eq. [1] and C; = 4), ft/sec 


Table 3 Calculated values of v’ 


1.2 1.6 2.0 2.5 3.0 3.5 
.05 0.083 0.123 0.165 0.211 0.276 
.0052 0.0070 0.0088 0.0113 0.0142 0.017 
3 1.2 iad 1.00 0.88 0.80 
.80 0.99 1.16 1.22 1.24 1.31 
22 2.14 3.30 4.30 5.16 7.03 
.52 0.70 1.02 1.32 1.86 2.71 


has been employed for computation of v’. Values of v’ have 
also been calculated from the experimental quantities from 
Equation [4] assuming in one instance the /2 predicted by the 
Dryden relation, and, in a second, that /, has the value ob- 
tained from Equation [1] assuming C; = 4. The results of 
these calculations are shown in Table 3 where they are com- 
pared with values of v’ predicted by the relation of Dryden 
for turbulence intensity downstream of screens. 

Equation [3] might be expected to apply fairly well for 
perhaps an inch or two downstream of the burner rim. Nei- 
ther Equation [3] nor Equation [4], but something in be- 
tween, should apply for the remainder of the distance to the 
flame tip.1! 

Although the Dryden relation predicts a slight decrease 
of v’ with distance downstream due to decay, all of the calcu- 
lated values of effective v’ based on measured values of (¥*)'/”2 
show an increase. It is thus suggested that additional ef- 
fective turbulence is being introduced into the flame over and 
above that generated by the screen. One explanation could 
be the appearance of flame generated turbulence. 


Comparison of Experimentally Determined Curves of 
(S7/Sz) vs. (Sit/l2) With Those Obtained From the Theory 


Scurlock and Grover in their theoretical development em- 
ployed Equation [1], and, attempting to take into account 
the effects on the flames of eddy diffusion, flame propagation 
and flame generated turbulence, derived an approximate 
relation for (Y*)'/? as a function of the dimensionless time 
S1t/le, the ratio v’/S, and the ratio of unburned to burned gas 
density pu/p». Results of their theory presented as (S7/Sz) 
vs. (Szt/lz) for several values of v’/Sz, for density ratios 
of eight and ten, and with and without the presence of flame 


" Taking l, = 0.05 ft, v’ = 1.2 ft per sec, = $l. = 0.025 ft, 
and Sz = 1.2 ft per sec; Jo’ for our experimental conditions is 
found to be 0.014 sec. For our flow velocity of 20 ft per sec, 
this corresponds to a distance of approximately 3.4 in. Flow 
divergence with the attendant reduction in flow velocity might 
reduce this value to about 3 in. 
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Fig. 11 Plots of S7r/Sz vs. Sit/l, under different conditions 
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generated turbulence are shown in Fig. 11. All of the theo- 
retical curves were computed using C; = 4. Superimposed 
for comparison are curves derived from the present experi- 
mental results. The first such curve is for the larger scale of 
turbulence obtained from Equation [1| and the assumption 
that C,; = 4. The second assumes that the effective scale is 
that predicted by the Dryden relation. If, as has already 
been indicated by the experimental results, the correct value 
of C; is much lower than 4, then the first experimental curve 
is the one most directly comparable to the theoretical curves 
(since both are based on C; = 4). Based on the screen 
generated turbulence and the experimental density ratio, the 
theoretical curves most appropriate for comparison are those 
for (v’/S,) = 1 and (pu/p») = 8. The general shape of this 
experimentally derived curve is in reasonably good agree- 
ment with the theoretical. This lends support to the under- 
lying physical picture upon which the theory is based. The 
experimental curve, however, rises well above those obtained 
from theory. The explanation of this deviation could be 
inadequacy of the theory in taking account of flame gener- 
ated turbulence. Additional experimental results under 
other conditions should shed further light on this point and 
prove useful in assigning a more accurate value to the 
constant C;. 


Conclusions 


1 Under the experimental conditions of this study, the 
flame velocity increase in the presence of approach turbulence 
is, within experimental error, entirely accounted for by flame 
wrinkling. Both this result and the appearance of the flash 
tube photographs showing the instantaneous flame front 
(as denoted by the change in concentration of suspended 
particles) to be a thin wrinkled continuous boundary between 
unburned and burned gases, lend support to the wrinkled 
flame concept of turbulent flame propagation. 

2 The turbulent flame velocity is found to increase rap- 
idly with distance above the burner rim, this increase being 
somewhat greater than would be predicted by the Scurlock 
and Grover theory. 

Comparison of theory with experiment indicates the pres- 
ence of effective turbulence in excess of that in the approach 
flow. It is suggested that such additional turbulence might 
be flame generated. There is also evidence that the effective 
scale of the turbulence disturbing the flame is several times 
greater than that generated by the screen. This could be 
explained on the basis either that flame generated disturb- 
ances are present and that these would not necessarily have 
an effective scale equal to that in the approach flow, or that 
the flow divergence in the gases approaching the flame acts 
to increase the scale. Although the experimental values of 
turbulent flame velocity are somewhat greater than predicted 
by the Scurlock-Grover theory, the general shape of the ex- 
perimental (S7/S1z) vs. (Szt/lz) curve is similar to that of 
the theoretical curves. This lends support to the general 
validity of the physical picture of the mechanism of turbulent 
flame propagation which underlies the theory. The experi- 
mental results suggest that the value of 4 previously assigned 
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to the constant C, for calculation of the theoretical curves is 
much too high. 

3 A useful technique has been developed for determi- 
nation for a turbulent flame of the mean flow field in the neigh- 
borhood of the flame, the mean flame front, the rms displace- 
ment of the mean flame front, the local turbulent flame 
velocity, and the position of the instantaneous flame front. 
The analysis of additional experimental results obtained over 
a range of conditions using this technique can be expected to 
provide further enlightenment regarding the mechanism of 
propagation of turbulent flames. 
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On-Off Control System for Attitude 


Stabilization of a Space Vehicle 
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This paper develops certain nonlinear analysis and synthesis techniques as applied to an on-off 


control system for attitude stabilization of a space vehicle. As an example of the possible applica- 
tions of the theory, a hydrogen peroxide system is discussed. However, the analysis holds equally 
true for other gas systems, such as nitrogen, hydrazine, ete. It is shown how the analysis effort is 
simplified when vehicle attitude and attitude rate are represented by a set of nondimensional posi- 
tion-velocity coordinates defining a ‘“‘phase plane.’’ The requirement for compensation (lead) 
networks to achieve stable limit cycle operation is established, and the effect of time delays on 
system performance is discussed. Since fuel consumption is shown to be dependent on lead net- 
work configuration, a synthesis method for optimization of the compensation network to result in 


minimum fuel consumption is presented. 


UTSIDE of the Earth’s atmosphere, the motion of a 

vehicle in coasting flight about its pitch, yaw or roll 

axes is determined by the respective angular rates existing at 

shutoff of the propulsion and control systems which put it 

into space; motion, about its body axes, of a vehicle in pow- 

ered flight could be the result of misalignments of the propul- 
sion system with respect to these axes. 

Because of this motion, attitude stabilization and control 
of a satellite, or a space vehicle in general, may be required 
in some instances for short periods of time. These require- 
ments, for example, may exist because of limitations on camera 
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and antenna look angles or in gimbal travel of inertial plat- 
forms used for navigational purposes. 

There are several known means of achieving the required 
results. The designer, for example, may select a propor- 
tional control system in which the applied correction force 
or moment is proportional to the attitude error sensed by some 
measuring device, or he may choose an on-off (‘‘bang-bang’’) 
system in which a correction is applied only when the error 
exceeds a previously established value (dead zone). 

Proportional systems have the disadvantage that they are 
elaborate and heavy because of the use of control vanes, 
gimballed thrust chambers or similar devices. 

On-off systems, on the other hand, weigh appreciably less 
and are simple and reliable because they consist of relatively 
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Fig. 1 Simplified block diagram of an on-off H,O, control system (pitch axis shown only) 


few components. For this application, hydrogen peroxide 
(H,0.) for example, is particularly well suited, although other 
gas systems, such as hydrazine, nitrogen, etc. are equally 
feasible. 

Fig. 1 shows a simplified block diagram of a single channel 
(pitch axis) of a typical H.O, on-off attitude stabilization and 
control system. The attitude reference system senses the 
pitch attitude of the vehicle, compares it with the reference 
attitude and feeds the difference or error signal into a phase- 
sensitive detector. Depending on the direction of the error, 
clockwise or counter-clockwise, the detector feeds the signal 
into the proper amplifier. When the amplifier output signal 
magnitude corresponds to the given dead zone value, the 
relay pulls in and the H,O2 control valve is electrically ener- 
gized, allowing for flow of H.O. to the proper jet motor. In 
the jet motor, the H,O. passes over a suitable catalyst which 
changes the H,O into steam; thrust is developed by expansion 
of the steam through a nozzle. Corrective action stops upon 
de-energization of the control relay when the attitude error 
becomes less than the dead zone value. 

Several analysis techniques for on-off control systems have 
been developed based on adaption of linear theory (see (1)? for 
example). However, it has been found that these methods 
do not lend themselves too well to the present application be- 
cause of the very low limit cycle frequencies required to 
minimize fuel consumption. 

In contrast, the phase plane technique, as will be shown in 
this paper, offers the designer the advantages of simplicity, 
flexibility and adequate accuracy in predicting the perform- 
ance of the system. 


Vehicle Dynamics* 

This paper will discuss single axis vehicle motion only. 
This approach, although not providing a complete solution to 
the coupled equations of motion, permits a straightforward 
evaluation of those system characteristics which are required 
for the control system design. 

Outside the Earth’s atmosphere, the equations of motion of 
a vehicle with respect to any one of its principal axes become 
extremely simple. 


* Numbers in parentheses indicate References at end of paper. 
* A list of definitions of symbols is given at end of paper. 
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During the ‘‘on” phase of the (or No, ete.) 
trol system, when a correcting moment is applied about the 
pitch axis (identical analyses hold true for the yaw and roll 
axes), the vehicle dynamics are given by 

rl 
= constant [1] 

If the initial conditions at time ¢ = ft are assumed to be 
6 = 6(0), where 0(0) > dead zone and 6 = 0, it immediately 
follows that at any time ¢ 


6 = #0) + [2| 


and 
6=al [3 | 


By elimination of time (t) between Equations [2 and 3}, 
Equation [2] may then be written 


= 40) + 


Equation [4] may be conveniently plotted on a position- 
velocity set of coordinates. This representation is called a 
phase plane (2). Inspection of Equation [4] shows that its 
phase plane trace will be a parabola with vertex at 6 = 6(0) 
and a latus rectum equal to 2a. 

During the “off” phase of the control system, the vehicle 
will rotate at an angular pitching velocity equal to that when 
the thrust of the control system was reduced to zero. The 
vehicle dynamics, therefore, are given by 


6 = 6(0) = constant [5 | 


On a phase plane plot, Equation [5] represents a straight 
line parallel to the displacement axis. 


Introduction to Phase Plane Analysis Methods 


If all quantities in Equation [4] are nondimensionalized, 
a considerable simplification in the analysis effort is accom- 
plished. A parabolic template then can be constructed for 
analyzing the effect of system parameter variations. 
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To nondimensionalize all quantities with dimensions con- 
taining time, a convenient time base, t, = 1 see is assumed. 
This defines a nondimensional time scale 7 = ¢/t,. Angular 
position may be nondimensionalized by dividing 9 by at,?. 
Therefore, a dimensionless displacement 2 may be defined 


6 
aly? [6] 
Also 
= at + at? 
dt dr 
6 
[7] 
and 
Q = - [8] 
Qa 


Inspection of Equation [8] shows that the dimensionless 
angular acceleration is of unity (@ = constant) or zero mag- 
aitude. When a correcting torque is applied, 2 may be 
positive or negative. Conventionally, a clockwise trajectory 
is assumed on a phase plane plot. Q then is positive when 
the dimensionless rate Q is increasing in the positive direction, 
and Q is negative when Q increases negatively. 

By use of Equations [6, 7 and 8], Equation [4] can be re- 
written in nondimensional form: 

In left-hand plane where Q= +1 


Q = M0) + Q? [9a] 
In right-hand plane where Q = —1 
Q = — (9b | 


Parabolic templates to satisfy Equations [9] are easily 
constructed. 

System dead zone and relay hysteresis may be shown on the 
phase plane as displacement boundaries at which relay pull-in 
and drop-out in either direction may occur. 

Assuming, for the time being, instantaneous thrust buildup 
and decay, and neglecting any time delays, such as relay 
armature closing time, ete., which may be present in the 
system, it can be easily shown on a phase plane plot (Fig. 2a) 
that, with relay hysteresis, any disturbance would build up 
with time. That is, the system would hunt at an increasing 
amplitude and frequency. Such a system is dynamically 
unstable. If the relay hysteresis is reduced to zero, the sys- 
tem would hunt at constant frequency and at the amplitude 
of the transient displacement (Fig. 2b). However, in any 
physical system, time delays are always present. The assump- 
tion of instantaneous thrust buildup and decay also cannot 
be considered valid. 

To include time delays in the analysis, it is necessary to be 
able to make time measurements on the phase plane plot. 
Along a constant velocity path [Q = Q(0)], the dimensionless 
time 7r(= t/ty) to move from one point Q; to another Q» is 
given by 

— 
T= 
|Q(0) | 

Along a constant acceleration (parabolic) path, the dimen- 
sionless time to travel from one point to another is given by 
the absolute value of the difference in ordinates Q of the start- 
ing and end points (see Appendix for derivation). 


Comparison of Figs. 2a and 2¢ shows the system amplitude 
and frequency to increase more rapidly when a nominal time 
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delay due to relay armature closing and opening time is con- 
sidered. 

Therefore, some means for system stabilization and fuel con- 
servation must be devised. Both of these requirements 
may be accomplished by incorporation of lead compensation 
in the system. 


System Operation With Lead Compensation 


Stable limit cycle operation can be achieved by incorpo- 
ration of compensation (lead) networks in the system. It has 
been found that a single stage lead crcuit, operating on the 
attitude error signal as represented by the amplifier output, 
is generally adequate when the angular acceleration during 
the system ‘“‘on”’ mode is less than approximately | radian per 
sec®. As will be shown later, fuel consumption as measured 
by limit eycle frequency, or more directly by thrust-on time 
per cycle, may be minimized by proper selection of lead cir- 
cuit parameters. 

The Laplace transform for the output /o(s) of a single lead 
circuit in response to an input /;(s) is given by 


| T's) 
+ (7/y)s! 
For a constant acceleration input, the complete solution 
for £)(t) is given by 


= Ei +-(y7 - DE (y l -—e 7 
Y 


whereas for a constant velocity input 


E,(t) = E + (Y 7 [13] 


The input to the lead circuit is proportional to the vehicle 
attitude (provided the amplifier is not saturated), that is 


E; = 


E\(s) = 


Equations [12 and 13], therefore, may be rewritten in 
terms of 6, i.e. 
K Y 
{12a} 


and 


= +0p = 0+ - —€ [13a] 
kK Y 
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where 


6p = value of @ corresponding to the relay drop-out point 
in the uncompensated system 

6p = value of 6 corresponding to the relay pull-in point in the 
uncompensated system 


Equations [12a and 13a] determine, respectively, actual re- 
lay drop-out and pull-in boundaries on either side of the zero 
displacement axis. 

In nondimensional form, Equations [12a and 13a] are 


and 

+Qp = 24+ (7 im) [15] 

For most practical applications, however, the exponential 

(transient) terms, as will be shown, may be neglected. When 


this is possible, the steady-state relay drop-out and pull-in 
boundaries become straight lines. 


Then 
+ Qp = 2 + m2 — [16] 
and 
+ = 24 mOQ 17] 
where 
m= 1) 
and 
T 
7: 


The validity of the above assumption becomes apparent 
when the time interval 7 between successive relay pull-in and 
drop-out or drop-out and pull-in points on the phase plane 
trajectory is considered to be large compared to the time 
constant c of the second break point of the lead circuit. In 
this paper, the effect of the transient terms will be neglected. 
However, if found necessary, the required correction to the 
steady-state pull-in and drop-out boundaries to account for 
the transient terms may be obtained by graphical methods 
developed by the author. 


| PULL a | @ | 
wean | 


\ INITIAL DISTURBANCE \ 


2 wit | | 2 2) 3 


\ RELAY pRop-our | \ \ 
|ASSUMPTIONS: \ 


| 
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AND DECAY. | | SA | 
| 


2. ZERO TIME DELAYS) \ 
Fig. 3. Phase plane pitch trajectory of a lead-stabilized vehicle 


for an initial angular impulse (steady-state limit cycle period: 
14.2 sec) 
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If thrust buildup and decay would be instantaneous and 
system time delays zero, then the intersections between the 
constant acceleration or velocity trajectory and the corre- 
sponding steady-state relay drop-out and pull-in lines would 
completely determine the system limit cycle. Fig. 3 shows a 
complete phase plane trajectory, defining the attitude and 
attitude rate of a hypothetical vehicle as a result of an initial 
angular disturbance impulse. The stabilizing effect of the 
lead circuit is clearly indicated. Steady-state operation is 
reached when successive trajectories fall on top of each other. 
The steady-state limit cycle trajectory is defined by the 
points 8, 9, 10, 11 and 12. 

The analysis method, outlined above, has assumed instan- 
taneous thrust build-up and decay and also zero system time 
delays. However, generally, these are not valid assumptions 
as will be shown. 

Actual thrust build-up and decay characteristics of typical 
motors were determined by means of ground and flight tests 
on the Viking program in the early 1950’s. They indicated 
the thrust buildup and decay essentially to be of exponential 
nature. After relay closure, when the system enters into its 
“‘on’’ phase, the vehicle dynamics, therefore, can be described 
by 

6 = — € [18 | 


If the initial conditions at time ¢ = t) (where f) corresponds 
to time of relay closure) are assumed to be 6 = @(0), and 
6 = 0 at any other time ¢ 


6 = alt — Tl — [19] 


and 
2 
6 = 00) + al + — — ra | [20] 


Inspection of Equations [19 and 20] shows that time ¢ can- 
not any longer be eliminated from the equation for 6, as had 
been done previously (see Eq. [4]). 

Similarly, after opening of the relay armature, when the 
system enters into its “off” phase, the vehicle dynamics are 
described by 


6 = 60) — aT(1 — € [21] 


assuming the thrust build-up and decay time constants to be 
identical. 

Phase plane analyses of systems whose dynamics explicitly 
contain time ¢ are difficult and laborious. However, a con- 
siderable reduction in analysis effort may be accomplished 
when it is assumed that the exponential thrust build-up and 
decay characteristics can be represented by simple time 
delays equal in length to the time constant 7). This assump- 
tion is reasonable, since it is based on the equivalence of total 
impulse and the low frequency power spectrum distribution 
of the actual and assumed thrust characteristics. 

By experiment, system designs based on phase plane anal- 
yses with assumed simple time delay thrust characteristics 
have been found to be conservative as far as actual fuel con- 
sumption and damping is concerned. Experimental data 
show that full thrust is generally not attained during steady- 
state operation; that is, relay drop-out and, asa consequence, 
fuel flow stoppage occur prior to achievement of full thrust. 
Actual fuel consumption, therefore, is somewhat less than 
predicted by the phase plane analysis method. 

Experiments have also shown that there may be periods 
when simultaneous operation of opposing thrust motors 
could oceur as a result of the exponential nature of the thrust 
buildup and decay. This results in a reduction of attitude 
rate; hence, transients are damped out somewhat faster than 
predicted by the phase plane analysis method. 

The phase plane analyses must also consider the effect of 
system time delays. The more significant ones include those 
resulting from relay armature motion prior to actual relay 
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pull-in or drop-out. All of these time delays may be included 
in an equivalent time delay T'y. 
Then, the total system time delay becomes 


To = To + Tu [22] 
or, in nondimensional form 
Tp = + Tu [23] 


Equations [16 and 17] which defined the steady-state 
relay drop-out and pull-in lines respectively, can now be re- 
placed by thrust-off and thrust-on lines, i.e. 


+ Qp = 24+ (m — Ec + Tp) - Q [24] 
and 
+ Qp = Q4+ (m [25] 


It is obvious that the thrust-on and thrust-off lines have a 
different slope than the relay pull-in and drop-out lines. 

Fig. 4 shows a complete phase plane trajectory determined 
with the aid of thrust-off and thrust-on lines and with an 
assumed total system time delay of 7/5 sec. The vehicle and 
lead circuit configuration as well as the initial conditions are 
identical to those used for the construction of Fig. 3 which 
demonstrated the use of relay pull-in and drop-out lines. 


Selection of Lead Network Parameters 


It has been shown that a stable limit cycle can be achieved 
by the insertion of a lead network into the control system. 
It is then quite apparent that the limit cycle characteristics 
strongly depend on the lead circuit parameters. On the 
other hand, total fuel consumption is a direct function of the 
limit cycle frequency and the thrust-on time per limit cycle. 
Therefore, to minimize fuel consumption, the limit cycle 
period should be as long as possible and the thrust-on time 
per limit cycle should be as short as possible. 

It can be easily shown that the nondimensional thrust-on 
time per cycle Ton is given by 


Ton = 4 lQz | [26] 


where Q7 = maximum positive or negative angular velocity 
during limit cycle operation (nondimensional). 
Similarly, the nondimensional thrust-off time per cycle Tors 
is given by 
Q,/ 
Tott = 4 pa [27] 


where Q, = angular displacement at end of each constant 
velocity portion of limit cycle (nondimensional). 
In the right-hand side of the phase plane plot 


Qe + + To) + (T?/2) 


9 


Qs [28] 
and the maximum positive angular velocity during limit cycle 
operation is given by 


Qe — (Qn + me + Tp) —(Tp?/2) 
2(m — Tp) [29] 


But Qp may be expressed as a percentage of Qp. (This is 
true because Qp is a direct function of relay hysteresis and 
Qp.) 

Therefore 


r= — [30] 
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THRUST-ON ASSUMPTION: To =.1 
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Fig. 4 Phase plane pitch trajectory of a lead-stabilized vehicle 
for an initial angular impulse (steady-state limit cycle period: 


11.05 sec) 
Then 
m— Tb m(c + Tp) —(Tp?/2) 
[31] 
and 
| Qr(1 + 
m(c + Tp) — (Tp?/2) ] (32) 


Qp(1 — r) 
+ Tp) —(Tp?/2) 


Tott = 4(m — Tp) [ 


or, by letting 


Tp 
m(c + Tp) — ¢, To) = f 


and 
(m — Tp) = g(m, To) = 9 
— r 
and 
+ 7) 
f 
= ¢ 2: 
Toff 4g =?) [3 
f 
Also 
1 
| [33] 
f 


Figs. 5, 6and 7 show the significance of Equations [32a] and 
[33]. 

Inspection of Figs. 5 and 7 indicates that the ratio of Ton to 
Tort May be minimized (except for r equal to unity, i.e., zero 
relay hysteresis) by proper selection of lead circuit param- 
eters. This minimum is given by 


(=) — r)? 
Tott/min [(T/¥)(¥ — 1) — to]? +17) (8 — 


[34] 
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Generally, Qp, r, and rp are known. In most cases, there 
is also an upper limit on the value of y, the lead circuit atten- 
uation factor, because of signal-to-noise ratio considerations in 
the design of electronic hardware. In a typical design, the 
maximum value of y is usually taken to be equal to 10. 

By use of the graphs or of Equation [34], it is therefore 
possible to minimize fuel consumption by selection of an opti- 
mum value for the lead circuit time constant 7. Or, for a 
given set of system parameters, the limit cycle characteristics 
may be obtained directly without the plotting of phase plane 
trajectories or the use of other techniques. 


Appendix Time Measurements Along a Constant 
Acceleration Path on a Phase Plane Plot 


The equation of the parabolic phase plane trajectory de- 
termines the dimensionless time 7 required to move from one 


Fig. AS a function of circuit parameters 
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point (2, Q;) to another (Qe, 
The equation of the parabolic path is given by Equations 


[9a and 9b] or 
2 
2\dr 


Therefore, when the acceleration is negative 


dt 1 


A-? | 

— Q| 
or when the acceleration is positive 

dt 


— 2(0)] 

Since the quantities under the radical signs in Equations 
[A-2 and A-3] are always positive, the time 7 to go from Q(()) 
to another point, (%, 2:), may be written as 

1 dQ 


T=— 
v2 2) V |2 — 20) | 
or 
2[V/2 20) |] 


But |Q — Q(0)| is equal to 92/2. Therefore 


= 


= | [A-5 | 


Since the time 7 required to move from (Q), Q;) to (Qs, 2.) 
must be positive and is equal to the difference in time to go 
from Q(0) to the points in question 


TQ) = Q, = = [A-6 
Also, since 
t t (1 ) 
Te-=- = 1 sec 
1 
tam, = |AQ| (see) [A-7] 


Fig. 7 ton/Tott VS- Toft 
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Nomencla 


= 
6, 4 = 
0), 00) = 


f(t) = 
= 
BE; = 
Es) 


ture 


pitch attitude error, rad 

first and second time derivatives of 6, respectively 

initial pitch attitude and rate, respectively 

Fi/I, angular acceleration about the pitch axis, 
rad /sec? 

applied control (correcting) force, lb 

moment arm, ft 

pitch moment of inertia, slug-ft? 

time, sec 

time of relay pull-in or drop-out 

1, time base, sec 

t/ty, nondimensional time scale 

nondimensional angular displacement; nondimen- 
sional form of 6 

nondimensional angular velocity and acceleration, 
respectively 

nondimensional form of 6(0); vertex of phase plane 
parabola or Q-axis intercept 

nondimensional form of 6(0) 

Laplace operator 

lead circuit output 

Laplace transform of 

lead circuit input 

Laplace transform of FE; 


Y 

Ei, E; 
K 


time constant of first break point of lead circuit, sec 
lead network attenuation factor, nondimensional 
first and second time derivatives of E;, respectively 
proportionality constant 


Op = value of 6 corresponding to relay drop-out point in 
the uncompensated system, rad 

Op = value of @ corresponding to relay pull-in point in the 
uncompensated system, rad 

Qp = nondimensional form of @p 

Qp = nondimensional form of 6p 

c = T/y, time constant of second break point of lead 
circuit, sec 

m = (T/y) (7 — 1) 

0 = thrust buildup and decay time constant, see 

Tp = total system time delay, sec 

Tw = equivalent system time delay due to other causes 
than thrust buildup or decay, sec 

Tp = nondimensional form of 7'p 

™ = nondimensional form of 7's 

Tie = total thrust-on time during each limit cycle, non- 

dimensional 

Qy = maximum positive or negative angular velocity dur- 
ing limit cycle operation, nondimensional 

Toff = total thrust-off time during each limit cycle, non- 
dimensional 

Qs = angular displacement at end of each constant ve- 
locity portion of limit evcle, nondimensional 

r = relay hysteresis coefficient =(@)/@p) = (Qp/2p) 

= me + rp) (rp?/2) 

g =m-—tp 
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The use of ionization gages and omegatrons in rockets and satellites at very high altitudes is con- 
sidered. Appropriate instrumentation employing these devices is believed to offer good possibilities 
for the measurement of atmospheric pressure, temperature, density and composition. 


HE MEASUREMENT of atmospheric pressure, tem- 

perature and density, in regions of the high atmosphere 
where the mean free path of typical particles is long compared 
with the dimensions of the instrument, has interested upper- 
air researchers for many years. Early attempts (1 to 4)* using 
V-2 rockets were largely inconclusive because of instrument 
contamination by rocket gas or because the altitudes attained 
were insufficient. Some significant measurements to 120 km 
were made, but could not be considered as entirely fulfilling 
desired conditions of a long mean free path environment. 
No further significant measurements were made until the 
Viking rocket, and still later the Aerobee-Hi rocket, became 
available (5 to 9). 

The methods employed in these measurements have made 
use of local density determining devices which permit in- 
vestigation of the pressure or density at selected points on the 
surface of a rocket. This allows interpretation in terms of 
ambient conditions, provided there is knowledge of rocket 
velocity, altitude, air composition and other such factors. 

Advances in instrumentation, the availability of rockets 
with an increased capability for high altitudes, and the im- 
minent availability of satellite vehicles have stirred a com- 
pelling new interest in these measurement areas. 

Thus this paper considers the employment of ionization- 
type devices for performing measurements of ambient pres- 
sure, temperature, density and in addition, composition, in the 
light of the newer developments. 


Measurements With a Single Ionization Gage 


Each gas in a mixture of gases exerts a partial pressure upon 
the walls of the confining vessel in proportion to the gas con- 
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centration. Since pressure is defined in terms of the rate of 
momentum transfer, these partial pressures are also deter- 
mined by the mass, velocity and, hence, temperature of the 
various constituents. 

The presence of particles of different mass implies a velocity 
spectrum and thus suggests that some form of mass spec- 
trometry approach is applicable to the measurement of struc- 
tural parameters at altitudes where the mean free paths of 
various particles are large in comparison with the dimensions 
of an intruding object. 

Consider a body moving in space at velocities comparable to 
the random molecular speeds of the constituent gases. As- 
sume also that the body’s velocity and orientation are known. 
If an orifice to a body-borne chamber is exposed to the out- 
side, an interchange of particles will take place in a manner 
that is dependent upon the particle velocity distribution func- 
tion, the body velocity, its orientation, the characteristic tem- 
perature, the composition and other effects. 

If the body velocity V, is reduced to zero, then the inter- 
change will be altered only in degree. The relationship be- 
tween the internal and external pressure will be governed by 
the thermal transpiration effect as investigated by Maxwell 
(10), Reynolds (11) and Knudsen (12). If the body velocity 
is not zero, then the pressure relationship is less easily stated, 
but readily so upon assumption of a velocity distribution func- 
tion. Many experimenters (7,13,14,15) have derived a relation 
which is seen to be the transpiration equation modified by a 
factor dependent upon the ratio of body to particle velocity. 
Thus 


P; 
+ (1 — erf s)] (1] 


where 


f(s) = + sv/m (1 — erfs) 


P; = internal pressure 

P, = ambient pressure 

7; = internal temperature 
T, = external temperature 
s =V./V, 


component of body velocity normal to gage opening 


most probable particle velocity 


P; was defined only as the pressure in the chamber, implying 
total pressure. It can be considered as the sum of the partial 
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pressures of the component gases; thus 


P; = p(N2) + p(N) + p(Oz) + p(O) + p(A) + (2) 
p(He) + p(H) + p(He) +... 


considering some of the possibilities, including dissociation. 
Similarly, P, may be considered as the sum of the external 
partial pressures, and accordingly an equation like [1] can be 
written for each component. 


Since 
V, = 
Mo 


k = gas constant per molecule, 1.381 X 10~!* erg per deg K 
<i temperature, K 
Mo = mass per particle, gm 


where 


an f(s) exists for each component. Thus, the internal partial 
pressure ratios will differ from the external ratios, the dif- 
jerence being a function of the body velocity component 
normal to the chamber opening. 

Thermal equilibrium among gas particles can be assumed 
externally (16) as well as internally with small error. In- 
ternally, however, it is assumed that the particle tempera- 
tures will be set by the wall temperature, adopting a unity 
accommodation coefficient.? A Maxwellian velocity distribu- 
tion in the region being sampled follows the equilibrium as- 
sumption. 

These considerations suggest the possibility of composition 
measurement by observation of the change in pressure in a 
chamber as a function of Vz, the body velocity component 
normal to the chamber opening. As the relative velocity in- 
creases (or decreases), the total pressure will, in principle at 
least, increase (or decrease) gas by gas, so to speak, as re- 
quired by the above considerations. The rate of change will 
be governed by the distribution function as well as by the 
normal acceleration component which presumably can be 
known. The maximum relative velocity must appreciably 
exceed the velocity of most of the particles to make this 
measurement possible. Expected relative values of the most 
probable external particle velocity corresponding to the 
various components are readily determined upon assumption 
of external temperature and temperature equilibrium. Typi- 
cal values are presented in Table 1, which has been computed 
from Equation [3] for three temperatures. Since m can vary 
relatively from 1 (atomic hydrogen) to as much as perhaps 
40 (argon), V, can have a range of about 6.5 to 1 for a given 
temperature. For purposes of evaluating the feasibility of the 
experimental approach, the minimum expected particle ve- 
locity can be taken as 1300 fps (argon at approximately 400 K), 


7 See (15) for discussion of exceptions to this assumption. 


and the maximum as 15,000 fps (hydrogen atom at 1200 K) 
for any particle except the electron. Considering both rockets 
and satellites as possible research vehicles, satellite velocities 
are adequate for all particles, whereas rocket velocities of 
present upper-air experimentation vehicles are adequate only 
for the more massive particles. Thus, in general, this pro- 
cedure is discussed from the point of view of an orbiting 
vehicle. 

Consider application to a satellite, where the velocity may 
be taken as 26,000 fps. Thus, if the chamber opening is op- 
posite to the direction of motion, only a few of the least mas- 
sive components, as determined by the distribution function, 
could enter. On the other hand, if the opening is oriented in 
the forward direction, then all particles can be overtaken, 
again excepting a few of the lighter components. It follows 
that particular intermediate positions enable partial collection 
of the heavier atoms. 

In the expected situation of rotation that would cyclically 
orient the opening in the direction of motion and then opposite 
to the direction of motion, velocity “scanning” would occur. 
The measured total pressure would, in principle, vary in 
accordance with the concentration of the various constituents 
and their associated velocities, thus offering the possibility of 
a particle velocity measurement. 

The s function (Fig. 1) effect upon the internal partial pres- 
sure of a gas as compared with the external pressure is to re- 
duce the internal pressure to zero as s approaches —1 or to in- 
crease the internal pressure significantly as s assumes larger 
positive values. Since s is greater for the more massive par- 
ticles, contributions to total pressure from No, N, O2and O can, 
under high velocity conditions, tend to overwhelm the lighter 
and less plentiful components out of proportion to their ambi- 
ent composition stature. For negative velocities, the effect is 
opposite, possibly enabling detection of the lighter compo- 
nents. This is illustrated in Fig. 2, which shows the ratios of 
internal to external partial pressures of nitrogen and hydrogen 
as functions of the angle of attack which is a measure of the 
normal velocity component. The curves are valid only for 
equal concentrations of the gases, a condition assumed for the 
sake of illustration. Some of the following data were used in 
computing the curves: 


Altitude: 300 km 
Pressure: 2.0 X 10°* mb 
Density: 5 X 10° kg/m? 
Temperature: 1000 K 
Body velocity: 26,000 fps 
Rotation: 1 rps with rotational axis perpendicular to plane 
of orbit 
Internal temperature: 300 K 
Most probable velocity: 
atomic nitrogen at 1000 K = 3. 
atomic hydrogen at 1000 K = 1 


56 X 10° fps 
.33 X 104 fps 


Table 1 Most probable velocity of various atmospheric components as a function of the kinetic temperature 
——300 K——— 600 K -1200 K— 
Component em/sec fps em/see fps cm/sec fps 

air 4.17 X 104 1.36 x 105 5.9 X 104 1.93 X 103 8.34 X 10! 2.75 X 10? 

N2 4.2 1.38 5.95 1.95 8.4 220 

N 5.95 1.95 8.4 rr if 11.9 3.9 

Op 3.95 1.3 5.58 1.83 8.0 2.62 

O 5.56 1.82 7.9 2.6 1.11 X 105 3.65 

He 1.12 X 105 5 a 1.58 X 105 5.2 2.24 7.39 

H 2.23 7.3 3.15 1.04 x 10! 4.46 1.47 X 104 

H, 1.58 5.2 2.23 7.3 X 10° 3.01 9.9 

A 3.52 X 104 1.16 5:0 XX 10* 1.65 7.04 X 104 2.3 

e 5.18 X 107 1.7 X 106 7.3 X 10° 2.4 X 10 1.03 X 108 3.4 X 108 
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Fig. 1 Curve of f(s) vs. s 


Because the chamber wall temperature is likely to be sub- 
stantially less than the external temperatures, the internal 
pressures can be less by as much as a factor of 2. However, 
the reduction will take place for all components and thus has 
a bearing on sensor sensitivity considerations rather than on 
relative partial pressures. 

Although velocity scanning occurs as a result of rotational 
motion and in principle enables differentiation between the 
various component gases, as a practical matter there is a 
considerable overlap of the distribution functions primarily 
because the region is populated with gases of comparable 
mass. As a consequence, it is probably not possible to deter- 
mine relative concentration by this method. It may be 
feasible, however, to determine the most probable velocity and 
hence the characteristic temperature. If the assumption of 
external thermal equilibrium is valid, or nearly so, the charac- 
teristic temperatures of the various gases should be the same. 
Thus determination of the most probable velocity for any 
component indicates the external temperature. Knowing the 
temperature, it is then possible to return to Equation [1] and 
compute the external pressure. 

Summarizing, a pressure (or density) measuring device, such 
as an ionization gage in a body which is moving at an appro- 
priate and known altitude and velocity, with a prescribed rota- 
tional rate, makes pressure and density measurements possi- 
ble. Preferably, the maximum velocity should exceed, by an 
adequate margin, say 2:1, the most probable velocity of the 
lightest component to be measured, thus indicating the need 
for a vehicle of satellite capabilities, if the lighter gases are to 
be considered. 

With such an arrangement, it appears that it may be feasible 
to measure the temperature directly, which would enable a 
check of the directly measured ambient pressure and density. 


Ionization Gage Plus Omegatron 


The addition of another instrument which permits an ad- 
ditional and independent measurement of the chamber gas 
offers considerably greater promise of fruitful experimenta- 
tion. Assuming that a small mass spectrometer-type device 
is placed within the ionization chamber, so that the same 
gases measured by the ion gage are measured by the spec- 
trometer. A less than ideal but more practical arrangement in 
which the spectrometer would be independently mounted 
within the same sphere is also suitable. An appropriate in- 
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Fig. 2 Ratio of external to internal pressure for nitrogen and 
hydrogen as a function of angle of attack 


strument is the omegatron, a device reported by Sommer, 
Thomas and Hipple (17) who were concerned primarily with 
measurements of charge to mass ratio. A version of the de- 
vice applicable in the sense required here was later developed 
and reported by Alpert and Buritz (18) who suggested its 
usefulness as a pressure-measurement tool.’ A “Diatron’’ is 
also suggested.® 

The addition of such a device, which appears quite feasible, 
would enable determination of (a) the composition, (b) the 
shape of the distribution function and (c) the relative con- 
centration of the constituent gases, provided the associated 
sensor had adequate sensitivity to detect the less plentiful 
gases. Thus, the “total” pressure determined by the ion gage 
could be resolved in terms of the various partial pressures. 

Considering Equation [1] in partial pressure form 


Ip. 

Pai = Poo f(s) [4] 
where 
f(s) = + svn (1 + erf s) [5] 

s = [6] 

V.= 2k 17 

Ma 

P.; = internal partial pressure of gas a 


Pao = external partial pressure of gas a 


and 


s=V, 
Substituting in Equation [4] 
Pai = Pao (vag) j 
rs V, 1 4 erf V, | ma ) 
KT, V x7. 


8 Possible usefulness at high altitudes was suggested by H. 8. 
Sicinski. 
® Consolidated Electrodynamics Corp. 


[8] 
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This is the general form for any value of V; and hence any 
value of s. In simplification, as 


s—>2, erf—> 1 


and 


Ma 
(s) & 3.55V, 
Hie) 


Thus Equation [8] becomes 


Pai = Poo (3.55 (9] 
AKT, 


which, by combining, becomes 


Pat = 3.55 (10] 


It is considered here that the only unknowns are P,,. and 
7, which appear as a ratio. Consequently, since the pressure 
can also be defined as 


Poo = tiakT [11] 
where na = number of particles of gas a, one can substitute 


2 
I av 


— 
[12] 


nak = 


Equation [10] becomes an expression for the number density 
in space of the gas a as a function of the internal partial pres- 
sure, the normal component of the body velocity, and the in- 
ternal chamber temperature. Thus 

= 
V0.58kmaT: V2 

It follows that the ambient number density of any constitu- 
ent can be determined from the internal partial pressure. In 
addition, knowledge of the external temperature, so obtained 
from a characteristic velocity measurement with an ion gage 
or similarly with an omegatron, will permit computation of 
the external partial pressures. 

Thus through employment of an ion gage and an omegatron, 
which admittedly involves redundancy, quantities can be 
evaluated, including: Pressure, temperature and density of 
ambient gas; type and concentration of the constituent 


gases. 


Velocity Distribution Function 


The preceding comments have been based upon the as- 
sumption of a Maxwellian distribution of velocities among the 
various gas particles. Although there is considerable evidence 
that this is true or very nearly so (15), it is of course desirable to 
obtain experimental verification, which the presence of the 
omegatron, or a similar instrument, makes possible. 

Consider again the omegatron in the body whose orienta- 
tion pattern is such as to bring the chamber opening cyclically 
into and away from the “stream.” Assuming diffusive equi- 
librium between the chamber and the outside, the internal 
number density of a particular gas, which is the quantity 
measured by the omegatron, would vary as a function of the 
normal velocity component and the distribution function. 
Since the body velocity is known, determination of the dis- 
tribution function becomes possible. 

A plot of the change in internal number density (or partial 
pressure) as a function of angle of attack, or normal velocity 
component, should appear like Fig. 2, provided, of course, 
that the distribution is Maxwellian. Other distributions 
would alter the shape and would thus be definable. 
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Atmospheric Velocity 


All previous considerations have assumed that the relative 
velocity between the body or chamber and the gas was due to 
(a) that characteristic of the gas and (b) that of the instru- 
ment-carrying body itself. Accordingly, the assumption is 
implied that the atmospheric velocity is negligible. This is 
not the case, for if one considers that the atmosphere moves 
with the Earth, the velocity in the equatorial plane at an alti- 
tude of 300 km is 1.6 X 10 fps. 

This value is clearly significant in comparison with the most 
probable velocity of certain components (see Table 1). The 
effective value, insofar as possible measurements are con- 
cerned, is of course dependent upon the location of the plane 
of the orbit, and thus is maximum and nearly constant for 
passage near the equator and variable in magnitude and direc- 
tion for polar orbits. In any case it can be considered as a 
definable perturbation on the body velocity. 


Sensor Sensitivity 


To employ fruitfully the above mentioned devices, it is of 
course necessary to arrange a sensor that will ideally have 
adequate sensitivity and time characteristics so as to permit 
detection and transmission to the recording device of signals 
corresponding to the least plentiful components. 

The major constituents one expects at the higher altitudes 
are NO, Oo, O, No and N, in both neutral and ion states. 
Speculation leads to consideration of hydrogen and helium as 
well, but in very small quantities. Detection of the various 
gases ‘s dependent primarily upon their absolute abundance, 
and upon the sensitivity of the current-measuring instrument 
employed. D’fferent elements in like concentration effect dif- 
ferent responses in an ionization gage. In general, the ap- 
parent sensitivity is a function of the total number of electrons 
possessed by the molecule (19) rather than its ionization po- 
tential or atomic weight, and thus the ratio of sensitivities of 
helium and nitrogen is less than their mass ratio. The overall 
result, then, insofar as ion gage (or omegatron) sensitivity is 
concerned, is that the relative concentration is the more 
powerful factor. 

Accordingly, in regard to gage sensitivity (S) to dissociated 
gases, one would expect the ion current resulting from disso- 
ciated or recombined gases to vary with the number of elec- 
trons per particle. 

Ion gage sensitivities can be considered in terms of the ratio 
of collected ion current to ionization current at a given 
density. A typical value for present thermionic gage capabil- 
ity is S = 100 wa per microbar for an ionization current of 10 
ma. Thus at this pressure each 100 electrons, on the averrge, 
emitted by the cathode produces one ion for collection. The 
density is such that it is unlikely one electron can generate 
more than one ion pair during transit. Thus the sensitivity 
factor indicates that an electron has, on the average, one 
chance in a hundred of producing an ion pair at one microbar 
of pressure. 

Further, as long as the electron has sufficient energy to over- 
come the particle ionization potential, the energy level for 
typical thermionically-produced electrons is not of great con- 
sequence. 

Similar sensitivity factors are appropriate for at least some 
radioactive ion gages, such as a radioactive ionization gage 
using a tritium source (20). However, the resulting ion cur- 
rents are small compared with a thermionic gage, because the 
number of ionizing source electrons is very small. This factor, 
a limit imposed by practical considerations in the use of radio- 
active ionizing sources, in part precludes use of radioactive 
ion gages at very low densities. An additional factor that 
limits the usefulness of these devices is a dark current which 
overwhelms the ion current at densities that are relatively 
high in relation to this report. The dark current is a mani- 
festation of collector-electron emission resulting from x-ray 
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emission caused by the high energy radioactive source elec- 
trons, ions or nuclei, as the case may be. 

Phillips cold-cathode ionization gages are not useful at pres- 
sures less than 10~* to 10~7 mb Hg for the reason that the 
gaseous conduction which underlies their operation cannot 
be sustained at the lower density values. Thus the thermionic 
ionization gage appears to be the only useful type at pressures 
to be encountered at very high altitudes. 

At 300 km, and with an assumed pressure of 2 X 10~* mb, 
the ion current, assuming ground air composition, can be de- 
termined from 


10-3? SP [14] 
where 
7 = amp 
S = sensitivity in wa per microbar 
P = pressure in millibars 


Thus, for this example one should expect 


100 X 2 X 10° X 1078 
2 X 10-* amp 


Currents of this magnitude do not pose difficult problems 
for instrumentation design, for currents as low as 10~!? amp 
are presently being successfully measured in upper atmosphere 
research program experiments. 

The measurement of the ion current produced by an omega- 
tron poses a somewhat more difficult problem. Although 
omegatron and ion gage sensitivities are comparable, the use 
of a relatively small ionizing beam, rather than the gross 
emission from a relatively long filament, as employed in ion 
gages, results in greatly reduced ion current. Generally speak- 
ing, the extent to which the beam current and resulting beam 
size may be increased is limited because of the difficulty of 
extracting the innermost ions from the beam. Excessive space 
charge is also a problem. 

In addition to straightforward ion-current measurement, 
there are at least two other ways in which the low level output 
signal of the omegatron can be satisfactorily utilized. One of 
these involves use of an electron multiplier at the collector 
which can provide considerable amplification of the desired 
signal. The other method, also promising, involves detection 
of the resonant ions by measurement of the energy absorbed 
from the applied RF field (21). 

Accordingly, because of these and other sensitive detecting 
devices certain to be developed, it is reasonable to conclude 
that the thermionic ion gage and the omegatron have ade- 
quate sensitive for use in measurements of the type discussed 
in this paper. 


Difficulties To Be Anticipated in Use of 
Thermionic Ion Gages and Small Spectrometers 


The problem of measuring very low pressures at very high 
altitudes is not different in principle from performing the 
measurements in the laboratory, with the major exception of 
unattended operation. Thus perturbations caused by out- 
gassing of gage elements and the vehicle, pumping action by 
gage cathode, adsorption and absorption by gage elements, 
etc., must be considered. As the surfaces of the carrier be- 
come heated by solar radiation, they will produce quantities 
of gas which can perturb the measurement. Similarly, as they 
cool, gas will be sorbed. Data presented by Russian ex- 
perimenters as measured on Sputnik ITI are of interest in this 
regard. Ata recent conference” it was stated that the mini- 


20 CSAGI, Moscow, Aug. 1958. 


mum pressure measured the first day in orbit was 10~? mm 
Hg, the second day, 10~* mm and the third, 10-® mm. 

The hot cathode of an ion gage or omegatron is a major 
offender in regard to absorption and particularly in regard to 
chemical reactions with the gas being measured. Thus steps 
to reduce cathode temperature and perhaps ‘‘on time” would 
be desirable from this standpoint. Specifically, it appears de- 
sirable to employ oxide-coated rather than pure tungsten 
cathodes to enable lower temperatures and to conserve power. 
The potential assumed by the vehicle can likewise con- 
ceivably cause perturbations upon the assumed environment. 
Some studies bearing on this problem are being made at pres- 
ent in connection with the application of Langmuir probe 
techniques to the measurement of various ionosphere char- 
acteristics. 


Concluding Remarks 


This paper has presented a brief discussion of a measurc- 
ment technique which seems to the authors to offer reasonab) 
assurance of satisfactory measurements of pressure, tempera- 
ture, density and composition at very high altitudes. Con- 
siderable development work leading first to use in vertica! 
rocket soundings will of course be necessary prior to satellite 
use; however, the apparent promise of significant data will 
provide sufficient incentive for the development effort  re- 
quired. 
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Technical Notes 


Cireular Orbit Stability in the Oblate 
Earth’s Equatorial Plane 


MARVIN H. HEWITT! 
The Martin Company, Denver, Colo. 


In the oblate Earth’s field, circular equatorial orbits are 
stable for small planar perturbations. Disturbances caus- 
ing orbital precession are not considered here since their 
effects require a consideration of the Earth’s noncentral 
force field, which is outside the scope of the method used in 
this paper. 


'PYHE POTENTIAL function for the gravitational attrac- 
tion of a solid homogeneous oblate spheroid upon a distant 
unit particle is, approximately? 


MoM 


r 


2 — p?) 
P;(cos +... (1) 


r3 
where 


P2 (cos 6) = Legendre coefficient 


M = mass of the Earth 

r = distance from center of Earth to satellite 
a = semimajor axis of the Earth 

b = semiminor axis of the Earth 


In the equatorial plane the potential function reduces in 
approximation to 


V=- 
r 10 


[2] 


The corresponding acceleration g is given by the expression 


M 3M (a? — b?) 


[3] 


It is seen that this acceleration represents a central non- 
Newtonian foree, and it is well-known that a circular orbit is 
always possible for central interaction. The following sta- 
bility analysis will be restricted in application to perturbed 
motion confined to the equatorial plane. We define stability 
here in the sense of linearized perturbation theory. Perturba- 
tions that cause orbital precession will not be considered here. 
Satellite motion in a plane inclined to the equator involves 
a consideration of the Earth’s noncentral force field which 
is outside the scope of the method used in this paper. 

Let u = 1/r, where r is the orbital radius. Then the at- 
traction on unit mass is a function of u, and we may write it in 
the form 


P = ku®f(u) [4] 


Received Jan. 15, 1959. 
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third ed., chapter VIII. 


where 
k = —GMm [5] 
3 
=1+ 10 10 (a? — b?)u? [6] 


m being the unit mass. 
The differential equation of the orbit is 


k 
dg? tum [7] 


where h denotes, as usual, the constant moment of the velocity 
about the center of attraction, and @¢ is the orbital polar angle. 
Let c be a constant given by the equation 


k 
¢ = pet) [8] 
and note that vw = ¢ satisfies Equation [7]. The orbit is now 
perturbed, and we set 


u=c+r [9] 
where x is initially small. 
Substituting c + x for u in Equation [7], we get 


@x 


[10] 


For small x, we can expand the right-hand side of Equation 
[10] by Taylor’s theorem and delete all powers of x higher than 
the first. Thus 


dr k 
[f(e) + [11] 
Since ¢ satisfies Equation [8], this last equation gives 
d*x 
ae TO 
[12] 
af'(c) 
or 
d?x 
+2[ [13] 


If the coefficient of x is positive in Equation [13], this takes 
the form 


dd? + a’r = 0 [14] 
the solution of which is 
x = A cos (aod + B) [15] 


Therefore, if A is small, x is always small, and the orbit is 
therefore nearly circular. But if the coefficient of x in Equa- 
tion [13] is positive, the equation has the form 


— —}?r=0 [16] 
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the solution of which is 


xz = Ce? + [17] 


and thus z is large when ¢ is large (unless C is exactly zero) 
until the approximate Equation [13] is no longer applicable. 

The conditions, therefore, that a satellite in a circular orbit 
should be stable, undetermined, or unstable are, respectively 


> 
f(c) 


Substituting parametrically, we obtain Equation [18] in the 
form 


[18] 


[(6/10) (a? — b*)c?] 


~ + G/10) (a? — bye] [19] 


where satellite stability exclusively is realized. In the case 
where the equality holds, the corresponding orbital radius is 
given as 


— 
ae aa = 134 nautical miles = r [20] 


We see that this orbit lies inside the Earth. And it is ap- 
parent that the unstable case requires a smaller orbit than 
that given by Equation [20]. 

Inequality [19] holds then for actual satellite orbits; we 
conclude that perturbed satellite orbits in the equatorial 
plane are stable for small perturbations in that plane. 
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Particle Impaction on Spheres at High 
Mach Numbers! 


S. K. FRIEDLANDER? 


The Johns Hopkins University, Baltimore, Md. 


At high Mach numbers a detached shock wave is estab- 
lished upstream of the forward stagnation point of a blunt 
body. The equations of motion of a small particle in such 
a flow are set up, and deposition at the forward end by im- 
paction is considered. An expression is derived for the im- 
paction parameter, defined as in subsonic flow, at which 
the efficiency vanishes. This is applied to the high Mach 
number flow around a sphere, and it is found that the criti- 
cal impaction parameter varies from 0.153 to 0.29 for den- 
sity ratios from 4 to 14. A sample calculation is given. 
Possible applications are to high altitude sampling of par- 
ticulates, to prediction of erosion and to the development 
of more efficient laboratory samplers. 


HEN a gas flows at a high Mach number over a blunt 
nosed body, a detached shock (“bow”) wave is estab- 
lished upstream of the forward stagnation point. Between 
the shock and the stagnation point, the gas density is high 
and the flow subsonic. In this paper, the equations of motion 
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of a small particle in such a flow are set up, and deposition at 
the forward end by impaction is considered. These calcula- 
tions may have application to high altitude sampling, to pre- 
dictions of erosion and to the development of more efficient 
laboratory samplers. 

The calculation of the complete impaction efficiency curv: 
would require a numerical solution to the equation of particle 
motion with varying drag coefficient in a manner similar t: 
that for subsonic impaction (1). At this point the uncer- 
tainty in the velocity distribution at high Mach number; 
would not warrant such an effort. In the discussion whic] 
follows, a general limiting law is derived for the lower cut-of! 
point of the impaction efficiency. This is then applied to th: 
sphere, the blunt body for which the most theoretical and ex 
perimental information is available. Effects such as partic 
disintegration in the shock or melting at the surface are no! 
considered. It is interesting to note that impaction calcula 
tions at high Mach numbers will in many cases be relatively 
simple (once the gas velocity distribution has been estab 
lished), since the starting point for numerical calculations i- 
well-defined compared with subsonic flows. 


Lower Cut-off Point for the Impaction Efficiency 


The equation of motion for a particle traveling along thx 
stagnation streamline is given by 


1 
Cp prA |v — vp| (v — vf) [1| 
Immediately after the shock, the particle possesses its initial 
(supersonic) velocity while the gas velocity becomes subsonic. 
At low impaction efficiencies, however, the velocity of the 
particles near the forward stagnation point is small and the 
drag is given by Stokes law 


d 
m “4 = —f(v — vp) [2] 


The fluid velocity along the stagnation line can be expanded 
in a Taylor series around the stagnation point, r = a 


dog dup (r a)? 
(@ ) a+ (Sr) 


Substituting in Equation [2] with »v = dr/dt = dy/dt, there 
results 


d 
my + fo = f(bry + by? +...) [4] 


where the coefficients b, represent the derivatives and numeri- 
cal coefficients in Equation [3]. It has been shown by 
Poincaré (see, for example (2)) that the singularities of 
Equation [4] are the same as those of the first approximation 

dv fby — fv 
dy mv [5] 


The behavior of the v-y curve in the neighborhood of 0,0 can 
be determined (2) from Equation [5] by examining the roots 
of the characteristic equation 


S? + fS — mfb, = 0 [6] 
These are 


2 


When the discriminant f? + 4mfb, is negative, i.e., (mV/fa) = 
N,; > —(V/4ab:), the roots are complex conjugates, and the 
system has a focal or spiral point at 0,0 (Fig. 1a). The first 
intersection of the spiral with the v axis corresponds to finite 
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velocity at the stagnation point. The rest of the spiral has no 
physical significance, since the particle cannot pass the surface 
of the sphere. 

When the discriminant is positive, VN; < —(V/4ab,), the 
roots are real, and both are negative, since 6; is a negative 
number. The system has a nodal point at 0,0 (Fig. 1b), which 
corresponds to zero particle velocity at the forward stagnation 
point and zero impaction efficiency. If the discriminant 
vanishes, the roots are equal, and the system again has a 
node as its singularity. Thus 


= [8] 


represents a lower limit below which the impaction efficiency 
vnishes. 

This analysis is an extension of one originally suggested by 
‘Jayvlor (3) for subsonic flows. It neglects boundary layer 
e fects and is probably best applied when the particle diameter 
is larger than, or of the order of, the boundary layer thickness. 
Tor most real (viscous) flows, the first derivative of the 
\clocity component normal to the surface with respect to the 
rormal vanishes because of the no-slip boundary condition 
end the continuity relation. In this case, the equation of 
particle motion cannot be put into the form of Equation [5], 
and an analysis of the type developed above is not per- 
rissible, 


Application to Flow Over a Sphere at High Mach Numbers 


A number of investigators have chosen the sphere for 
theoretical or experimental studies of the fluid mechanics of 
high Mach number flow around blunt bodies. Lighthill (4) 
has developed a theory for the inviscid flow outside the laminar 
boundary layer assuming that near the forward stagnation 
point the shock front is spherical and concentric with the 
sphere, and that the flow between the shock and the surface is 
incompressible; this limits the analysis to the region within 
about 10 deg of the stagnation point. Although inviscid and 
incompressible, the flow is not irrotational because of the 
vorticity produced at the curved shock front. Using the 
strong shock approximation, Lighthill derived an expression 
for the vorticity at the shock surface. From the defining ex- 
pression for the vorticity, he was able to calculate the stream 
function for the flow between the shock and the surface 


= Wr) sin? [9] 


where 


30K 


[scx 1) — 5K(K — 4) @) + 2AK—1\(K 


The position of the surface relative to the shock (and the de- 
tachment distance) is obtained by solving for the greatest root 
of Y(r) = 0 which is less than c. 

The radial fluid velocity along the stagnation line is given 
by 


and 


n= dr |... (@=0) [11] 


b 


Fig. la Spiral point Fig. 1b Nodal point 


This expression can be combined with V(a) = 0 to give 


ne K(K »| [13] 


Combining Equation [14] with [8], the critical impaction 
parameter becomes 


K 
(K — 1)? [15] 


This is plotted in Fig. 2 as a function of K. The ratio of the 
detachment thickness to the radius of the sphere 6/a was taken 
from Lighthill’s Fig. 6. 


Sample Calculation 


As an example of the application of these results, one may 
calculate the diameter of the smallest particle collectible by 
a sphere of radius 1 ft traveling at a velocity of 8000 ft per sec 


at an altitude of 36,000 ft. Under these conditions, according 
to the tables of normal shock parameters prepared by Huber 
(5), temperature and pressure before the shock are 217 K 
and 0.224 atm, and after the shock, 2600 K and 18.5 atm. 
The upstream Mach number is 8.25 and K = 6.88. From Fig. 
1, the critical impaction parameter is 0.208. Assuming that 
the particles are spherical 


= “— = 0,208 [16] 
a 


The viscosity of the air can be calculated from Equation [10] 
and Fig. 12 of the review of the properties of high tem- 
perature air by Hansen and Heims (6). At 2600 K, the effect 


a ay") 
- | — 1) — 6(K — 1)(K — 6) ] [12] 
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of dissociation on viscosity appears to be negligible, and p is 
calculated to be 7.13 (10)~4 gm/em sec. For particles of unit 
density the minimum collectible diameter is 


[171 
8000 


= 5.78 microns 


It is doubtful that particles this large exist at 36,000 ft, and 
thus no deposition would occur on the sphere. For a sphere 
one hundred times smaller, i.e., about } in. in diameter, the 
minimum collectible particle would have a diameter of about 
0.6 micron. 


Nomenclature 

a = radius of sphere 

A = particle cross-sectional area 

b,» = Taylor series expansion coefficients 

c = radial position of shock 

Cp = drag coefficient 

d = particle diameter 

f = coefficient of resistance in linear drag law 

K = ratio of gas density downstream of shock to density up- 
stream of shock 

m = particle mass 

N; = impaction parameter 

n = radial position from center of sphere 

S = roots of the characteristic equation 

i = time 

v= particle velocity along stagnation line 

vy = gas velocity along stagnation line 

V = velocity of gas at infinity relative to sphere 

y =r-a 

6 = shock detachment distance 

@ = angle measured from forward stagnation point 

“ = gas viscosity 

p = particle density 

pr = gas density 

y = stream function 

WW = radial dependence of stream function 
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Simple Analytical Approximation to the 
Equation of State of Dissociated Air' 


W. L. BADE? 


Aveo Research and Advanced Development Division, 
Wilmington, Mass. 


NUMBER of investigators (1 to 4)* have calculated th« 

thermodynamic properties of equilibrium air to 10,000 K 
or higher. At temperatures exceeding about 2000 K, the 
density is found to be a rather complicated function of tem- 
perature and pressure. The compressibility factor Z(p, T) = 
pM)/pRT increases rapidly with temperature in regions where 
QO, or Ne molecules dissociate into atoms, and varies slowly in 
intermediate regions. The equation of state has been available 
only in tabular and graphical forms and in terms of fairly 
complicated approximate analytical representations. 

A simple analytical approximation to the equation of state 
of equilibrium air would be useful in many types of calcula- 
tions pertaining to boundary layers, aerodynamic heating, 
skin friction, inviscid flow problems, etc. Such a simple ap- 
proximation may be found by noting that, in regions where 
the density changes rapidly with temperature owing to occur- 
rence of dissociation, the specific enthalpy h also changes 
rapidly with temperature because of increase of enthalpy of 
dissociation. This crude argument makes it plausible that 
density should be a simpler function of enthalpy and pressure 
than of temperature and pressure. Such is indeed the case. 
Fig. | is a logarithmic plot of p/p vs. h/RT» for the four pres- 
sures 0.1, 1, 10 and 100 atm. The data are from Hansen‘ (4). 
It is seen that the points for each pressure are fitted rather 
well by a straight line, that the four lines corresponding to the 
different pressures can be drawn so as to pass through a single 
point at 


Il 


31.9RT> = 1080 Btu/Ib = 600 cal/gm 

= 0.124p. = 0.0100 lb/ft? = 1.60 X 10-4 gm/em3 

and that the slopes of the lines vary in a regular way with the 

logarithm of the pressure. A simple curve-fitting operation 

leads to the formulas 

p/p: = (p/pr)(h/hi) ~* [2] 


x = 0.70 + 0.04 logic (p/p:) [3] 


> 
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* Hansen’s tables provide the most convenient basis for seeking 
a correlation of this type, because they give the thermodynamic 
functions with pressure as an independent variable from 0.0001 
to 100 atm. 
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Fig.1 Density (p) of equilibrium air as a function of pressure (p) 
and enthalpy (h/ 


where p; and hy; are given by Equation [1], and p,; = 1 atm. 
Within the range of conditions 


31.9 < h/RT, < 480 O.atm< p< 100atm [4] 


Equation [2] agrees with Hansen’s results to within 5 per cent. 
In the same range, Equation [2] also’ agrees with the results 
of Logan and Treanor (3) to within 5 per cent, except in the 
region 0.1 < p < 1 atm, 60 < h/RT < 100, where the 
error rises to 8 per cent. In addition, Equation [2] is valid 
to within 5 per cent at pressures considerably greater than 100 
atm over the full enthalpy range indicated in Equation [4], 
and at pressures smaller than 0.1 atm over a reduced enthalpy 
range (120 < h/RT, < 400 down to 0.001 atm). 

Below the enthalpy h; = 31.9 RT», the ratio p/p is nearly 
independent of pressure and can be represented readily in 
various ways as a simple function of temperature or enthalpy. 
lt may be convenient, for some applications, to represent the 
density in this low temperature region by a function similar 
in form to Equation [2] and joining onto [2] at the point 
h=h,. Such a function is 


p/p: = (p/p.i)(h/h,) 4 (5) 


which agrees with the results of Hilsenrath et al. (5) to within 
7 per cent in the range 3.8 < h/RT) < 31.9,0.latm < p < 100 
atm. 

The dependence of p upon p through the logio(p/p,) term in 
the exponent of (h/hi) in Equation [2] is weak, so that, for 
processes not involving excessively large changes in pressure, 
the exponent can be replaced, to a close approximation, by a 
fixed value = x(p) = constant, where p is a mean or typical 
pressure. For example, by differentiating Equation [2], with 
xc = constant, and employing the equation 


DQ = dh — dp/p = 0 (6) 


defining an adiabatic process, one readily obtains an approxi- 


5 It was considered necessary to check Equation [2] against 
the accurate tabulation of Logan and Treanor as well as against 
Hansen’s tables, because the latter are based upon an approximate 
method which leads to errors (<5 per cent) of the same general 
size as those being sought in Equation [2]. 
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mate expression for the velocity of sound in dissociated air 


op Ss 1 hy h 


Within the range of enthalpy and pressure defined by Equation 
[41, Equation [7] agrees with Hansen’s results to within 7 per 
cent. 

The constant-exponent approximation to Equation [2] can 
also be used in the integration of Equation [6] to obtain a 
real air analog 


i) 
p exp [ - = constant [8} 
Pi hy 


of the ideal gas adiabatic equation of state p/TY‘7-) = 
constant. Equation [8] leads to an analog of the Bernoulli 
equation, for the compressible flow of dissociated air 


l1—zp 2h 


2 2h 6 


where 7, is the stagnation pressure and (p, p, h) is the state of 
the air at a point where the velocity is u. 

A further application of Equation [2] is in calculations of 
aerodynamic heat flux. For example, Fay and Riddell’s (6) 
expression for the laminar heat transfer rate at the stagnation 
point of a blunt body in dissociated air may be written, for the 
case of a strongly cooled surface, in the form 


Ps 


F(p,, hy) (du,/dx),'/*(hs {10a} 


where h,, represents the wall enthalpy, (du./dx), is the stag- 
nation point velocity gradient, and F(p,, h,) depends upon the 
Prandtl and Lewis numbers and the density, viscosity, stag- 
nation enthalpy and dissociation enthalpy of the air at the 
edge of the boundary layer. Equation [2], together with a 
crude fit to the temperature of equilibrium air as a function 
of enthalpy, provides an approximate analytical expression 
for this function F 


F(p., hy) ~ 


Ps [0.19 + 0.016 logio(ps/p:)] 
7.68 X 1074 {10b] 
1 


in lb/ft? see'/?. Equation [10] agrees with the Fay-Riddell 
formula to within 5 per cent in the range of conditions 100 < 
h,/RT, < 400,01 atm < ps; < 100 atm. 

Acknowledgment is made to Martin Zlotnick, who read this 
Note in manuscript and made helpful comments, and to 
Leonard Comeau, who carried out many of the numerical 
calculations. 


References 


1 Gilmore, F. R., ‘Equilibrium Composition and Thermodynamic 
Properties of Air to 24,000°K,”’ U. S. Air Force Project RAND, Research 
Memo 1543, Aug. 1955. 

2 Hilsenrath, J. and Beckett, C. W., ‘‘Tables of Thermodynamic Proper- 
ties of Argon-Free Air to 15,000°K.,’’ Arnold Engi ing Develop t 
Center, AEDC-TN-56-12, Sept. 1956. 

3 Logan, J. G., Jr. and Treanor, C. E., ‘‘Tables of Thermodynamic 
Properties of Air from 3000°K to 10,000°K at Intervals of 100°K,"’ Cornell 
Aeronautical Laboratory, Report no. BE-1007-A-3, Jan. 1957. 

4 Hansen, C. F., ‘‘Approximations for the Thermodynamic and Trans- 
port Properties of High-Temperature Air,’’, NACA TN 4150, March 1958. 

5 Hilsenrath, J. et al., ‘‘Tables of Thermal Properties of Gases,’’ Nat. 
Bur. Standards, Circular 564, Nov. 1955. 

6 Fay, J. A. and Riddell, F. R., ‘Theory of Stagnation Point Heat 
Transfer in Dissociated Air,’’ J. Aeron. Sci., vol. 25. 1958, pp. 73-85. 


299 


| 

1 

[9} 

| 

| 
| 

the 

10n 

[2] ; 
[3] 

stic 
per. 
ing 
mic 


Minor Circle Flight for Boost Glide 
Vehicles 


W. H. T. LOH! 


Chance Vought Aircraft, Dallas, Texas 


YPERSONIC gliders designed to strike targets located 
halfway around the Earth will fly great circle paths in 
reaching these targets and returning to the launch area. Great 
circle paths are also used in striking targets at shorter ranges. 
However, the possibility of reaching these shorter range tar- 
3 gets by flying paths other than great circles (minor circle 
paths) may be of interest. Recognizing that the great circle 
path is the shortest path between any two given points and 
recognizing that the great circle path is also the glider’s 
natural automatic gliding path (in contrast to the minor circle 
path which needs continuous accurate aerodynamic control to 
force the glider on the desired path) nevertheless, one might 
use a minor circle path simply for strategic reasons. The 
approximate range equation of such a glider flying along a 
great circle path, derived recently by Eggers and Allen? is 


2 D | 1 — (V,?/gRo) 
or 


2/\D L1 — (Vi?/gRo) 


On deriving the equation, they used the following assumptions: 


1 Flight is in planes containing the great circle are be- 
tween takeoff and landing. 
2 Small inclination angle to the horizontal, i.e., cos 6= 1, 
and sin 
3 Constant gravity acceleration, i.e., R, = Ro. 
4 Constant L/D flight. 
5 Rotation of Earth is neglected. 
6 Since (L/D) 6< 1, and V? d6/ds < 1/2(L/D)(dV?/ds), 
both small terms are neglected in the basic equations. 
7 Contribution to range of powered phase of flight is 
neglected. 
> Similar assumptions for minor circle flight may be written 
here: 


1 Flight is in planes containing the minor circle are be- 


L 1L 
g sin (1 + (3% sin a + V 


with Earth radius is negligible. Therefore R, = Ro; r=r, 
and A = Apo (see Figs. 1 and 2). 

4 Constant L/D flight. 

5 Rotation of Earth is neglected. 

6 Since (L/D) @sin a <1, and V? d@/ds < 1/2(L/D) > 
(dV2/ds) sin a, both small terms are neglected in the basic 
equation (provided @ is not too small). 

7 Contribution to range of powered phase of flight is 
neglected. 

With these assumptions and referring to Figs. 1 and », 
one obtains immediately the following three equations of 
motion normal and parallel to the direction of flight 


L cos a = mg cos X [le] 


L sin a — mg sin A cos 0 = —m — 
Te 


dv 1 
—D + mg sin Xd sin 6 = 


d¥ cos# 1 
= - [1d 
ds r To 


dV? 
ds 


[1c 


ds 


Combining Equations [lb, le and 1d] and using the ap- 
proximations as listed, one obtains after simplification 


{2a 


dV? 


~ Resin 


Since (L/D)6@ sin a < 1, and V%(d0@/ds) < 1/2(L/D) x 
(dV2/ds) sin a, they may be dropped out. Equation [2a) 
therefore simplifies to 

29 sin Xo 
(L/D) sina 


av? 2 
| ] [2b 
ds Ro (L/D) sin Xo sin a 


Combining Equations [la and 1b] and using the approxima- 
tions as listed, one obtains after simplification 


v2 
gRo sin? Xo 


sin a => sin Xo 
gRo 


Substituting Equation [3] into [2b] and integrating it between 
initial and final conditions, one obtains 


1 Vr? 2 Vr? 
) sin Xo sin (22 ) 


1 — ———_ 
Ry L gRo sin? Xo 
R= (5) sin do In ( V2 

1 


sin? Xo 


tween takeoff and landing. 
E 2 Small inclination angle to the minor circle plane hori- 
Be zontal, i.e., cos @ = 1, and sin 6 = @. 


3 Constant gravity acceleration, i.e., Rj = Ro. Note that 
this implies here that the altitude of the glider when compared 


[4a} 


ve \2 (2) 
-2(—)+1 
) sin Ao sin 
As maximum range is approached L/W = 1, Vr?/gRo sin? do 
then is small compared to one, and it may be neglected in Equa- 
tion [4a] without much loss of accuracy. In’ this case, the 


maximum range for the glider flying along a minor circle 
path becomes simply 


r 
(2) sin don V;? 


gRo sin? x 1 


‘ Received Jan. 27, 1959. 
1 Staff Engineer to Assistant Chief Engineer, Technical. 
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[4b] 


| V;2 V;?2 2 
sin Ao (: V2 ) 
sin? Ao 


2 Eggers, A. J. and Allen, H. J., ‘‘A Comparative Analysis of the 
Performance of Long-Range Hypervelocity Vehicles,’?” NACA 
TN 4046, Oct. 1957. 
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Fig. 1 Three-dimensional view on minor circle flight 


—Minor Circle 
Earth Surface 


LSinc 
Fig. 2 Two-dimensional view on minor circle flight 


When Xo = (2/2), the minor circle becomes a great circle, and 
the present range of Equations [4a and 4b] are reduced to the 
range equation derived for great circle paths by Eggers and 
Allen. Consequently, the great circle range equation may be 
considered as a special case of the range equation derived here. 
The reduction of range on minor circle path is mainly due to 
the fact that part of the dynamical energy possessed by the 
glider is now used to control and force the glider in and within 
the minor circle path being particularly specified; conse- 
quently, this dynamical energy is simply dissipated without 
accomplishing any useful range. 


Nomenclature 

D drag 

g acceleration due to force of gravity 
L lift 

m mass 


distance shown in Figs. 1 and 2 
distance shown in Figs. 1 and 2 


re = radius of curvature of flight path, shown in Fig. 2 
R, = radius of Earth 

R, = glider distance from center of Earth 

R = flight range, measured along surface of the Earth 
s = distance along flight path 

= time 

V = velocity 

W = weight 

6 = angle of flight path to horizontal 

¢ = partial range angle, rad 

® = total range angle, rad 

WwW = remaining range angle, rad 

Xo = angle defined in Fig. 1 

\ = angle defined in Fig. 1 

a = angle defined in Fig. 1 

Subscripts 

i = initial conditions 

F = final conditions 
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Aerodynamic Moment on Bodies 
Moving at High Speed in the 
Upper Atmosphere 


W. A. GUSTAFSON! 


Lockheed Aircraft Corp., Sunnyvale, Calif. 


HE ATTITUDE control of vehicles traveling in the upper 

atmosphere is of considerable importance and requires a 
knowledge of the aerodynamic characteristics of various body 
shapes. At sufficiently high altitudes and correspond- 
ingly low densities, the aerodynamic and _ gravitational 
moments acting on a body may be of similar magnitude, 
thereby making an analytical] expression for the aerodynamic 
moment useful for an investigation of the dynamics of the 
vehicle about its center of gravity. 

In the following discussion, only altitudes above 100 statute 
miles are considered, and consequently free molecular flow 
theory is applicable. It is assumed that the surface tem- 
perature is primarily controlled by thermal radiation phe- 
nomena, and thus 7’, can be determined by a separate analysis. 

The three body shapes of greatest practical interest are the 
cylinder, the cone or cone frustum, and the spherical segment. 
Applying free molecular theory to these shapes for the pur- 
pose of determining the aerodynamic moment results in in- 
tegrals which apparently must be evaluated numerically ex- 
cept for the case of the cylinder. However, analytical ex- 
pressions for the moment on a cone and on a spherical seg- 
ment can be obtained by neglecting the free stream thermal 
motion. This is a good approximation as long as the molecu- 
lar speed ratio is not too low. In the following paragraphs, the 
static aerodynamic pitching moment is determined about an 
arbitrary point on the body axis and for any angle of attack. 
It is assumed that the free stream strikes the body surface 
and gives up its total momentum, and that diffuse reflection 
then takes place with the atmospheric particles leaving the 
surface possessing a Maxwellian distribution corresponding 
to the surface temperature. It is also assumed that the mass 
of particles striking a unit of surface area per unit of time is 
equal to the mass leaving that area in a unit of time. Based on 
these assumptions, the forces acting on an element of area can 
be described and the moment arm determined. Thus, inte- 
grating over the body surface exposed to the free stream pro- 
vides the desired result. Since the method is straightforward, 
only the final results are given below. A more detailed de- 
scription of the derivation of the results is presented in (1).? 


Received Dec. 22, 1958. 

1 Associate Research Scientist, Aerodynamics Department, 
Missile Systems Division. 

2 Numbers in parentheses indicate References at end of paper. 
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Cone 


The moment M, about the point C is given below, cor- 
responding to the notation in Fig. 1. However, when the 
angle of attack is greater than the semi-vertex angle, part of 
the cone is shielded from the free stream. Hence, this area 
is not included in the integration. The final result is given 
in two parts. 


(a) B <6 

M. 
Gy Fa Bsn B42 tan B cos B — | + (1) 
(b) 6 

M. tan 


(1/2) pU? = 1?) LA sinB — - (ue! yi3)(A sin B + Bcos@ tan 5) + WV xT,/T; 


where 


O for a cone frustum 


A = cos B sin 6 [x — cos~ (cot B tan 5)] + sin B cos 6 sin [cos~! (cot 8 tan 6) | 
L 
B = cos B sin 6 sin [cos~ (cot B tan 4) | +5 B cos 6 — (cot B tan 6) — sin [2 (cot B tan 
Sphere 


On a complete sphere, the total force is only a drag forces 
so that the moment of a sphere about any point is easily de- 
termined. However, in some cases it may be required to cap a 
cylinder or a cone frustum by a hemisphere or a spherical seg- 
ment, respectively. According to Fig. 2, 6, is the angle which 
specifies the size of the spherical segment. If a spherical cap 
for a cone frustum is required, then @, is determined by the 
cone geometry, since 6 = 7/2 — @,. In a manner similar to 
that for the cone, the analysis is done in two parts, since for 
angles of attack greater than 7/2 — 6, part of the surface is Cc 
shielded, and the limit of the ¢ integration then depends both 
on 6, and 


Fig. 2. Geometry of spherical segment 


(a)6B 6 


M. 1 
(1/2) Tr? sin B cos B sin? 6, (l + r cos 6;) g E 3 cos 8; (sin? » [3] 


to 


(b) 8B > 1/2 6 


sin B cos B - + 2r sin? 6; cos 0, + 5 008 B + sin? + 3 sin B cos? BQ — sin B cos B VQ + 


l 1 
sin? B + ; sin B cos B cos~! (cot B cot E sin? 0, + 2r cos? 6, — 3r cos 6; (1 cos? 


l 
9 sin? B cos 0, VP — 4 sin B cos? B sin! [1 + 2 cot? B — 2 esc? B sin? 6,| + : sin? 6(2 sin? B — cos? B) ~/P3 + [4] 


l 
3 sin? tan! ~/cos? 6,/(sin? 6, — B) — 5 sin B cos B sin-! [—2 cos? B ese? 6, + uf 


2lr? 2 1 = 1 
Wer. /T; sin B [1-5 — = cos 8, (sin? + g sin? B cos B + 3 sin B 


sin B cos? B VQ + sin B cos (: ; cos? a) cos~' (cot B cot a 
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+ 
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In Equation [4] 
P = sin? 0, ese? B — cot? B 


1 — cot? B cot? 6; 
cot? B ese? 6; 


Q= 


Cylinder—Free Molecular Theory 


The analysis of the cylinder based on free molecular flow is 
given below using the notation and local coordinate system of 
Talbot (2), except that his T-axis is reversed in direction. The 
components of the pressure and shear forces acting on an ele- 
ment of area in the N and T directions for a= ¢ = o’ = 1 
are 


(1/2) pU? sin B sin @ + VE) sin @ e + sin? + sin? B sin + sin B x 


TN 


(1/2) pU? 


(1/2) pU? 


The principal values of the are sine terms between — 7/2 and 
x ‘2 are to be taken in the preceding expression. The negative 
root of the radical in the are tangent terms should be used, 
and the corresponding angle between —7/2 and 0 should be 
selected. The positive root of all other radicals is retained. 
}-quation [4] is valid only to 8 = 7/2, although an equation 
for B > 2/2 can be derived in a similar manner. 


{1 + erf (Esin d)] [5] 


— sin B cos? sin? o+ Esin [1 + erf (Esing)]} [6] 


1 = 
fe + War Esin [1 + erf(E sin (71 


where 
= SsinB 
B = the angle of attack of the cylinder 


Hence, the pitching moment about a point on the cylinder 
axis a distance L from the origin of the y coordinate (see Fig. 3) 
is given by 


2 
[8] 
Evaluating the integrals gives 


Fig. 3. Geometry of cylinder 
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tr2y. sin B cos B [9] 


where J(£?/2) and J,(&#/2) are modified Bessel functions of the 
first kind and of zero and first order, respectively. The ends 
of the cylinder have been neglected in determining Equation 


(9). 


Nomenclature 
p = atmospheric density 
U = flight velocity 


7; = ambient temperature 


T, = surface temperature 

R = gas constant _ 
S = molecular speed ratio = U/~V2RT; 
References 


1 Gustafson, W. A., ‘‘The Newtonian-Diffuse Method for Computing 
Aerodynamic Forces,’ Technical Memorandum LMSD-5132, Lockheed Air- 
craft Corp., Aug. 28, 1958. 

2 Talbot, L., ‘‘Free Molecular Flow Forces and Heat Transfer for an 
Infinite Circular Cylinder at Angle of Attack,’’ 7. Aeron. Sci., vol. 24, 1957, 
p. 458. 
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Book Reviews 


Turbojet Fundamentals, by Howard E. 
Morgan, McGraw-Hill Book Co., Inc., 
New York, 1958, 104 pp. $4.00. 

Reviewed by P. Roy CHoupHURY 
University of Southern California 


The very elementary nature of this book 
makes it suitable for those aircraft techni- 
cians and mechanics who have little pre- 
liminary fundamental technical training. 

Having had occasion to lecture on after- 
burners and flame stabilization to groups 
of, high school students who had no pre- 
vious training in thermodynamics, fluid 
mechanics,/etc, the reviewer fully appre- 
ciates the problems faced by the author in 
writing the book. Asa matter of fact, the 
author has done a commendable job in 
explaining the operating principles of the 
various components of a jet engine. The 
use of numerous appropriate sketches and 
figures makes the subject matter easily 
readable. In some instances, however, 
the author fails to stimulate the interest of 
a somewhat brighter student by not adding 
a note of explanation to various topics. 
The omission of the variable stator type of 
compressor is rather unfortunate. The 
author defines isentropic as “inactive or 
static condition,’ and because of the un- 
orthodox nature of this definition an ex- 
planation is certainly justified. 


Jet Engine Manual, by E. Mangham and 
A. Peace, Philosophical Library, Inc., 
New York, 1955, 133 pp. 

Reviewed by P. Roy CHoupHURY 
University of Southern California 


Unlike “Turbojet Fundamentals,” “Jet 
Engine Manual”’ gives a detailed descrip- 
tion of the components of a jet engine. 
The authors have discussed the operation 
of the fuel, lubrication, starter and control 
systems of turbojet and turboprop engines. 
In addition to these, the construction and 
functional details of the two types of 
engines are discussed. This manual will 
be very useful to those who are concerned 
with the maintenance and operation of 
turbojet and turboprop engines. 

The lack of adequate sketches and fig- 
ures, in many sections, makes it rather 
difficult to visualize the working of various 
small components. The addition of ap- 
propriate figures will greatly improve the 
utility of this manual. 


Pumps, Fans, and Compressors, by A. de 
Kov4ts and G. Desmur, Blackie and 
Son, Ltd., Glasgow, 1958, xii + 327 pp. 
45s. 

Reviewed by R. C. BinpER 
Purdue University 


This book is an English translation (by 
R. 8. Eaton of Manchester College) of the 
original French work ‘Pompes, Ventila- 
teurs, Compresseurs’’ (published by 
Dunod, Paris, 1953). This book discusses 
the design, construction and operation of 
centrifugal and axial flow pumps, fans and 
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Ali Bulent Cambel, Northwestern University, Associate Editor 


compressors. It has a practical industrial 
viewpoint and assumes some knowledge of 
basic fluid mechanics. In making design 
studies, some engineers rely heavily on 
purely theoretical relations. As a con- 
trast, this book by Kovdts and Desmur 
does not follow this purely analytical ap- 
proach; rather it relies directly on experi- 
mental data and empirical relations in 
many cases. 

The book is divided into seven parts 
with a total of 34 chapters. Some of the 
chapters are short, but the authors man- 
age to bring out pertinent main topics in a 
book of moderate size. Part 1, on funda- 
mentals, treats incompressible flow in 
channels, the flow around blades and 
through cascades, and fluid power losses. 
Part 2, on the design and construction of 
impellers and diffusers, covers performance 
characteristics, velocity triangles, specific 
speed, energy losses, cavitation, operation 
of a pump as a turbine, problems affecting 
the choice of a machine, impeller design, 
diffuser, volutes, and axial thrust and 
balancing devices. Part 3, construction 
of pumps and fans, describes mechanical 
details, such as casings, impellers, shafts, 
bearings and stuffing boxes. Part 4 de- 
scribes pumps for such special duties as 
boiler feed, condensate, sewage, dredging 
and chemical applications. Part 5 gives 
some examples of the calculation of the 
principal dimensions of various types of 
pumps and fans. The various numerical 
examples and illustrations are very helpful 
to one interested in learning something 
about design. Part 6 covers briefly the 
operation of pumps and fans, along with 
the selection of type. Parts 1 through 6 
are primarily for incompressible flow. 
Part 7, on turbo-compressors, includes a 
treatment of basic thermodynamics, the 
calculation of losses, the selection of com- 
pressor type, compressor characteristics, 
the mechanical details of compressors, and 
an example of the thermodynamic calcula- 
tion of a compressor. E 

This book is more of the practical ref- 
erence type of book rather than a college 
text. It should be useful to the engineer 
who may be concerned with practical 
problems involving the design, construc- 
tion and operation of these machines. 


Fundamentals of Gas Dynamics, edited 
by Howard W. Emmons, Vol. III of 
“High Speed Aerodynamics and Jet 
Propulsion,” Princeton University 
Press, Princeton, 749 pp. $20.00 
Reviewed by A. Gross 
Fairchild Engine Division 


This book, the sixth to be published in 
what will be a twelve volume series, is a 
valuable addition to the aerodynamics 
literature. Gas dynamics is treated here 
as a branch of physics and applied mathe- 
matics, and this volume reflects the recent 
rapid growth and increasing breadth of 
this subject. It contains eight major 
sections each written by persons well- 
known for their original contributions to 


the subject. Because each section of this 
book strongly reflects its author, each is 
treated separately in this review. 

A Equations of Gas Dynamics; H. §. 
Tsien, pp. 3-63. This section presents « 
classical continuum approach to the basi: 
equations of gas dynamics. Equations 
are developed for quite general flows, and 
kinetics are discussed only as they effec 
relaxation times. The equations of mo- 
tion are expressed finally in general orthog- 
onal coordinates. It is unfortunate thi: 
work was completed before the current in 
terest in magnetohydrodynamics, sinc 
this subject may well be considered as : 
domain in gas dynamics, 

B One-Dimensional Treatment of 
Steady Gas Dynamics; L. Crocco, pp 
64-349. This section is nearly a book ir 
itself. It treats many interesting and im- 
portant subjects from an engineering point 
of view. Although correction coefficients 
are frequently employed, the difficult sub- 
jects which are discussed require occa- 
sional empirical treatment. Some of the 
topics covered are: Thermodynamics of 
real substances, shock waves and “‘pseudo- 
shocks” (boundary layer-shock interaction 
in one dimension), nozzle flow, supersonic 
diffusers, friction in long ducts, heat ex- 
change and combustion, flows with piece- 
wise uniformity, and mixing. The analy- 
sis and physical descriptions are coupled 
with frequent graphs and 67 references. 

C One-Dimensional Treatment of 
Nonsteady Gas Dynamics; A. Kantro- 
witz, pp. 350-415. This section reflects 
activity in this field up to about 1952. 
The fundamental equations are again de- 
rived. Plane isentropic waves of large 
amplitude and weak and strong shocks are 
discussed. The formation of shocks in 
channels and nozzles are analyzed along 
with their stability. Intense explosions 
and implosions are briefly treated, and, 
finally, a very short section on the applica- 
tion of pressure waves to heat engines is 
presented. 

D The Basic Theory of Gasdynamic 
Discontinuities; W. D. Hayes, pp. 416- 
481. This section presents a theoretical 
introduction to gasdynamic discontinui- 
ties from a quite general point of view, and 
appears to contain much original work. 
The basic conservation laws are presented. 
The existence of various solutions is ex- 
plored using thermodynamic and stability 
arguments. Normal shocks, exothermic 
discontinuities and Chapman-Jouguet con- 
ditions are examined. The internal struc- 
ture of a wave is examined for stability 
with regard to its ability to jump to an- 
other condition. The Navier-Stokes 
structure of shocks and exothermic dis- 
continuities are treated, and, finally, some 
physics of shocks, such as radiation, ioni- 
zation and_ relaxation are briefly 
mentioned. Fifty-one references are cited. 

E Shock Wave Interaction; H. Pola- 
chek and R. Seeger, pp. 482-525. Shock 
waves are again treated as discontinuities 
with emphasis on normal and oblique re- 


(Continued on page 308.) 
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MISSILE 


This huge antenna is part of an electronic and 
optical system that RCA installed and is 
operating on the S.S. American Mariner. The 
purpose of the equipment, for which the ship 
has been refitted, is to provide the most precise 
data yet obtained at sea on missile flights over 
a range extending from Cape Canaveral, Fla., 
to the area of Ascension Island. The project is 


Tmk(s) ® 


sponsored jointly by the Advanced Research 
Projects Agency, Department of Defense and 
the Army Ordnance Command. A scientific 
staff—most of them RCA personnel—will 
operate the equipment and report on missile 
performance from descent from space to final 
plunge, the data to be shared by all branches 
of the armed services. 
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STORABILITY 


Just one of the many outstanding attributes of 


DIMAZINE 


unsym-Dimethylhydrazine 


DIMAZINE has many characteristics that 
make it excel as a storable fuel. 


’ Its remarkable stability has been proven in 
tests of more than two years’ duration, even 

= in mild-steel vessels. Its high thermal sta- 
i bility and low freezing point permit wide 
variations of storage temperatures. It is not 
shock sensitive and is unaffected by many 

types of accidental contamination. It is high- 


Putting 


FOOD MACHINERY 
AND CHEMICAL 
CORPORATION 
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ly compatible with most metals and appro- 
priate sealant materials. 

Excellent storability, high performance, ease 
of handling, ample availability . . . these 
and many other advantages are combined in 
DIMAZINE, the outstanding storable fuel. 
We will be pleased to assist in your evalu- 
ation of DIMAZINE and to supply detailed 
data on its properties and handling. 


Work 


FOOD MACHINERY AND CHEMICAL CORPORATION 
Westvaco Chlor-Alkali Division 


General Sales Offices: 
161 E. 42nd STREET, NEW YORK 17 
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revolutionary in design 
and engineering as the technological transition 
from flying machine to space vehicle... Achieveda 
quarter of acentury ago, the only regulators which 
could be used for the first missile experiments 
in the Navy Experimental Station at Annapolis 
... Still today the standard for the critical control 
of high pressure fluids... Now on the threshold 
of space, Grove is uniquely prepared—through 
advance research and development—to meet 
the challenge of even higher pressures | | 


GROVE POWREACTOR DOME REGULATOR MODEL GH-408—50-6000 PSI INLET...5-3000 PSI REDUCED PRESSURE 


GROVE VALVE and REGULATOR COMPANY 


6529 Hollis St., Oakland 8, California » 2559 W. Olympic Blvd., Los Angeles 6, California 


Offices in other principal cities Subsidiary of Walworth 
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ENGINEERS - SCIENTISTS 
<> 
Are you getting 
these advantages 
in your 
present position? 


Advanced studies at Boeing include 
boost-glide vehicles, anti-ballistic 
missile area defense systems, super- 
sonic missiles, interplanetary sys- 
tems, and manned and unmanned 
space vehicles. 


Space-Age 
Assignments 


Advanced Research 
Facilities 


Boeing’s research, development 
and engineering facilities are 
among the most varied, extensive 
and complete in the industry. Your 
progress is faster. And you work 
with men who have “written the 
book” on the state-of-the-art. 


Opportunity for 
Advancement 


At Boeing, you work with small, 
integrated teams where initiative 
and ability get plenty of visibility. 
A personal merit review every six 
months assures continuing oppor- 
tunity for individual recognition. 
And you are eligible for advance- 
ment to a higher level at any time. 


C] Better Family Living 


There’s a real bonus of better liv- 
ing in the heart of the evergreen 
Puget Sound country. Four full 
seasons of comfortable moderation 
in an area where scenery and out- 
door recreation are nation-famed. 
To give you an idea: per capita 
ownership of sail and power boats 
in Seattle is the highest in America. 


| Write today for your free 
| copy of the 24-page book 
“Environment for Dy- 
namic Career Growth.” 
Indicate your degree(s) 
and field of interest. Ad- 
dress: Mr. Stanley M. 
Little, Boeing Airplane 
Co., P.O. Box 3822 - JPA, 
Seattle 24, Washington. 


flections. Some excellent photographs of 
such waves are presented. Some mention 
of detonation, refraction of shocks and 
particle models of continuum are made. 
Sixty-one references are cited. 

F Condensation Phenomena in High 
Speed Flow; H. G. Stever, pp. 526-573. 
Condensing vapors are first treated as an 
equilibrium problem. The kinetics of 
condensation are theoretically treated and 
then applied to high speed flow. Some 
experimental techniques are outlined; 
their results presented and compared with 
theory. It is concluded that the phe- 
nomena are well understood, but thermo- 
dynamic data and more precise and usable 
calculations are needed. 

G Gas Dynamics of Combustion and 
Detonation: 

1 Aerothermodynamic Problems of 
Combustion; T. von pp. 574- 
584. The case of one-dimensional flow of 
a gas mixture in which chemical reactions 
take place is formulated. Approximations 
are made and the eigenvalue flame speed 
problem is discussed. 

2 Flow Discontinuities Associated 
with Combustion; H. W. Emmons, pp. 
584-621. The theoretical formulation of 
the flame speed problem is presented along 
with a discussion of its inadequacies. 
The flame front is then treated as a dis- 
continuity in one- and two-dimensional 
flows. Flame generated vorticity is dis- 
cussed. Finally, the stability of flame 
discontinuities is treated analytically. 

3 Gas Dynamical Aspects of Deto- 
nation; G. I. Taylorand R.S. Tankin, pp. 
622-686. The Chapman-Jouguet theory 
of detonation is presented along with sev- 
eral numerical examples. The possibility 
of strong and weak detonations is discus- 
sed in some detail, and the transformation 
from deflagration to detonation is ex- 
amined for some particular cases. Inter- 
esting sections on spherical detonation, 
spinning detonation and explosions in 
solids and liquids are presented. 

H_ Flow of Rarified Gases; S. Schaaf 
and P. Chambré, pp. 687-739. This sec- 
tion gives a modern summary of this sub- 
ject including free molecular flow, heat 
transfer, aerodynamic forces, nonsteady 
flows, slip flow and experimental results. 

This book despite its size, does not at- 
tempt to treat all of gas dynamics. Tran- 
sonics, mixed flow problems and boundary 
layers are not included in this volume. It 
suffers somewhat from repetition, resulting 
from the numerous contributors. The fine 
photographs are often located far from 
their associated text. These are minor 
distractions however, and this undoubt- 
edly is an outstanding work in gasdynamic 
fundamentals and will be a standard ref- 
erence for future research efforts in this 
fascinating and still growing field. 


Principles and Applications of Random 
Noise Theory, by Julius 8S. Bendat, 
John Wiley and Sons, Inc., New York, 
1958, ix+ 422 pp. $11.00. 

Reviewed by G. J. Murpuy 
Northwestern University 


According to the author, this book was 
written for advanced senior students and 
graduate students in engineering, applied 
mathematics, or physics. The author’s 


purpose is to “present a systematic de- 
velopment of topics deemed to be most 
fundamental in random noise theory, and 
show how to apply the theory in certain 
areas of far-reaching importance which are 
not covered elsewhere.” 

His primary objective in the first four 
chapters is “to explain principles of ran- 
dom noise analysis and optimum filtering 
techniques in understandable engineering 
language, emphasizing physical meanings 
and mathematical expressions.’”’ In the 
opinion of this reviewer he has partially 
succeeded, in that he does discuss, in 
“engineering language,” the principles cf 
random noise analysis and optimum filte’- 
ing techniques, and he does emphasize 
physical significance. 

However, his attention to mathematic: | 
detail is far from adequate. For exampl:, 
on page 21, in making the transition fro1 
(1-60) to (1-69) he neglects to mentioy 
that the conditions under which one ma 
change the order of differentiation and in- 
tegration are restricted. As a result of 
hig failure to recognize this restriction, h: 
concludes that for all stationary randon: 
processes the slope of the autocorrelatio:: 
function is zero for zero argument. Ther 
on the next page he presents an example ii: 
which this conclusion is contradicted, and 
on pages 204 through 206 he attempts a» 
great length to force an (unjustifiable 
agreement of that example with the incor- 
rect conclusion. 

In the remaining six chapters of the 
book, the author’s objectives are ‘“‘to show 
how: (a) to formulate certain difficult noise 
problems; (b) to obtain their mathematica! 
solutions, and (c) to draw proper physica! 
designs and interpretations.’ Particu- 
larly noteworthy inclusions are the work of 
S. O. Rice on the distribution of a fixed 
sine wave plus random noise, a chapter on 
envelope detection, a chapter on analogue 
computer techniques, and a liberal discus- 
sion of time-varying filters. 

The book is enhanced by the use of 
chapter summaries, immediately following 
the table of contents, giving the reader a 
concise statement of the scope of each 
chapter, and by the effective use of modern 
techniques of layout, printing and binding. 

This book is recommended to both the 
university student and the research worker 
interested in random noise theory. 


Nuclear Rocket Propulsion, by R. W. 
Bussard and R. D. De Lauer, McGraw- 
Hill Book Co., New York, 1958, 370 pp. 


Reviewed by Donatp H. LouGuripGE 
General Motors Corp. 


The title of this book will doubtless at- 
tract much attention, since the subject 
matter has at present so much technical as 
well as political glamor. Potential read- 
ers, however, should be reminded that due 
to existing military secrecy, a complete re- 
view of present U. 8S. knowledge of the sub- 
ject cannot be expected. Accordingly, 
only overall general discussions of the state 
of the art are presented, with no direct dis- 
cussion of the KIWI or SNAP projects. 
Chapters on basic rocket performance and 
system analysis contain much engineering 
data and formulas which have been col- 
lected and are useful for the design en- 


(Continued on page 311.) 
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ELECTRONICS 


Inertial Guidance & Navigation + Digital 
Computer Development «+ Systems Engi- 
neering « Information Theory « Telemetry- 
SSB Technique « Doppler Radar « Counter- 
measures « Radome & Antenna Design 
* Microwave Circuitry & Components «+ 
Receiver & Transmitter Design + Airborne 
Navigational Systems * Jamming & Anti- 
Jamming + Miniaturization-Transistoriza- 
tion Ranging Systems Propagation 
Studies * Ground Support Equipment 

A new $14,000,000 Research Cen- 
ter—to be completed this year —is 

part of Republic’s far-ranging R&D 
programs aimed at major state-of- 
the-art breakthroughs in every flight 
regime and environment. 
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ENGINEERS-SCIENTISTS 


Engineers and Scientists at Republic Aviation—with a 
turn for dry professional humor — say that they’re asked 
to solve every problem of upper atmosphere and 

space flight with the factor 


V/ > 


This human application of the relativistic theory 
of space-time relationships is predicated on the conviction 
that the creative content of technological thinking 
can be immeasurably expanded and enriched — 
given a propitious environment. 


That the environment at Republic is propitious, 
is evident from the results. 


Engineers and scientists like it. They thrive on it. 
Technical ideas of a “revolutionary” character — 
rather than “evolutionary” — are appearing at a rate 
that exceeds the norm of even 5 years ago. 


In every profczzional area — research, development, 
experimental engineering —the goal is the same: 


PING TIME 


in terms of technological progress 


+ in exotic propulsion systems for 
space operation — plasma propulsion...advanced 
nuclear power applications... 


* in integrated electronic systems for 
flight vehicles to operate at every altitude 


> in supersonic and hypersonic weapons systems, 
both manned and unmanned 


For an invigorating environment, where 
new ideas flourish, look at the range of 
opportunities at Republic 


eee 


THERMO, PLASMA PROPULSION 


AERODYNAMICS 

Theoretical Gasdynamics + Hyper- 
Velocity Studies « Astronautics Pre- 
cision Trajectories » Airplane/ Missile 
Performance « Air Load and Aeroelas- 
ticity « Stabitity and Controls + Flutter 
& Vibration + Vehicle Dynamics & 
System Designs «+ High Altitude 
Atmosphere Physics « Re-entry Heat 
Transfer Hydromagnetics Ground 
Support Equipment 


Plasma Physics + Gaseous Electronics 
* Hypersonics and Shock Phenomena + 
Hydromagnetics + Physical Chemistry 
« Combustion and Detonation « Instrumen- 
tation + High Power Pulse Electronics 


NUCLEAR PROPULSION 

& RADIATION PHENOMENA 
Nuclear Weapons Effects + Radiation 
Environment in Space « Nuclear Power & 
Propulsion Applications + Nuclear Radia- 
tion Laboratories 


Send resume in complete confidence to: 
Mr. George R. Hickman, Engineering Employment Mgr, Dept. 10D 


Farmingdale , Long Island, New York 
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STANT 


. . - Requires no refrigeration . .. Easier to store than liquid oxygen NO, Oxidant for liquid 


Performance of Nitrogen Tetroxide as an oxidizer for many fuels is comparable 
to that of liquid oxygen. In convenience it is far superior to liquid oxygen. Molecular weight 
N204 can be shipped easily and stored indefinitely without refrigeration in Bolling Point 
ordinary carbon steel containers. It is a dense, mobile liquid that is noncorrosive Freezing Point 
if kept dry. The quantity of oxygen contained per unit volume of N2Ox is 1.01 Vaporization 
Kg/liter at 20°C. Critical Temp. 
is available in tank car quantities from Allied’s Hopewell, Virginia plant. 
For experimental purposes, N2OQ. is 
available in 125 lb. cylinders and 2000 Ib. _ Density of Liquid 
containers. Write for technical data, and &"sit of Gas 


Basic To information on prices and delivery. Vapor Pressure 
America’s 
Progress”’ NITROGEN DIVISION 


Dept. NT 4A-12-5, 40 Rector St., New York 6, N. Y. 
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rocket propellants 


92.02 
21°C 
—11.3°C 


99 cal/gm @ 21°C 
158°C 
99 atm 
0.36 cal/gm 

—10 to 20°C 
1.45 gm/ml at 20°C 
3.3 gm/liter 

21°C, at 1 atm 
2 atm at 35°C 
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gineer working in this field. References to 
original papers are given at the end of each 
chapter. However, for the technical 
reader who may be interested in more basic 
concepts and their interrelationships, such 
a purely engineering approach is never 
very satisfying. A chapter on heat trans- 
fer and fluid flow and the following one on 
materials constitute about a third of the 
book and proved, to the present reviewer, 
to be the most useful and convincing part 
oi the text. The following chapter on 
nucleonies depends upon a previous know]l- 
edge of elementary reactor theory, and 
without such, a reader would find little 
usefulness in this cursory survey. The 
treatments of system control and testing 
could just as well apply to any power re- 


Technical Literature Digest 


actor as to the particular application to nu- 
clear rockets. The last chapter on ad- 
vanced and exotic systems gives a short 
discussion of gaseous-fission, thermome- 
chanical and accelerator systems which are 
under such wide study at present. 

It is to be regretted that with such a 
great deal of work as has evidently been 
expended by the authors in collecting the 
wide range of engineering data here pre- 
sented, more direct reference to the work of 
the Los Alamos and Livermore Labora- 
tories in the field of nuclear rockets was 
not apparently cleared for publication. It 
is to be hoped that within a year or so a 
second edition of this book can be pro- 
duced, the contents of which will more 
closely conform to the title of the book. 


M. H. Smith, Associate Editor, and M. H. Fisher, Contributor 
The James Forrestal Research Center, Princeton University 


Jet Propulsion Engines 


TRIDAC Flight Simulator for Missiles, 
by J. J. Gait, Pusées et Recherche Aeron., 
vol. 3, April 1958, pp. 10-19. (In French.) 

British Rockets and Missiles, by J. 
Pellandini, Fusées et Recherche Aeron., vol. 
3, April 1958, pp. 29-33. (In French.) 

A Theoretical and Experimental Inves- 
tigation of the Aerothermopressor Process, 
by Alve J. Erickson, MIT’, Dept. Mech. 
Engng., Gas Turbine Lab., May 1958, 171 
pp. 

A Simplified Method for Estimating the 
Performance of Supersonic Ram-jets, by 
B. P. Beeton, Aeron. Quart., vol. 9, Aug. 
1958, pp. 195-212. 

Theoretical Investigation on the Fluido- 
dynamics of ‘‘Wave-engines’’ Operating 
According to the Principle of the ‘‘Com- 
prex,’? by Dino Dini, L’ Aerotecnica, vol. 
38, Feb. 1958, pp. 8-18. (In Italian.) 

DeHavilland Dual Spectra Rocket, 
Aviation Week, vol. 69, Aug. 25, 1958, p. 
25. 

Engineering Method of Ram-jet Thrust 
Determination Based on Experimentally 
Obtained Combustor Parameters, by H. 
Rudolph Dettwyler and Maxime A. Faget, 
NACA Res. Mem. L53E21, Aug. 1953, 44 
pp., diagrs., tab. (Declassified from Con- 
fidential by authority of NACA Res. 
Abstr. 128, 7/28/58.) 

Military Studies Rocket Test Chamber, 
by Michael Yaffee, Aviation Week, vol. 
69, Aug. 4, 1958, pp. 51, 53, 55. 


Aerodynamics of Jet 
Propelled Missiles 


Effects of Boundary-layer Displacement 
and Leading-edge Bluntness on Pressure 
Distribution, Skin Friction, and Heat 
Transfer of Bodies at Hypersonic Speeds, 
by Mitchel H. Bertram and Arthur Hen- 


Eprror’s Note: Contributions from Pro- 
fessors E. R. G. Eckert, J. P. Hartnett, T. 
F. Irvine Jr. and P. J. Schneider of the 
Heat Transfer Laboratory, University of 
Minnesota, are gratefully acknowledged. 
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derson Jr., NACA TN 4301, July 1958, 33 
pp. 

Laminar Incompressible Boundary 
Layer Flow over the Cylindrical Portion of 
Slender Bodies of Revolution at Small 
Angles of Attack, by John C. Martin and 
Nathan Gerber, Aberdeen Proving Ground, 
Ball. Res. Labs., Rep. 1038, Nov. 1957, 54 
pp. 

An Experimental Investigation of the 
Stability of the Hypersonic Laminar 
Boundary Layer, by Anthony Demetri- 
ades, Calif. Inst. Tech., Guggenheim 
Aeron. Lab., Hypersonic Res. Proj., Mem. 
43, May 1958, 85 pp. 

An Experimental Investigation of the 
Turbulent Boundary Layer behind a For- 
ward Facing Step in Supersonic Flow, by 
Dave Robert Moore, Univ. of Texas, De- 
fense Res. Lab. CF 2707, (DRL-425), 
April 1958, 27 pp., 3 tabs., 25 figs. 

Dynamic Stability of Vehicles Travers- 
ing Ascending or Descending Paths 
through the Atmosphere, by Murray To- 
bak and H. Julian Allen, NACA TN 4275, 
July 1958, 35 pp. 

Inviscid Hypersonic Flow Near the 
Stagnation Point of Oblate Ellipsoidal 
Noses, by Marcel Vinokur, J. Aero/- 
Space Sci., vol. 25, July 1958, pp. 469- 
470. 

Flight Measurements of Average Skin- 
friction Coefficients on a Parabolic Body 
of Revolution (NACA RM-10) at Mach 
Numbers from 1.0 to 3.7, by J. Dan Lopo- 
ser and Charles B. Rumsey, NACA Res. 
Mem. L54G14, Aug. 1954, 32 pp. (De- 
classified from Confidential by authority 
of NACA Res. Abstr. 129, 7/17/58.) 

Flow Around a Cone at Angle of Attack 
at High Supersonic Speeds, by A. L. 
Gonor, Akademii Nauk SSSR, Tzvestiia, 
Otdelenie Tekhnicheskikh Nauk, no. 7, July 
1958, pp. 102-105. (In Russian. ) 

An Estimate of the Fluctuating Surface 
Pressures Encountered in the Reentry of a 
Ballistic Missile, by Edmund E. Calla- 
ghan, NACA TWN 4315, July 1958, 18 pp., 
diagrs. 

Wind Tunnel Investigation at Mach 
Numbers of 2.0 and 2.9 of Several Con- 
figurations of a Supersonic Ram-jet Test 


Our Advanced Design Group 
is at work on 


NUCLEAR 
ROCKET ENGINES 


and needs the help of 
the following men 


Senior Engineer. MS or PhD. He 
will study nuclear engine applications, 
determine the most suitable missions, 
correlate engine thrust requirements 
with vehicle performance, and evalu- 
ate vehicle and engine system design 
problems peculiar to the thermal, 
nuclear, and space environment of 
nuclear missile systems. 


Senior Engineer. MS or PhD. He 
will make system dynamic studies of 
nuclear rocket engine systems to 
determine system stability and opera- 
tional characteristics, establish con- 
trol system requirements to obtain 
required accuracy and response char- 
acteristics, and synthesize new meth- 
eds of control required by thermal 
and nuclear radiation environment. 


Senior Design Engineers or Spe- 
cialists. BA or MS or equivalent. 
These men should be conceptual 
rocket engine designers with experi- 
ence in turbopump or high speed 
rotating machinery layouts and famil- 
iarity with combustion devices. They 
should be creative designers, oriented 
toward preliminary design. 


Please write to Mr. A. D. Jamieson, 
Engineering Personnel Dept., 6633 
Canoga Avenue, Canoga Park, Calif. 


ROCKETDYNE F2 


A DIVISION OF NORTH AMERICAN AVIATION, INC. 
FIRST WITH POWER FOR OUTER SPACE 
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AEROJET... 
for engineering careers 


STRESS ANALYSTS 
1. Experience and/or training in 
analytical stress computations, with 
emphasis on structural and pressure 
vessel applications. B.S. required; 
graduate work desirable. 


2. Experience and/or training in 
experimental stress analysis. Familiar- 
ity with advanced experimental tech- 
niques and expedient methods of stress 
analysis. B.S. required; graduate work 
desirable. 


GAS DYNAMICIST 

M.S. or Ph.D. in engineering, physics 
or applied mechanics. Five to ten years’ 
experience in analytical and experimen- 
tal work in the fields of boundary layer 
studies, supersonic and hypersonic 
flows, or application of non-ideal-gas 
theories. 


AERODYNAMICIST 

M.S. or Ph.D. in aeronautical or 
mechanical engineering or applied 
mechanics. Five to ten years’ experi- 
ence in the field of advanced aero- 
thermodynamics, missile aero- 
dynamics, heat transfer in rockets, 
missiles, etc. 


MISSILE SYSTEMS ANALYST 
B.S., M.S. or Ph.D. in engineering, 
physics or applied mechanics. Five to 
ten years’ experience in analysis and 
evaluation of interactions between pro- 
pellant systems and other missile com- 
ponents. Capable of developing 
optimum vehicle configurations from 
the standpoint of the propulsion sys- 
tem. 

U.S. Citizenship Required 

Resumes cordially invited. Write: 

E. P. JAMES 
AEROJET-GENERAL CORPORATION 
P.O. BOX 1947 
SACRAMENTO, CALIFORNIA 
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Vehicle, by J. Richard Spahr and Robert 
A. Robinson, NACA Res Mem. A50C20, 
Aug. 18, 1950, 61 pp. (Declassified by 
authority of NACA Res, Abstr. 126, 6/24/ 
58.) 


Heat Transfer and Fluid 
Flow 


Flow of Conducting Fluids Past Magne- 
tized Bodies, by A. Kulikovskii, Doklady 
Akad. Nauk SSSR, vol. 117, no. 2, Nov. 
11, 1957, pp. 199-202. (In Russian.) 

Some Investigations on Arc Discharges 
in Helium and Hydrogen at High Pres- 
sures, by A. Michels, H. Tjin, A. Djie, H. 
DeKluiver and C. A. Ten Seldam, Phys- 
ica, vol. 23, no. 12, Dec. 1957, pp. 1115- 
1125. 

Plasma Jet: Research at 25,000 F, 
by J. W. Reid, Machine Design, vol. 30, 
no. 3, Feb. 6, 1958, pp. 22-24. 

On the Visualization of Laminar Bound- 
ary Layer Oscillations and the Transition 
to Turbulent Flow, by H. Bergh and B. 
Van Den Berg, Zeitschrift fur Angewandte 
Mathematik und Physik, vol. 9, 1958, pp. 
97-104. 

Influence of Free-Stream Turbulence 
on Heat Transfer by Convection from 
Isolated Region of Plane Surface in Par- 
allel Air Flow, by A. Edwards and B. N. 
Furber, Proc. Inst. Mech. Engrs., vol. 
170, no. 28, 1956, pp. 941-954. 

Heat Transfer in a Laminar Boundary 
Layer with Constant Fluid Properties and 
Constant Wall Temperature, by A. Smith 
and D. Spalding, J. Roy. Aeron. Soc., vol. 
62, Jan. 1958, pp. 60-64. 

Heat Transfer from Surfaces of Non- 
uniform Temperature, by D. Spalding, 
- Fluid Mech., vol. 4, May 1958, pp. 22- 
32. 


A Boundary-layer Treatment of Lam- 
inar-film Condensation, by E. Sparrow 
and J. Gregg, ASME Paper 58-SA-2, 
June 1958. 

On the Aerodynamic Heating of Blunt 
Bodies, by E. Van Driest, ZAMP, vol. 
9, 1958 pp. 233-248. 

Prediction of Flashing Water Flow 
through Fine Annular Clearances, by A. 
Agostinelli and V. Salemann, ASME 
Paper 57—A-78, Dec. 1-6, 1957, 5 pp. 

Heat-transfer and Pressure-drop Char- 
acteristic of Internal Finned Tubes, by E. 
C. Brouillette, T. R. Mifflin and J. E. 
Myers, ASME Paper 57—A-—47, Dec. 1-6, 
1957, 9 pp. 

Subcritical and Critical Flow through 
Strain-through Elbow, and Tee A-N 
Fittings and Sharp-edged Orifices at 
Elevated Temperatures, by A. L. Ducoffe, 
J. R. Bennett and C. G. Ray, ASME 
Paper 57—A-60, Dec. 1-6, 1957, 8 pp. 

Local Laminar Heat Transfer in Wedge- 
shaped Passages, by E. R. G. Eckert, T. 
F. Irvine Jr. and J. T. Yen, ASME 
Paper 57-A-133, Dec. 1-6, 1957, 6 pp. 

Effect of Heat Transfer on Flow Field at 
Low Reynolds Numbers in Vertical Tubes, 
by T. J. Hanratty, E. M. Rosen and R. L. 
Kabel, Ind. & Engng. Chem., vol. 50, no. 5, 
May 1958, pp. 815-820. 

Heat Transfer in Fully Developed Flow 
between Parallel Plates with Variable 
Heat Sources, by A. L. Loeffler Jr., Nu- 
clear Sci. & Engng., vol. 2, no. 5, Sept. 
1957, pp. 547-566. 

Non-Newtonian Flow, by W. L. Wilkin- 
son, Indust. Chemist, vol. 34, Feb. 1958, 
pp. 79-84. 

Investigation of Effective Thermal Con- 
ductivities of Powders in Various Gases, 


by R. G. Deissler and J. S. Boegli, ASME 
Paper 57-A-110, Dec. 1-6, 1957, 7 pp. 

Prismatic Fin with Non-linear Heat Loss 
Analyzed by Resistance Network and 
Iterative Analogue Computer, by M. 
Hutcheon and D. Spalding, Brit. J. 
Appl. Phys., vol. 9, 1958, pp. 185-191. 

Heat Transfer from Rods Normal to 
Subcooled Water Flow for Non-boiling 
and Surface-boiling Conditions up to and 
including Burnout, by 8S. P. Kezios and 
R. K. Lo, Argonne Nat. Lab. (U.S. AEC) 
ANL-5822, Jan. 1958, 12 pp. 

On the Uniqueness of the Solution of 
Stephan’s One-dimensional Problem for 
a Heat Conducting Medium in One-phase 
Initial State, by L. I. Rubinshtein, Dov- 
lady Akad. Nauk SSSR, vol. 117, no. , 
Nov. 21, 1957, pp. 387-390. (In Ru-- 
sian.) 

Heat and Mass Exchange in Forced 
Flow around Bodies of Various Shapes, by 
O. Krischer and G. Loos, Chemie-Ingenieu - 
Technik, vol. 30, no. 2, Feb. 1958, pp. 69-7.:. 
(In German. ) 

Mass Transfer in Liquid Metal Systems, 
II: Isothermal Transfer, by J. W. Tayloi, 
and A. G. Ward, Nuclear Power, vol. 3, no. 
23, March 1958, pp. 101-105. 

Measurement of Intense Beams of 
Thermal Radiation, by A. Broida and A. 
B. Willoughby, J. Optical Soc. Amer, 
vol. 48, no. 5, May 1958, pp. 344-350. 

Some High-pressure, High-temperature 
Apparatus Design Considerations: Equip- 
ment for Use at 100,000 Atmospheres and 
3000°C, by H. T. Hall, Rev. Sez. Instr., 
vol. 29, no. 4, April 1958, pp. 267-275. 

Apparatus for Precision Flash Radiog- 
raphy of Shock and Detonation Waves in 
Gases, by H. T. Knight and D. Venable, 
Rev. Sci. Instr., vol. 29, no. 2, Feb. 1958, 
pp. 92-98. 


Combustion, Fuels and 
Propellants 


Combustion—an Aeronautical Science, 
by H. W. Emmons, (The 1958 Minta 
Martin Aeron. Lecture), Inst. Aeron. Sci., 
SMF Fund Paper FF19, 1958, 36 pp. 

Aerodynamics of Flame Stabilization; 
Experimental Procedures, by A. Sherman, 
Princeton Univ., Dept. Aeron. Engngqg., 
Rep. 426, July 1958, 23 pp. 

Effect of Pressure on Turbulent Burning 
Velocity, by B. Fine, Combustion and 
Flame, vol. 11, June 1958, pp. 109-116. 

Separation of the Products of Cool 
Flame Oxidation of Propane, by C. R. 
Yokley and R. E. Ferguson. Combustion 
and Flame, vol. 11, June 1958, pp. 117- 
128. 

The Detonability Envelope of Gaseous 
Fuming Nitric Acid-n-Heptane-Nitrogen 
Mixtures as a Function of Initial Pressure, 
by R. L. Potter and D. H. Wayman, 
Combustion and Flame, vol. 11, June 1958, 
pp. 129-136. 

The Hydrazine Flame, by M. Gilbert, 
Combustion and Flame, vol. 11, June 1958, 
pp. 137-148. 

Kinetics of Hydrazine Decomposition 
in a Laminar Non-isothermal Flow, by M. 
Gilbert, Combustion and Flame, vol. 11, 
June 1958, pp. 149-156. 

Gaseous Nitric Acid - Hydrocarbon 
Flames, by J. Mertens and R. L. Potter, 
Combustion and Flame, vol. 11, June 1958, 
pp. 181-192. 

Some Recent Developments Concerning 
Free Radicals (A Review), by G. J. Min- 
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CHEMISTRY 


Expanding the 


Research and Development facilities 
in the Stanford Industrial Park at 


Frontie rs of Palo Alto, California, provide the 


latest in technical equipment. 


Space Technology 


(above) Synthesizing high-energy 
propellant in development of new 
and exotic fuels in propulsion 
chemistry laboratory. 


(left) Vacuum systems used in 
study of gas metal reactions at 
elevated temperatures — one of 
many studies in materials 
research program. 


At Lockheed, chemistry plays a highly important cells, and explosives. In electrochemistry, solid state 
role in materials, research, metallurgy, solid state work is directed to related fuel research in electrode 
electronics, ordnance, propulsion and test. Appli- reaction kinetics and materials synthesis by electro- 
cations are for advanced missile, spacecraft and chemical methods and vacuum deposition. 
nuclear projects. The test laboratories assay raw materials, eval- 
Special programs include: research in high and uate materials and components and assist in design 
low temperature materials; radio-active and toxic and development. 
materials; pure metals; properties of beryllium and Scientists and engineers of outstanding talent 
the alloys and the development of refractory and and inquiring mind are invited to join us in the 
special-purpose metals, cermets and ceramics. nation’s most interesting and challenging basic 
Other chemical studies embrace reaction kinet- research programs. Write: Research and Develop- 
ics and thermodynamics, and the development of ment Staff, Dept. D-25, 962 W. El Camino Real, 
new and exotic fuels, high-energy batteries and fuel Sunnyvale, California. 


“The organization that contributed most in the past year to the advancement of the 
art of missiles and astronautics.” | NATIONAL MISSILE INDUSTRY CONFERENCE AWARD 


SUNNYVALE, PALO ALTO, VAN NUYS, 


Lockheed SANTA CRUZ, SANTA MARIA, CALIFORNIA 
MISSILES AND SPACE DIVISION CAPE CANAVERAL, FLORIDA 


ALAMOGORDO, NEW MEXICO 
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koff, Combustion and Flame, vol. 11, June 
1958, pp. 193-207. 

Theory of Detonations, I: Irreversible 
Unimolecular Reaction, by Joseph 0. 
Hirschfelder and Charles F. Curtiss, //. 
Chem. Phys., vol. 28, June 1958, pp. 1130- 
1146. 

Theory of Detonations, II: Reversible 
Unimolecular Reaction, by B. Linder, (. 
F. Curtiss and J. O. Hirschfelder, //. 
Chem. Phys., vol. 28, June 1958, pp. 1147- 
1151. 

Specific Heat of Liquid Ozone, ly 
Robert I. Brabets and Thomas E. Wate:- 
man, J. Chem. Phys., vol. 28, June 195s, 
pp. 1212-1214. 

Determination of Hydrogen Atoms in 
Rich, Flat, Premixed Flames by Reaction 
with Heavy Water, by C. P. Fenimore an | 
G. W. Jones, J. Phys. Chem., vol. 62, June 
1958, pp. 693-697. 

Approximate Solutions of Transient and 
Two-dimensional Flame Phenomena: 
Constant-enthalpy Flames, by D. 1: 
Spalding, Proc. Roy. Soc., Series A., vo. 
245, June 24, 1958, pp. 352-372. 

Experimental Thermal Conductivities 
of the N,0,—2NO, System, by K. P. 
Coffin and C. O'Neal Jr.. NACA TN 
4209, Feb. 1958, 22 pp. 

Spontaneous Ignition of Pentaborane 
Sprays in a Hot-Air Stream, by Erwin A 
Lezberg and Albert M. Lord, NACA Res. 
Mem. E55F 29a, Nov. 1956, 26 pp., diagrs 
(Declassified by authority of NACA Res. 
Abstr. 128, 6/24/58.) 

The Heats of Forma- 
tion of Hexacyclo (7:2:1:07'> 
DRINKING FAUCET co. :0°"2) dodecane; Two Techniques for 

. lids, by Ward N. Hubbard, F. R. Frow 
1443 FOURTH STREET - BERKELEY 10, CALIFORNIA Chen. 
—+_————— vol. 62, July 1958, pp. 821-823. 
The Differential Thermal Analysis of 
: Perchlorates, II: The System LiCI1O,- 
SOUTHWEST — LiNO,, by Meyer M. Markowitz, J. 
Phys. Chem., vol. 62, July 1958, pp. 827- 
829. 

Measurement of Flame Speeds by a 
Nozzle Burner Method, by Carl Halpern. 
J. Res., Nat. Bur .Standards, vol. 60, June 
1958, pp. 535-546. 

The Initiation of Explosion in Solid 
Secondary Explosives, by G. B. Cook, 
Proc. Roy. Soc., London, vol. A 246, July 
29, 1958, pp. 154-159. 


RID THE BODY OF 
CAUSTICS and CHEMICALS 
Instantly! Thoroughly! 


ACCIDENTS with caustic chemicals 

strike with shocking swiftness—and | 

<4 Haws Drench Showers are instant- | 
oe ly ready to deliver relief just.as fast! | 
He: A solid downpour washes away de- | 
structive materials—saving seconds | 

until medical aid arrives, possibly | 

averting serious injury and exces- © 

sive compensation claims. Haws | 

Drench Showers can help you! Write | 

for details and illustrated literature. | 


MODEL 8935— Drench shower augmented by 


Haws eye-wash fountain. A complete safety 
station —always ready! 


J The Initiation of Liquid Explosives by 
Shock and the Importance of Liquid 
7 Be - Break-up, by G. H. Johansson and others, 
ae =. Proc. Roy. Soc., London, vol. A 246, July 
PATENTED U.S.A. 29, 1958, pp. 160-166. 
e World Rights Reserved 
$ CHARACTERISTICS 
: e ANALYSIS RECOMMENDED USE Materials of Construction | 
: For types operating under high | 
Stointess Stel Ball ond {temperature (800-1200 degrees F.). Compatibility of Metals with Liquid | 
Chrome Alloy Steel Ball For types opera 13080 695 00 at  / Fluorine at High Pressures and Flow 
and Race ultimate loads 
For types operating under norm 1 Velocities, by Harold W. Schmidt, NACA 
Steel Ball { with minimum friction requirements. Res. Mem. E58D11, July 1958, 15 pp. | 
Compressive Strength and Creep of 
Thousands in use. Backed by years of service life. Wide variety of 17-7 PH Stainless-steel Plates at Elevated 
Plain Types in bore sizes 3/16” to 6” Dia. Rod end types in similar Temperatures, by Bland A. Stein, NACA 
size range with externally or internally threaded shanks. Our Engi- TN 4296, July 1958, 33 pp. 
neers welcome an opportunity of studying individual requirements A Review of Seal Materials for Guided i 
and prescribing a type or types which will serve under your demand- Missile Applications, by Robert J. Berg, 
ing conditions. Southwest can design special types to fit individual Calif. Inst. Tech., Jet Prop. Lab., Progr. 


Rep. 20-340, March 1958, 16 pp. 

Dynamic Response of Cylindrical Tanks, 
by F. L. DiMaggio, Columbia Univ., Dept. 
1-58 XZ 554 i), rmed Forces 
SOUTHWEST PRODUCTS CO. Special Weapons Proj. 1075), May 1958, 


1705 SO. MOUNTAIN AVE., MONROVIA, CALIFORNIA 24 pp., 17 figs 
Plastics in the Space Age, by John H. 


specifications. As a result of thorough study of different operating 
conditions, various steel alloys have been used to meet specific 
needs. Write for Engineering Manual No. 551.Address Dept. JP-59. 
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Lux, J. Franklin Inst., vol. 266, July 
1958, pp. 21-26. 

Wires and Cables: Vital to Missiles, 
by Peer Fossen, Raymond M. Nolan, 
Missiles and Rockets, vol. 4, Sept. 8, 
1958, pp. 23-25. 

Magnesium: Will It Find Biggest 
Market in Missiles, by Carl N. Morten- 
son, Missiles and Rockets, vol. 4, Sept. 8, 
1958, pp. 46-49. 

Effect of Temperature on Dynamic 
Modulus of Elasticity of Some Structural 
Alloys, by Louis F. Vosteen, NACA TN 
4348, Aug. 1958, 19 pp., diagrs., tabs. 

Influence of Not-working Conditions on 
High-temperature Properties of a Heat- 
resistant Alloy, by John F. Ewing and J. 
W. Freeman, NACA Rep. 1341, 1957, 42 
pp. (Formerly 3727.) 

Surface Tension at Elevated Tempera- 
tures, IV: Surface Tension of Fe-Se and 
Fe-Te Alloys, by W. D. Kingery, J. 
Pi ys. Chem., vol. 62, July 1958, pp. 878- 
879. 


Instrumentation, Data 
Recording, Telemetering 


Theory of Construction of High-speed 
Automatic Regulating Systems, by M. V. 
Meerov, Automatika i Telemekhanika, vol. 
19, no. 7, July 1958, pp. 621-633. (In 
Russian. ) 

On the Speediest Transient Process in 
Servo Systems with Power, Torque and 
Rate Limitations of the Power Unit, by E. 
A. Rosenman, Automatika i Telemekhan- 
ika, vol. 19, no. 7, July 1958, pp. 633-653. 
(In Russian. ) 

New Principles for the Design of Tel- 
emetering Systems with Time and Width 
Pulse Modulation, by B. A. II’in and A. I. 
Novikov, Automatika i Telemekhanika, 
vol. 19, no. 8, 1958, pp. 757-761. (In 
tussian. ) 

Noise Stability of Telemetering Signal 
Transmission along a Channel, with 
Fluctuation Noise, by V. A. Kashiriu and 
G. A. Shastova, Automatika i Telemek- 
hanika, vol. 19, no. 8, 1958, pp. 762-775. 
(In Russian. ) 

A Constant Liquid Flow Device, by G. 
H. Laycock, Brit. J. Appl. Phys., vol. 9, 
Aug. 1958, pp. 333-335. 

Pressure Measurement in Armament 
Research, by A. J. Cable, NATO, AGARD, 
Rep. 177, March 1958, 10 pp., 6 figs. 

Simple Surface Thermometer, by C. F. 
Rothe, Rev. Sci. Instr., vol. 29, no. 5, May 
1958, pp. 436-437. 

Experimental Studies on the Initiation 
of Electric Arcs, by T. Yamaguti, J. Phys. 
Soc. Japan, vol. 13, no. 5, May 1958, pp. 
520-525. 

A Method for the Measurement of Heat 
Conductivity at 1100°C, by K. Schafer and 
F. W. Reiter, Zettschr. f. Elektrochemie, 
vol. 61, no. 9, 1957, pp. 1230-1235. (In 
German. ) 

Stabilization of Energy Transfer at Long 
Distances, by Rene Perret, France, Minis- 
tere de l’ Air, Publications Sci. et Tech., 
Note Tech. 74, 1958, 104 pp. (In French.) 

A Thermistorized Apparatus for Dif- 
ferential Thermal Analyisis, by Jack M. 
Pakulak Jr. and Guy William Leonard, 
NAVORD Rep. 5829 (NOT'S 1920), June 
1958, 15 pp. 

Telemetry Looks into the Future, 
Aviation Age, vol. 30, Aug. 1958, pp. 20- 
24. 


Design Digest: Missile Telemeter 
1959 


The challenge of outer space has created the need for a man. Space flight 
is his business. Recognition is his goal. He'll find it at Martin. He will 
make decisions, move up as quickly as his abilities take him. If your inter- 
ests lie in Astrophysics, Space Mechanics, Systems Dynamics, Magneto 
Hydrodynamics or Aerophysics, this man could be you. Write or phone 
William Spangler, Manager, Professional Employment, Dept.A-4, The Martin 
Company, Baltimore 3, Md. 
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Nuclear 
Scientists 


Steadily advancing 
programs at 
this leading Nuclear 
have created 
several entirely new 
positions on the 
professional staff. 


PhD PHYSICIST 
Salary range: $10,000 to $15,000 


—to develop theory, ana- 
lyze phenomena and pre- 
dict performance and char- 
acteristics of nuclear re- 
actors, and/or recommend 
experimental research pro- 
grams. 


PhD MATHEMATICIAN 


Salary range: $10,000 to $15,000 


—to deal with advanced 
problems arising in the de- 
sign of nuclear reactors and 
powerplants. 


PhD METALLURGIST 


Salary range: $10,000 to $15,000 


—to initiate and evaluate 
experiments in materials 
and process development; 
study irradiation effects. 


Inquiries will be held in 
confidence. Forward com- — 
plete résumé, including sal- 
ary requirement, to: 


Box 1001 
AMERICAN ROCKET SOCIETY 
500 Fifth Ave. 

New York 36, New York 


Mates PDM-FM with FM-PM, Aviation 
Age, vol. 30, Aug. 1958, pp. 126-127. 

New Missiles and Spacecraft Challenge 
Telemetry Technology, Aviation Age, vol. 
30, Aug. 1958, pp. 128-133. 

Lockheed Switches to PAM-FM, Avia- 
tion Age, vol. 30, Aug. 1958, p. 138. 

Light Modulator Records Airborne 
Radar Displays, by Leo Levi, Electronics, 
vol. 31, Aug. 1, 1958, pp. 80-83. 


Guidance Systems and 
Components 


Rocket Guidance as a Problem of En- 
slaving, by E. Perret, Fusées et Recherche 
Aeron., vol. 3, April 1958, pp. 21-25. 
(In French.) 


Servomechanisms and 
Controls 


Analytical Method of Synthesis of 
Linear Control Systems in the Presence of 
Noise, when Dynamic Precision is Speci- 
fied, by K. I. Kurakin, Automatika i Tele- 
mekhanika, vol. 19, no. 5, 1958, pp. 408- 
417. (In Russian.) 

On Improving the Quality of On-Off 
Control, by A. I. Cherepanov, Automatika 
t Telemekhanika, vol. 19, no. 5, 1958, pp. 
471-474. (In Russian.) 

Stable Receiving Circuits for Remote 
Control, by S. J. Neshyba _— F. E. 
Brooks Jr., Electronics, vol. 31, Aug. 1, 
1958, pp. 74-76. 


Flight Vehicle Design and 
Testing 


Stationary Conditions for Problems 
Involving Time Associated with Vertical 
Rocket Trajectories, by Angelo Miele, J. 
Aero/Space Sci., vol. 25, July 1958, pp. 
467-469. 

Guided Weapons and Aeronautics, by 
R. Cockburn, J. Roy. Aeron. Soc., vol. 62, 
Aug. 1958, pp. 539-558. 

Solving Missile Problems, Ordnance, 
vol. 43, July-Aug. 1958, p. 149. 

Talos Land-based System, by D. B. 
Holmes, Ordnance, vol. 43, July-Aug. 
1958, p. 156. 


Space Flight 


Study of the Vital Activity of Animals in 
Flight in Rockets in the Upper Atmosphere, 
by Pokrovski, Fusées et Recherche Aeron., 
vol. 3, April 1958, pp. 5-9. (In French. ) 

Regulation of Sovereignty in the Very 
High Atmosphere and the Law of Space 
Ships, by E. B. Winandy, Fusées et 
Recherche Aeron., vol. 3, April 1958, pp. 
34-39. (In French.) 

Problems Posed by the Law of Astro- 
nautics and Atomic Power in Space Flight, 
by E. B. Winandy, Fusées et Recherche 
Aeron., vol. 3, April 1958, pp. 40-43. 
(In French.) 

Air Drag Effect on a Satellite Orbit De- 
scribed by Difference Equations in the 
Revolution Number, by R. E. Roberson, 
Quart. Appl. Math., vol. 16, July 1958, 
pp. 131-136. 

Earth-Moon Rocket Trajectories, by 
Louis Gold, J. Franklin Inst., vol. 266, 
July 1958, pp. 1-8. 

The Human System in Space, by C. 
Wright Reininger, Missiles and Rockets, 
vol. 4, Sept. 8, 1958, pp. 33-34. 


Soviets Certain They’ll Be First on 
Moon, by Dr. Victor P. Petrov, Missiles 
and Rockets, vol. 4, Sept. 8, 1958, pp. 50, 
55-58. 

An Analytical Solution for Flight Time 
of Satellites in Eccentiric and Circular 
Orbits, by F. M. Perkins, Astronautica 
Acta, vol. 4, no. 2, 1958, pp. 113-134. 

Orbit Determination from the Vector of 
Orbit Velocity, by M. Vertregt, Astro- 
nautica Acta, vol. 4, no. 2, 1958, pp. 135- 
137. (In German.) 

Moon Refueling for Interplanetary 
Vehicles, Aviation Age, vol. 30, Aug. 
1958, pp. 22-24. 

Design Progress: Missile-launched 
Manned Space Capsule, Aviation Ag, 
vol. 30, Aug. 1958, pp. 172-173. 

How to Find Thermal Equilibrium in 
Space, Aviation Age, vol. 30, Aug. 195%, 
pp. 174-176. 

Army Explorer IV Satellite Reports 
Corpuscular Radiation Intensity, by Eve: 
Clark, Aviation Week, vol. 69, Aug. 4, 
1958, pp. 30-31. 

Ground Testing May Ease Space Prob- 
lems, by Richard Sweeney, Aviatio, 
Week, vol. 69, Aug. 4, 1958, pp. 59, 61. 


Aerophysics, Astrophysics 


High-altitude Atmospheric Density, b) 
Theodore E. Sterne, Phys. Fluids, vol. 1 
May-June, 1958, pp. 165-170. 

Interaction of the —, Wind with the 
Geomagnetic Field, by E. N. Parker, Phys 
Fluids, vol. 1, May -June 1958, pp. 171 
187. 

The Statistics of Meteors in the Earth’s 
Atmosphere, by Gerald 8S. Hawkins and 


Richard B. Southworth, Smithsonian 
Contrib. to Astrophys., vol. 2, no. 11, 1958, 
pp. 349-364. 


Granulation and Oscillations of the Solar 
Atmosphere, by Charles Whitney, Smith- 
sonian Contrib. to Astrophys., vol. 2, no. 
12, 1958, pp. 365-376. 

A Cosmic Ray Increase Related to Solar 
Activity, by J. Katzman, Canadian J. 
Phys., vol. 36, July 1958, pp. 807-815. 

Analysis of Meteoric Body Doppler 
Radar Records Taken during a Geminid 
Shower Period, by M. Srimama Rao, 
Canadian J. Phys., vol. 36, July 1958, pp. 
840-848. 

On the Fluorescence of Air, Excited by 
Fast Electrons: Light Yield as a Function 
of Pressure, by A. FE. Griin, Canadian J. 
Phys., vol. 36, July 1958, pp. 858-870. 

Photon Correlation in Coherent Light 
Beams, by E. Brannen, H. I. 8. Ferguson 
and W. Wehlau, Canadian J. Phys., vol. 
36, July 1958, pp. 871-879. 

Cosmic-ray Modulation by Solar Wind, 
by E. N. Parker, Phys. Rev., vol. 110, 
June 15, 1958, pp. 1445-1448. 

Nature of the Current Reduction in the 
Primary Cosmic-ray Intensity, by M. A. 
Pomerantz, 8S. P. Agarwal and V. R. Pot- 
nis, Phys. Rev. Letters, vol. 1, no. 2, July 
15, 1958, p. 65. 

Large-amplitude Hydromagnetic Waves 
Above the Ionosphere, by A. J. Dessler, 
Phys. Rev. Letters, vol. 1, no. 2, July 15, 
1958, p. 68. 


Atomic Energy 


New Batteries for the Space Age, by D. 
Linden and A. F. Daniel, Electronics, 
Engng. Ed., vol. 31, July 18, 1958, pp. 59- 

Stability of the Pinch, by S. Chandra- 
sekhar, A. N. Kaufman and K. M. Watson. 
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Proc. Roy. Soc., London, vol. A 245, July 
8, 1958, pp. 485-455. 

Nuclear Fuels for Hydromagnetic Shock 
Engine, by C. F. Johnson, Aviation Age, 
vol. 30, Sept. 1958, pp. 118-124. 

Stabilized Pinch and Controlled Fusion 
Power, by S. A. Colgate and H. P. Furth, 
Science, vol. 128, Aug 15, 1958, pp. 337- 
343, 29 refs. 

The Ion Rocket Engine, by R. H. Boden, 
\. A. Aviation, Rocketdyne Div., Rep. R- 
645, Aug. 26, 1957, 40 pp. (AFOSR TN- 
57,573; ASTIA AD 136,558. ) 

Outer Space Propulsion by Nuclear 
Energy: Hearings before Subcommittees 

. Jan. 22, 23, Feb. 6, 1958, Congress, 
8 th, 2nd Session, Joint Comm. on Atomic 
Energy, Ww ashington, Govt. Print. Off., 

1158, 232 pp. 

Neutron Generation from Straight 
Pinches, by R. E. Dunway and J. A. 
Phillips. J. Appl. Phys., vol. 29, Aug. 
1958, pp. 1137-1143. 

Supersonic Motion of Vacuum Spark 
oe along Magnetic Fields, by D. 

Finkelstein, G. A. Sawyer and T. F. Strat- 
ton, Phys. Fluids, vol. 1, May-June 1958, 
pp. 188-192. 

Effect of Charge Separation on Plasma 
Diffusion in a Strong Magnetic Field, by 
A. N. Kaufman, Phys. Fluids. vol. 1, May- 
June 1958, p. 252. 

Containment of a Fully-ionised Plasma 
by Radio-frequency Fields, by H. A. H. 
Hoot, S. A. Self, and R. B. R. Shersby- 
Harvie, SERL Tech. J., vol. 8, no. 3, May 
1958, pp. 107-131. (Gt. Brit. Services 
Electronics Res. Lab.) 

The Thermodynamic Properties of 
Hydrogen and Water as Possible Working 
Fluids for Nuclear Rockets, by I. Sanger- 
Bredt, Forschungsinstitut fiir Physik der 
Strahlantriebe, Stuttgart, Mitt. 16, May 
1958, 58 pp. (In English.) 

Stochastic Theory of Transport Phe- 
nomena in a Reacting Plasma at Extreme 
Temperatures, by H. J. Kaeppler, For- 
f. Phys. der Strahlantriebe, 
Stuttgart, Mitt. 15, March 1958, 115 pp., 
(In E nglish. ) 


Russian Technical 
Articles* 


Perturbation Method for Problems of 
Strong Detonations, by E. I. Andrianikin, 
levestiia Akad. Nauk SSSR, Otdelenie 
Tekhnicheskikh Nauk, no. 12, Dec. 1958, 
pp. 5-14. (In Russian.) 

Absorption and Dispersion of Infinites- 
imal Force Plane Waves of Very Small or 
Very Large Frequency under the In- 
fluence of Radiative Heat Transfer, by 

. A. Prokoph’ev, Izvestiia Akad. Nauk 
SSSR, Otdel. Tekh. Nauk, no. 12, Dee. 
1958, pp. 15-23. (In Russian. ) 

On the Stability of Permanent Rotations 
of a Body Around a Fixed Point, by V. V. 
Rumyantsev, Prikladnaya Matematika i 
Vekhanika, vol. 21, no. 3, 1957, pp. 339- 
346. (In Russian.) 

Calculation of Some Gas Flows for Sonic 
Velocity, by P. I. Chuskin, Priklad. Mate- 


* The editors of Technical Literature 
Digest are making a systematic search for 
pertinent Russian articles, both in the orig- 
inal Russian and in translated form. For 
the balance of 1959, these will appear in a 
separate section with this heading, mainly 
to draw the readers’ attention to the new 
listing. After that, they will appear item 
by item under the proper subject headings. 


Aprit 1959 


MISSILE HARDWARE 


NEWBROOK 


Specializing in 
eMOTOR CASES 


Solid and Liquid 
Propellants 


JATO CASES 
© NOZZLES 4 
¢PLENUM CHAMBERS 
°BLAST TUBES 
© FUEL INJECTORS 


We have the 
“KNOW-HOW” 

We have developed new 
techniques, new methods, 
new processes that effect 
production economy so 
necessary to a successful 
missile program. 

Here at Newbrook you 
will find men with experi- 
ence gained from doing... 
a modern plant with up- 
to-date equipment... pre- 
cision inspection to meet 
your most exacting quality 
control requirements. 
And most important, New- 
brook specialization re- 
sults in strict reliability! Let 
us help you with your Mis- 
sile Hardware problems. 


Finest Welding Facilities 
Certified Welders 


Machining Nozzles 


X-Ray Inspection 


Below: Machining Motor Cases 


MACHINE CORPORATION 
45 MECHANIC ST. PHONE 45 SILVER CREEK, N. Y. 
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PRODUCERS OF 


turbines 


J69-T-29A 
1,700 Ib.-thrust model for 
drone applications. The 
new engine has 60% more 
thrust with only a 6% in- 
crease in weight. It is pres- 
ently powering the Ryan 
Q-2C target drone which 
recently underwent suc- 
cessful flight tests. 


J69-T-25 


Latest power plant for the 
Air Force T-37A twin-jet 
trainer manufactured by 
Cessna, the J69-T-25 
has increased thrust to 
1,025 Ibs. 


CAE has more than 10 years of turbine engine experience in the medium 
and low power range. The Air Force MA-1A trailer-mounted turbine 
air generator manufactured by CAE is on duty at air bases the world 
over, demonstrating outstanding performance and reliability. 


CONTINENTAL AVIATION & ENGINEERING CORPORATION 


12700 KERCHEVAL, DETROIT 15, MICHIGAN 


SUBSIDIARY OF CONTINENTAL MOTORS CORPORATION 


mat. t Mekhan., vol. 21, no. 3, 1957, pp. 
353-360. (In Russian.) 

Zempleri’s Theorem in Magnetohydro- 
dynamics, by S. V. Iordansky (translated 
from Doklady Akad. Nauk SSSR, vol. 121, 
Aug. 1, 1958, pp. 610-612). Morris D. 
Friedman, Inc., Needham Heights 94, 
Mass., Translation no. I-117. 

The Stability of an Elastic Plate in a 
Supersonic Flow, by Tung Ming-teh, 
Doklady Akad. Nauk SSSR, vol. 120, June 
1, 1958, pp. 726-729. (In Russian.) 

The Spread of Perturbations in Media 
with a Non-linear Dependence of Stresses 
on Strains and Temperatures, by G. J. 
Galin, Doklady Akad. Nauk SSSR, vo. 
120, June 1, 1958, pp. 734-738. (In 
Russian. ) 

Aerodynamics of High Speeds, Grazi- 
danskaia Aviatsiia, no. 6, June 1958, p} 
17-20. (In Russian.) 

Some Remarks on the High Speed 
Motion of Bodies in a Weak Gravitationa! 
Field, by K. P. Stanukovich, Doklady Akad 
Nauk SSSR., vol. 120, May 11, 1958, pp 
277-280. (In Russian.) 

Hypersonic Flow Around a Cone at ar 
Angle of Attack, by A. L. Gonor (trans: 
lated from IJzvestiia Akad. Nauk SSSR 
Otdel. Tekh. Nauk, no. 7, 1958, pp. 102- 
105). Morris D. Friedman, Inc., Need- 
ham Heights 94, Mass., Translation no 
G-154, 1958, 6 pp. 

Three-dimensional Streamline Flow 
Around Slender Bodies, by M. D. Khas- 
kind (translated from Priklad. Matemat. : 
Mekhan., vol. 20, no. 3, 1956, pp. 203-210) 
Morris D. Friedman, Inc., Needham 
Heights, 94, Mass., Translation no. K- 
165, 1958, 12 pp. 


On Minimum Drag Bodies at Hyper- 
sonic Speeds, by A. L. Gonor and G. G. 
Chernyi (translated from Jzvesttia Akad. 
Nauk SSSR Otdel. Tekh. Nauk, no. 7, 
1957, pp. 89-93). Morris D. Friedman, 
Inc., Needham Heights 94, Mass., Trans- 
lation no. G-141, 1958, 6 pp. 

Certain Peculiarities of the Flow around 
Bodies at High Supersonic Velocities, by 
G. L. Grodzovskii (translated from Jzves- 
tiia Akad. Nauk SSSR, Otdel. Tekh. Nauk, 
no. 6, 1957, pp. 86-92). Morris D. Fried- 
man, Inc., Needham Heights 94, Mass., 
Translation no. G-140, 1958, 11 pp. 


An Experimental Investigation of the 
Specific Heat of Water and Steam at 
High Parameters, by D. 8S. Rasskazov 
and A. E. Sheindlin, Doklady Akad. Nauk 
SSSR, vol. 120, June 1, 1958, pp. 771-774. 
(In Russian.) 

Motions with Homogeneous Deforma- 
tion in Magnetohydrodynamics, by A. G. 
Kulikovsky, Doklady Akad. Nauk SSSR, 
vol. 120, June 11, 1958, pp. 984-986. 
(In Russian. ) 

An Investigation of the Oxygen Dis- 
sociation Process Behind a Strong Shock 
Wave, by S. A. Losev, Doklady Akad. 
Nauk SSSR, vol. 120, July 21, 1958, pp. 
1291-1293. (In Russian.) 


Plane Problems in Magnetohydrody- 
namics, by G. 8. Golytzin, Zhur. Eksper. 
Teor. Fiziki, vol. 34, March 1958, pp. 
688-693. (In Russian.) 
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THERMODYNAMICISTS 
with ME, ChE or Nuclear Engineering Degrees (Bs, MS) 


There’s a new order of complexity for the 


in developing the nuclear propulsion system for aircraft 


The requirements for small-size, coupled 
with high-power density in a reactor for 
flight propulsion, raise problems in heat 
transfer never faced before. 


As an indication of the technical difficulty 
of these problems, consider the fact that 
the reactor now being developed by Gen- 
eral Electric for flight propulsion must 
generate more power than a submarine 
reactor, though it is only a fraction of the 
size. 


If problems “of this order” attract you, 
and you have had 3 to 8 years relative 
experience, your inquiry on any of the 
following assignments will be welcomed: 


@ EXPERIMENTAL FLUID DYNAMICS & THERMO- 
DYNAMICS ANALYSIS OF REACTOR PERFORMANCE. 


@ HEAT TRANSFER AND FLUID FLOW DESIGN CAL- 
CULATIONS FOR NUCLEAR REACTORS. 


@ CALCULATIONS ON NUCLEAR DESIGN PARA- 
METERS OF EXISTING OR PROPOSED REACTOR COM- 
PONENTS AS BASIS FOR CRITICAL EVALUATION. 


@ DEVELOPMENT TESTING OF HEAT TRANSFER & 
FLUID FLOW DESIGN OF REACTOR COMPONENTS. 


@ DEVELOPMENT TESTING OF REACTOR CORE 
COMPONENTS. 


@ MECHANICAL DEVELOPMENT OF POWER PLANT 
COMPONENTS. ALSO OF REACTOR STRUCTURAL 
COMPONENTS. (Includes devising and conduct- 
ing laboratory tests.) 


Assignments are at several levels 


Write in confidence, including salary requirement. 
Mr. P. W. Christos, Div. 35-MD 
AIRCRAFT NUCLEAR PROPULSION DEPARTMENT 


GENERAL @ ELECTRIC 


P.O. Box 132 


CINCINNATI 15, OHIO 
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CHEMIST S 


ANALYTICAL PhD or equivalent to perform research on 
new analytical methods and techniques. 

© INORGANIC—ORGANIC PhD or equivalent for synthetic 


PHYSICAL—PhD or ep Supervisory experience and back- 
e ground in experimental and ly and 

cybernetics required. To supervise research programs in 
e —, of physical and ballistic properties of liquid and solid 


PHYSICAL—B.S. with education or experience emphasis on 
thermodynamics, particularly heat capacity 


CHEMICAL ENGINEERS 


DEVELOPMENT ENGINEERS with experience design 
analysis, and operation of bench scale units and pilot tg 
© SENIOR PROCESS ENGINEER with broad rather than 
specialized pilot plant and production experience including 
e studies. Technical supervisory experience 
require 


OTHER 


APPLICATIONS RESEARCH CHEMIST or ENGINEER with 
industrial experience and interest in several technological fields. 
bo new industrial products or applications of existing 


PROPELLANT TESTING SPECIALIST—Chemist, Physicist 
or Engineer, B.S., M.S. or PhD with a minimum of 3 years’ 
e — ience in small scale testing and evaluation of liquid and 
propellants. To organize and coordinate programs for 
e st utilization of high energy fuels. 
Permanent positi with prod of Boron chemicals and high 
energy fuels located 25 miles north of Pittsburgh's Golden Triangle. 
Please send resume including salary requii ts in confid to: 


PERSONNEL MANAGER 


LLE: RY Callery 


CHEMICAL COMPANY Pennsylvania 


RESEARCH ENGINEERS 


Stimulating and creative research positions are open 
at ARMOUR RESEARCH FOUNDATION in the fields of: 


Hypersonic Flight Problems 
Hydrodynamics 

Cavitation Problems 

Gas Dynamics 


Boundary layer investigations and wake flow prob- 
lems of high speed flight are some of the programs 
and concepts you will work on. 

M.S. to Ph.D. in Aerodynamics or Mechanical Engineer- 
ing required with good theoretical background 
compressible flow, heat transfer and thermodynamics. 
We offer outstanding employment benefits including 
liberal vacation policy. If you would like to work 
in a well known research organization with some of 
the leading engineers in this field, please write to: 


E. P. Bloch 


ARMOUR RESEARCH 
FOUNDATION 


Illinois Institute of Technology 
10 West 35th Street 
Chicago 16, Illinois 
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